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ABSTRACT 


This  thesis  presents  a  high  -resolution,  computer 
simulation  model  of  aircraft  attrition  in  the  air  defense 
enviroment.  The  model  employs  extensive  pre-processing 
submodels  and  programs  in  order  to  efficiently  examine 
tactical  scenarios  and  reduce  program  execution  time.  The 
pre-processing  outputs  are  loaded  into  a  dynamic  simulation 
submodel  to  analyze  the  aircraft/air  defense  engagement 
sequence.  The  overall  simulation  model  is  modular  and  can  be 
easily  modified  to  satisfy  the  user's  particular  analysis 
objectives. 
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I .  INTRODUCTION 


A.  BACKGROUND 

In  modern  warfare,  predicting  the  attrition  of  high 
performance  aircraft  as  they  penetrate  the  hostile  area 
is  very  important.  In  general,  this  prediction  depends  on 
many  variables.  Factors  such  as  weather,  terrain,  aircraft 
performance  and  weapons  capabilities  affect  the  air  battle 
greatly.  Calculation  of  aircraft  attrition,  in  large-scale 
combat  models,  is  difficult  because  the  process  is  highly 
stochastic.  Current  modeling  has  analyzed  aircraft 
attrition  and  air  defense  effectiveness  by  aggregating 
many  factors. 

Two  large  high  resolution  combat  models,  CORDIVEM  and 
VECTOR-2,  have  used  aggregate  air  models  in  an  attempt  to 
evaluate  the  effect  of  air  defense  against  overflight 
aircraft.  The  model  designers  have  tried  to  represent 
warfare  in  a  realistic  and  accurate  way.  However  they  have 
ignored  or  oversimplified  factors  such  as  command  and 
control,  terrain  masking,  weather,  and  enemy  electronic 
countermeasures (ECM) .  The  command  and  control  issues  have 
been  oversimplified  by  assuming  that  radar  sites  and  firing 
units  act  autonomously.  Terrain  masking,  which  is  extremely 
important  for  aircraft  survivability,  has  been  totally 
ignored.  Most  terrain  models  require  many  computations, 
large  data  bases,  and  are  cumbersome.  ECM  capabilities  have 


been  disregarded  because  they  also  require  extremely  high 
resolution  modeling  to  yield  reliable  results.  High 
resolution  models  can  be  used  to  combine  these  factors  and 
add  realism  to  the  air  battle. 

Due  to  the  aggregation  methods,  sensitivity  analysis  is 
much  more  difficult.  The  model  builder  can  not  easily 
estimate  the  change  in  attrition  rate  when  critical 
parameters  such  as  aircraft  flight  level,  type  of  aircraft 
formation,  aircraft  speed  and  ECM  capabilities  are  varied. 

In  general,  the  current  air  models  have  sacrificed  detail  for 
computational  simplicity. 


B.  THE  PROBLEM 

A  high  resolution  (HIRES)  model  which  permits  realistic 
analysis  of  aircraft  attrition  must  be  very  complex  and 
contain  many  submodels  and  subprograms.  Furthermore,  the 
computer  time  required  to  execute  a  HIRES  model,  as  part  of  a 
high  resolution  wargame,  is  significant.  The  task  facing  the 
model  builder  is  to  design  a  HIRES  air  model  that  is 
efficient  and  produces  acceptable  results  in  a  timely  manner. 

Lee  and  Cho  introduced  a  framework  for  designing 
attrition  models  in  their  Naval  Postgraduate  thesis  entitled, 
A  Generalized  Approach  to  Air  Defense  Combat  Modeling, 
September  1985.  This  author  has  implemented  their  concepts 
by  designing  and  programming  a  high  resolution  aircraft 
attrition  model  which  improves  the  analyst's  capability  to 
account  for  the  uncertainties  in  the  air  attack/defense 
problem. 


C.  THESIS  OBJECTIVES 


The  objectives  of  this  thesis  are: 

1.  To  create  a  model  which  realistically  estimates  the 
attrition  of  an  aircraft  formation  during  its  route  across 
the  hostile  area. 

2.  To  program  the  model  so  that  most  calculations  are 
performed  in  an  off-line,  pre-processing  mode. 

3.  To  design  the  dynamic  simulation  segment  of  the  model 
to  take  advantage  of  the  pre-processing  output  and 
efficiently  simulate  the  stochastic  elements  of  the  air 
battle. 

4.  To  perform  sensitivity  analyses  of  the  critical 
variables  in  the  air  attack/defense  problem. 

D.  APPROACH 

The  air  model  will  be  developed  in  the  following 
sequence. 

1 .  Pre-Processing  Sub-Models 

The  deterministic  elements  of  the  air  battle  will  be 
accounted  for  in  this  phase, 
a.  Terrain 

The  terrain  module  will  be  developed  and  used  to 
calculate  a  minimum  altitude  of  detection  matrix  for  all 
radar  sites  contained  in  the  simulation  (Chapters  II  and 


b.  Detection  Geometry 

The  detection  module  will  be  developed  and  used 
to  determine  the  possibility  of  detection  based  on  line-of- 
sight  and  acquisition  radar  range  considerations  (Chapter  IV) 

c.  Aircraft  Movement 

The  aircraft  movement  module  will  be  developed  to 
simulate  the  flight  paths  of  all  the  aircraft  in  the  scenario 
(Chapter  V) . 

d.  Engagement 

The  detection,  acquisition,  interception,  and 
priorization  modules  will  be  developed  to  determine  the 
feasibility  of  the  engagement  sequence  (Chapter  VI) 

2.  Dynamic  Sub-Models 

The  stochastic  elements  of  the  air  battle  will  be 
simulated  in  this  phase.  These  sub-models  will  allow  for  the 
replication  of  scenarios  and  the  summarization  of 
probabilistic  events.  The  dynamic  portion  of  the  model  will 
use  the  outputs  of  the  pre-processing  submodels  for  each 
replication  (Chapter  VII  and  VIII). 

3 .  Scenario  Demonstration  and  Model  Execution 

In  Chapter  IX,  a  specific  scenario  will  be  described 
and  loaded  into  the  air  model.  This  scenario  will  be  used  to 
demonstrate  the  execution  of  the  computer  simulation  and 
generate  sample  outputs  of  a  typical  air  battle.  The  outputs 
will  be  statistically  summarized  and  analyzed.  Sensitivity 
analysis  will  be  performed  to  demonstrate  the  relationship 
between  aircraft  altitude  and  expected  attrition  rate. 


E.  SUMMARY 


The  thesis  will  demonstrate  that  aircraft  attrition  in  an 
air  defense  environment  can  be  efficiently  and  accuratly 
modeled  using  high  resolution  techniques  and  efficient  use  of 
off-line,  pre-processing  programs  and  routines. 


II.  TERRAIN  SIMULATION 


Terrain  simulation  must  be  developed  to  represent  the 
terrain  over  which  the  attacking  formation  flies. 

The  attacking  aircraft  a/c  uses  terrain  masking  in  order 
to  avoid  enemy  radar  detection.  Since  radars  require  line- 
of-sight,  earth  detection  can  be  restricted  by  the 
interfering  terrain  form.  In  order  to  determine  if  a  radar 
has  line-of-sight  (LOS)  to  the  a/c  in  the  air  defense 
simulation,  it  is  necessary  to  analyze  the  terrain  which  lies 
between  the  radar  and  the  a/c  formation. 

This  terrain  simulation  developed  in  this  thesis  is  based 
on  a  grid  square  method.  The  area  to  be  analyzed  is  divided 
into  uniform  squares.  Other  shape  such  as  hexagon  could 
also  be  used.  The  length  of  the  side  of  each  square  depends 
on  the  accuracy  of  the  terrain  representation  required.  The 
smaller  the  length,  the  greater  the  accuracy  of  the  terrain 
model.  If  great  accuracy  is  desired,  extensive  calculations 
must  be  anticipated.  This  chapter  will  describe  the  terrain 
model  and  the  factors  which  affect  the  terrain  simulation 
accuracy. 

A.  ACCURACY 

The  most  important  question  concerning  terrain  modeling 
is  determining  the  appropriate  degree  of  accuracy.  A  high 
resolution  terrain  representation  will  require  extensive 
computer  processing.  Furthermore  the  formation  of  the 
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attacking  a/c,  the  number  of  a/c,  and  the  flight  altitude  are 
also  factors  which  affect  the  accuracy  question  and  must  be 
considered. 


1.  Formation 


c 
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Attack  formations  vary  greatly  and  depend  on  the 


tactics  and  the  mission.  For  example,  consider  the  two  kinds 


of  formation  shown  in  Figure  2.1. 


Figure  2.2  Formation  of  Five  A/C 


Both  formations  have  three  a/c,  but  formation  one  has 


dimensions  approximately  1/3  by  1/3  NM  while  formation  two  is 


2  NM  by  2  NM.  Because  formation  one  is  smaller,  it  is  more 


flexible  and  it  can  use  terrain  masking  better  than  formation 
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two.  Thus,  if  the  attack  formations  are  of  the  first  type, 
greater  terrain  modeling  accuracy  will  be  required. 


2 .  Number  of  A/C 

The  number  of  a/c  in  a  given  formation  also 
influences  the  terrain  accuracy.  The  formation  in  Figure  2.2 
above  is  similar  to  formation  two  in  Figure  2.1?  however,  the 
number  of  aircraft  has  been  increased.  In  this  case,  the 
five  aircraft  are  flying  in  a  formation  which  is  4  by  2  NM. 
This  type  of  formation  is  less  flexible  and  must  fly  higher. 
As  a  result,  less  terrain  simulation  accuracy  is  required. 

In  particular,  the  dimensions  of  the  grid  square  have  to  be 
approximately  equal  to  the  dimensions  of  the  formation. 

3 .  Flight  Level 

The  flight  level  used  by  the  attacking  formation  is 
also  critical  to  the  determination  of  the  terrain  simulation 
accuracy.  If  the  user  wants  to  analyze  attrition  against 
extremely  low  flying  aircraft,  the  terrain  simulation  must  be 
of  very  high  resolution  in  order  to  generate  accurate 
predictions. 

B.  CALCULATION  OF  THE  GRID  SQUARE  ALTITUDES 

The  first  step  in  the  procedure  used  here  is  to  calculate 
a  representation  altitude  for  each  grid  square.  To  find  the 
altitude  of  each  grid  square  corner  a  map  is  needed.  The 
desired  accuracy  of  the  terrain  simulation  determines  the 
scale  of  this  map.  In  general,  any  military  map  with  scale 
li 50,000,  1:250,000  or  1:500,000  is  appropriate. 


In  Figure  2.3  a  military  map  with  a  1:500,000  scale  is 
shown.  Consider  the  map  as  an  example  of  a  defended  area 
divided  into  grid  squares.  The  size  of  each  square  is  equal 
to  6  NM. 

This  map  contains  all  the  information  needed  for  the 
average  altitude  calculation.  It  contains  the  scale  of  the 
ranges  and  enough  contours  to  extract  the  heights  of  the  four 
corners  of  each  grid  square  in  the  following  manner. 

If  one  contour  passes  through  a  corner  of  the  grid 
square  assign  the  altitude  of  this  contour  to  this  corner. 

If  no  contour  passes  exactly  through  a  corner,  we  interpolate 
the  altitude  between  the  two  contours  which  contain  the 
particular  corner  and  assign  that  value  to  the  corner  in 
question. 

Figure  2.4  shows  an  example.  The  grid  square  A  B  C  D 
is  used  to  calculate  the  average  heights  in  this  particular 
area.  Contours  pass  through  the  corners  A  and  D  so  the 
altitude  2,000  can  be  assigned  to  these  two  corners.  For  the 
corners  B  and  C  we  have  to  interpolate  the  heights  by  using 
the  contours  2,500  and  3,000  feet  for  point  B  and  contours 
3,000  and  3,500  feet  for  point  C.  A  good  approximation 
height  is  2,700  and  3,250  feet  for  points  B  and  C 
respectively. 


Table  1  shows  the  altitudes  of  the  grid  square  corners  of 
the  terrain  shown  in  Figure  2.3. 

To  find  these  altitudes  is  cubersome  work  but  it  must 
only  be  done  once.  Later  average  altitudes  can  be  calculated 
for  each  grid  and  a  data  matrix  can  be  stored  to  be  used  for 
many  simulation  purposes.  Table  2  shows  the  average  altitude 
of  each  grid  square.  After  the  calculations  have  been  made 


in  this  particular  model ,  each  average  altitude  represents 
the  altitude  of  a  6  by  6  NM  area. 

TABLE  2 

AVERAGE  ALTITUDES  OF  GRID  SQUARES 
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Figure  2.5,  shows  the  plot  of  terrain  made  by  a  computer 
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by  using  the  data  in  Table  2  and  the  program  in  Appendix  A 
In  this  figure  we  can  see  Monterey  Bay  and  most  of  the 
contours  which  are  shown  in  Figure  2.3  above. 


C.  ADVANTAGES  AND  DISADVANTAGES  OF  GRID  SQUARE  METHOD 


1 .  Advantages 

This  method  is  very  simple.  It  does  not  require 
difficult  mathematical  formulas  on  extensive  calculations. 

All  computations  can  be  done  off-line  and  can  be  used 
in  many  different  types  of  simulations. 

The  computer  program  which  calculates  the  average 
altitudes  -  shown  in  Appendix  A  -  can  be  used  in  any  terrain 
situation. 

2 .  Disadvantages 

The  method  is  discrete.  When  the  average  altitude 
of  the  grid  square  is  calculated,  it  is  used  to  represent  the 
altitude  for  any  point  inside  this  area.  Figure  2.6  compares 
the  actual  terrain  with  the  grid  square  approximation.  The 
terrain  is  divided  into  columns  whose  bases  are  squares  6  by 
6  NM.  By  examinimg  the  points  a  and  b  in  column  4,  we  can 
see  that  this  method  may  sometimes  overestimates  (point  b)  or 
underestimates  the  actual  terrain  height  (point  a).  This 
disadvantage  can  be  eliminated  by  decreasing  the  length  at 
the  side  of  the  square  grid,  but  with  the  rsulting  increase 
in  computation. 


Ill 


CONSTRUCTION  OP  A  MINIMUM  ALTITUDE  GRID 


A.  INTRODUCTION 

The  possibility  that  a  radar  site  will  detect  a  target 
which  penetrates  its  area  depends  critically  on  its  ability 
to  "see"  in  the  space  where  the  target  flies.  Usually, 
especially  in  high  threat  areas,  the  attacking  aircraft  use 
terrain  masking  to  avoid  radar  detection.  On  the  other  hand, 
radar  sites  try  to  improve  the  ability  to  "see"  as  low  as 
possible  by  using  more  sophisticated  equipment.  Both  radar 
sites  and  aircraft  are  concerned  about  the  minimum  altitudes 
that  a  particular  radar  can  detect  a  target.  Radar  sites  can 
use  this  information  to  find  out  their  limitations  and  to 
design  the  area  defense  in  better  way.  The  attacking  a/c  can 
use  the  minimum  altitude  grid  to  plan  their  penetration  route 
through  the  enemy  area  and  take  advantage  of  terrain  masking. 
These  routes  must  be  planned  to  avoid  radar  detection,  to 
save  fuel,  and  to  decrease  the  possibility  of  navigational 
errors. 

This  chapter  will  develop  a  model  for  calculation  the 
minimum  altitude  grid.  This  minimum  altitude  grid  will  be 
used  in  the  next  chapter  to  determine  radar  line-of-sight  and 
the  possibility  of  detection. 


B.  DESCRIPTION  OF  THE  MODEL 


S 


Consider  the  average  altitude  grid  shown  in  Figure  3.1. 


Figure  3.1  Average  Altitude  Grid 


The  line-of-sight  from  site  1  (S)  to  the  center  of  grid 


(O)  square  (4,4)  will  be  analysed. 


This  line  starts  from  the  position  of  site  1,  passes 


through  the  grid  (2,2),  (3,2),  (3,3)  and  ends  at  point 


(4,4).  A  profile  of  this  line-of-sight  is  shown  in  Figure 


3.2. 


Figure  3.2 


Profile  of  Average  Altitude 


r: 
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The  profile  can  be  adjusted  for  earth  curvature,  as  shown 
in  Figure  3.3.  Drawing  a  ray,  AB,  from  site  1  through  the 
edge  of  the  first  block  illustrates  the  terrain  masking  that 
the  first  terrain  block  creates.  Based  on  the  first  grid 
square  masking,  the  minimum  altitude  required  for  radar 
line-of-sight  to  a  target  flying  over  grid  square  (4,4)  will 
be  calculated.  The  geometry  of  this  situation,  considering 
the  grid  (2,1)  position  and  the  grid  (4,4)  position,  is  shown 
in  Figure  3.4. 


Figure  3.3  Profile  of  Average  Altitude  Adjusted 


For  Earth  Curvature 
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l.  Calculation  of  Distances  D1  and  D2 


Re  *  Radious  of  earth 


D1  ■  Distance  (AG)  between  radar  site  1  and  the 


edge  of  the  first  grid  (2,1) 


D2  "  Distance  (AF)  between  radar  site  1  and  grid 


(4.4) 


S  ■  Average  altitude  of  radar  site  1 


h  «  Average  altitude  of  (2,1)  grid 


n 


iM 


al  *  Minimum  altitude  that  the  radar  site  1  can 


"see"  due  to  the  grid  (2,1) 

The  length  of  the  arcs  AG  and  AE  can  be  calculated  by 
using  the  formula: 

1/2 

) 

(1) 

1/2 

) 

Where: 

XI,  X2  =  the  cartesian  coordinates  of  radar 
site  1. 

XI,  Y2  »  the  cartesian  coordinates  of  grid 

(2,1) 

X3,  Y3  *  The  cartesian  coordinates  of  grid 
(4.4) 

The  angles  W1  and  W2  can  be  found  by  using  the 

formula: 


ARC ( AG )  -  ( (X2-X1 )  +  ( Y2-Y1 ) 


ARC ( AE )  =  ( (X3-X1 )  +  (Y3-Y1) 


ARC (AG)  ARC ( AE ) 

W1  - - - —  (in  rands)  ,  W2  =  - - -  (2) 

Re  Re 


The  distances  D1  and  D2  can  be  now  calculated  by 
using  the  law  of  cosines: 


2 

(Dl) 


2  2 

(OA)  +  (OG)  -  2(0A) (OG)COS(Wl) 


(3) 
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Substitution  of  OA*OF»R,  ARC (AG),  ARC(AE)  and  W1  from 


formulae  (1)  and  (2)  into  formula  (3)  yields: 


2  •  21/2 
( (X2  -  XI)  +  (Y2  -  Yl)  )  *  1/2 

D1  -Re[2  -  COS( -  )] 

Re  (4) 


In  the  same  way  the  distance  D2  can  be  found. 

2 .  Calculation  of  the  altitude  al 

In  order  to  calculate  the  altitude  al,  we  also  need 
the  distance  BC  and  the  angle  W3.  These  two  can  be  found  as 
follows: 


2  2  2 

(BC)  »  (OB)  +  (OC)  -  2(0B) (OC)COS(Wl)  (5) 


From  Figure  3.4  we  can  see  that  OB  *  R+S  and  OC  = 

R+h.  By  using  this  result  and  the  formula  (2)  we  can  rewrite 
formula  (5)  as: 


222  ARC ( AG ) 

(BC)  -  (R+S)  +  (R+H)  -  2 ( R+S ) ( R+h ) COS ( - )  (6) 


Now,  we  know  the  sides  OB  and  BC  of  the  triangle  OBC, 
so  we  are  able  to  find  the  angle  W3  by  using  the  law  of 
cosines  and  formulae  (5)  and  (6)  : 


2  2  2 
(OB)  +  (BC)  -  (OC) 

W3  -  ARC  C -  ]  * 

2 (OB) (BC) 


(R+S)  [1  -  ( R+H ) COS ( W1 ) J 


(R+S)  +  (R+h) 

-  2[R+S][R+h]C0S(Wl) 

We  also 

have: 

W4  = 

(W1+W2+W3) 

Finally 

,  by  using  the 

law  of  sines  we  can  find  the 

altitude  al  as 

follows : 

SIN ( W3 ) 

SIN(W4) 

SIN(W3) (OB) 

- = - ==>  OD  - -  =»> 

OD 

OB 

SIN(W4) 

SIN(W3) (R+S) 

Ke  t  al 

SIN ( W4 ) 

SIN(W3)[R+S] 

•  n « 

ax  — 

SIN(W4) 

Thus,  given  R,  S,  h,  D1  and  D2,  the  minimum  altitude 
required  because  of  terrain  masking  by  the  first  block  can  be 
calculated. 

Formula  (9)  may  give  us  negative  altitudes.  This  can 
happen  in  the  case,  when  the  height  of  the  radar  site  is  much 
higher  than  the  height  of  the  first  block. 

Negative  minimum  altitudes  mean  that  the  radar  site 
can  "see"  even  below  the  surface  of  the  earth.  It  is  obvious 


that  in  this  case,  we  have  to  assign  the  average  height  of 
the  target  grid  (4,4),  namely  700  feet  (see  Figure  3.1). 

To  summarize,  the  complete  algorithm  for  calculating 
the  minimum  altitude  al  is: 

If  al  is  greater  or  equal  to  zero,  then  store  the 
result. 

Otherwise  store  zero  and  continue  to  calculate  a  new 
al  because  of  the  second  block 
3 .  Calculation  of  the  maximum  altitude  al 

As  one  proceeds  across  grid  squares  along  the  line- 
of-sight  path  in  Figure  3.1,  four  grid  blocks  are 
encountered.  Therefore,  minimum  altitude  calculations  would 
be  performed  four  times.  These  four  altitudes,  al  through 
a4,  are  shown  in  Figure  3.5. 

The  minimum  altitude  entry,  a44  for  grid  square  (4,4) 
in  relation  to  site  1,  is: 

a44  -  maximum  (al , a2 , a3 , a4 ) .  This  procedure  must  be 
repeated  for  every  grid  square  located  within  the  max  radar 
range  of  site  1.  Then,  the  site  1  minimum  altitude  grid 
would  be  completed. 

The  program  "MINALT"  (MINimum  ALTItude),  described  in 
Appendix  B,  calculates  the  minimum  altitudes  for  a  given 
radar  site.  A  small  output  of  this  program  is  shown  in 
Figure  3.6.  The  upper  number  of  each  square  represents  the 
average  terrain  altitude  for  this  grid  and  the  lower  number 
represents  the  minimum  altitude  that  the  radar  at  site  1, 
located  in  grid  (3,3),  can  "see".  For  example,  the  minimum 


altitude  for  grid  (5,5)  is  876  feet  instead  of  562  feet 
because  of  the  terrain  masking  that  the  grid  (4,4)  creates 


b.  The  minimum  altitudes  have  to  be  calculated  once 
for  each  radar  site.  After  that,  these  altitudes  can  be  used 


for  finding  the  line-of-sight.  They  also  can  be  used  as 
inputs  to  other  models  in  a  more  complex  simulation. 


Figure  3.6  Minimum  Altitudes  for  Site  1  (3,3)  -  Elevation  - 
2.  Disadvantages 

In  accurate  line-of-sight  determination,  many 
calculations  are  necessary.  However,  if  we  consider  that 
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these  calculations  will  be  done  once  and  can  be  performed 
in  a  pre-processing  stage,  this  algorithm  poses  no  major 
problems. 


IV.  LINE-OF-SIGHT  AND  POSSIBILITY  OF  DETECTION 

A.  INTRODUCTION 

All  radars  which  operate  in  high  frequencies  have  some 
restrictions  which  must  be  satisfied  so  that  they  will  be 
able  to  detect  a  target  flying  across  their  area.  The  two 
most  important  restrictions  are  the.  existence  of  line-of- 
sight  and  the  requirement  that  the  target  flies  inside  the 
maximum  range  of  radar  site. 

This  chapter  will  develop  a  method  to  find  out  if  there 
is  a  line-of-sight  for  a  particular  radar  site.  It  will  also 
develop  a  model  to  determinate  the  possibility  of  detection 
of  the  target  given  that  it  flies  inside  the  maximum  range  of 
the  radar  site. 

B.  LINE-OP-SIGHT  DETERMINATION 

Line-of-sight  means  that  there  is  nothing  between  the 
radar  site  and  target.  Line-of-sight  may  be  prohibited  for 
two  reasons. 

First,  terrain,  such  as  hills  and  mountains,  prohibit 
line  of  sight.  Figure  4.1  below  shows  this  case. 

The  second  reason  is  the  curvature  of  the  earth.  If  the 
distance  between  radar  site  and  target  is  lengthy  then 
line-of-sight  may  be  resticted  due  to  earth  curvature. 


Figure  4.1  Line-of-Sight  Restricted  by  Hill 
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This  case  is  shown  in  Figure  4.2.  As  an  example,  when 


the  altitude  of  the  radar  site  is  zero  and  the  distance 


between  radar  site  and  target  is  50  NM  the  target  must  fly 
approximately  2200  feet  above  ground  in  order  to  be  seen 


by  the  radar  site 


On  ths  other  hand,  if  the  altitude  of  radar  site  is  1000 


feet  and  the  target's  altitude  is  2000  feet,  the  distance 


that  the  radar  site  will  be  able  to  detect  the  target  is 
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87  NM. 
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Figure  4.2  Line-of-Sight  Restricted  by  Earth  Curvature 

In  the  program  "CONTOUR"  described  in  Appendix  C,  the 
determination  of  radar  line-of-sight  uses  the  minimum 
altitude  look  up  table  described  in  Chapter  III.  Figure  4.3 
is  the  same  as  the  Figure  3.6  in  Chapter  III.  It  contains 
the  minimum  altitudes  that  a/c  over  other  grids  must  fly  if 
radar  site  1,  located  in  grid  (3,3),  is  to  see  a  target. 

Figures  4.4,  4.5  and  4.6  are  outputs  of  the  computer 
program  "CONTOUR".  The  data  used  by  this  program  are  the 
data  shown  in  Table  2  and  Figure  3.6  in  chapter  II  and  III 
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respectively.  Assume  that  the  target  flies  at  500,  2000  and 
6000  feet  above  the  surface  of  the  earth.  This  means  that 
the  target  uses  terrain  masking.  For  example,  if  the  average 
altitude  of  the  grid  (2,2)  is  1000  feet  and  the  target  flies 
500  feet  above  this  grid,  then  the  altitude  of  the  target  is 
1500  feet.  By  comparing  this  altitude  with  the  minimum 
altitude  that  the  radar  site  can  see  in  the  grid  (2,2)  we  can 
determine  if  there  is  line-of-sight • 


Figure  4.3  Minimum  Altitudes 


In  Figures  4.4,  4.5  and  4.6  we  can  see  how  the  ability 
of  the  radar  site  increases  as  the  altitude  of  the  target 
increases.  In  these  Figures  there  is  a  flight  path  which 
starts  from  grid  (2,1)  and  ends  at  grid  (11,9).  In  Figure 
4.4,  a  target  which  is  flying  across  this  flight  path,  at  500 
feet  over  the  ground  (using  terrain  masking),  will  not  be 
detected  by  the  radar  site  located  at  the  grid  (3,5). 

However,  when  the  target  uses  the  same  path  at  2000  feet, 
above  ground  level,  it  may  be  detected  over  the  grids 
(3, 2). (3, 3), (4, 3), (5, 4), (6, 5), (7, 6), (8, 6)  and  (8,7)  (see 
Figure  4.5).  If  the  target  uses  the  same  path  at  6000 
feet,  above  ground  level,  it  will  not  be  seen  over  the  grids 
(10,3)  and  (11,3). 

This  method  can  be  used  to  create  a  look  up  table  for 
each  site  in  the  air  defense  system. 

C.  DETERMINING  THE  POSSIBILITY  OF  DETECTION 

In  the  previous  chapter  we  described  the  method  for 
determining  the  existense  of  line-of-sight  from  the  radar  to 
the  target.  In  this  section,  we  will  describe  how  to 
determine  if  there  is  a  possibility  of  detection  of  the 
target  by  the  radar  site,  given  that  there  is  line-of-sight. 
We  use  the  term  possibility  of  detection,  because,  even  if 
there  is  line-of-sight  and  the  target  flies  inside  the 
maximum  range  of  the  radar  site,  radar  may  not  be  able  to 
detect  the  target.  For  example,  the  performance  of  the  radar 
could  be  reduced  by  heavy  precipitation.  In  some  cases,  the 
target  uses  ECM  which  can  degrade  radar  performance. 
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Figure  4.4  Visible  Area  for  Target  Flying  at  500  Feet 
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Figure  4.5  Visible  Area  for  Target  Flying  at  2000  Feet 


2.  If  radar  line-of-sight  exists,  calculate  the  slant 
range  X,  from  the  site  to  the  target.  If  the  slant  range  is 
more  than  the  maximum  acquisition  range  R  of  the  radar,  no 
detection  will  occur. 

If  the  slant  range  is  less  than  or  equal  to  the  maximum 
acquisition  range,  a  detection  is  possible. 

The  calculation  of  the  slant  range  can  be  performed  as 
followst  (see  Figure  4.7) 

Consider  a  target  at  time  t  with  the  following 
parameters! 

1.  [X2(t),  Y2(t)]  is  the  target  grid  cordinates  at  time  t 

2.  a(t)  is  the  target  altitude  at  time  t 

The  lateral  range  D1  of  the  target  from  the  site  at  time 
t  is! 

2  2  1/2 

Dl(t)  -  [(XI  -  X2( t ) )  ♦  (Y1  -  Y2 ( t ) )  ] 

Where  (XI, Yl)  are  the  grid  coordinates  of  the  radar  site. 

The  angle  W  between  ray  OA  and  OB  is  Dl(t)/Re,  measured 
in  radians. 

The  slant  range  X  (AB)  can  be  calculated  from  the  law  of 
cosines! 

22  2  Dl(t) 

X  -  (Re+S)  +  (Re+a( t ) )  -  2 (Re+S ) (Re+a ( t )COS ( - ) 

Re 

Where  S  *  the  altitude  of  the  radar  site  and  Re  *  radiaus 


of  the  earth 


Figure  4.7  Calculation  of  the  Slant  Range 


The  computer  program  "CONTOUR",  explained  in  Appendix  C, 
calculates  the  possibility  of  detection.  Figures  4.4,  4.5 
and  4.6  are  outputs  of  this  program.  This  output  will  later 
be  used  by  the  target  movement  simulation. 


D.  ADVANTAGES  AND  DISADVANTAGES 


1 .  Advantages 

Determination  of  radar  line-of-sight  and  possibility 
of  detection  for  a  given  radar  site  and  target  altitude  has 
to  be  done  only  once. 

Figures  like  4.4,  4.5  and  4.6  can  be  used  by  pilots  to 
plan  their  routes.  They  can  choose  a  flight  altitude  to 
minimize  possible  detection.  This  result  is  extremely  useful 
to  the  pilots  and  it  may  be  one  of  the  most  important 
spinoffs  of  this  thesis. 

This  procedure  can  be  employed  in  air  models  as 
well  as  other  military  models  which  are  interested  in 
determining  line-of-sight. 

2.  Disadvantages 

Since  the  model  uses  the  average  grid  square  altitudes 
given  by  the  computer  program  "TERRAIN"  described  in  Chapter 
II  it  is  not  an  accurate  calculation.  The  determination  of 
line-of-sight  is  therefore  an  approximation.  However  the 
approximation  becomes  better  and  better  as  the  dimensions  of 
the  square  grid  are  reduced. 


V.  A/C  MOVEMENT  MODEL 


A.  INTRODUCTION 

The  tactics  used  by  the  aircraft  to  penetrate  the  hostile 
area  is  one  of  the  most  important  factors  affecting  the 
probability  of  survival.  The  flight  path  chosen  across  the 
enemy  country  depends  on  the  threat.  For  example,  if  there 
are  radar  sites  along  the  aircraft  route,  the  penetration 
flight  path  may  consist  of  small  connected  paths  in  order  to 
minimize  detection.  On  the  other  hand,  if  there  is  no  radar 
threat,  long  straight  paths  are  chosen  in  order  to  minimize 
navigational  errors  and  the  time  it  takes  to  get  to  the 
objective. 

B.  DESCRIPTION 

The  simulation  of  aircraft  movement  can  be  performed  in 
several  ways.  It  is  possible  to  design  the  simulation  so 
that  many  of  the  detailed  calculations  can  be  accomplished 
off-line.  Like  the  other  models  in  the  previous  chapters, 
the  aircraft  movement  model  chosen  here  is  performed  totally 
off-line.  Aircraft  flights  can  be  planned  to  take  advantage 

Aircraft  flights  can  be  planned  to  take  advantage  of 
terrain  masking.  A  typical  flight  path  was  shown  in  Figure 
5.1. 

A  flight  path  can  be  described  as  follows: 

.  Start  time  of  enemy  area  penetration. 

.  Initial  position. 


Figure  5.1  A  Typical  Flight  Path 


1 .  Start  Time  of  Enemy  Area  Penetration 

In  actual  air  battles,  the  time  attacking  formations 
choose  to  enter  the  hostile  area  depends  on  the  target  and 
the  enemy  defense  system.  For  example,  if  the  enemy  uses  a 
missile  defense  system,  the  attacking  aircraft  may  attack 
during  daylight  in  order  to  use  terrain  masking. 


In  this  model/  the  start  time  of  the  first 
penetrating  attack  formation  is  arbitrary.  For  example,  if 
the  first  a/c  enters  the  hostile  area  at  07:00  o'  clock  in 
the  morning,  then  this  is  considered  to  be  time  zero.  If  a 
a  second  formation  crosses  the  enemy  area  at  07:30  then,  its 
time  start  t  is  30  minutes.  For  computation  purposes,  the 
time  steps  from  07:00  to  07:30  are  filled  with  zeros,  so  that 
the  air  combat  model  knows  that  this  particular  formation 
enters  the  battle  at  07:30.  In  Figure  5.2  the  fixed  time  is 
07:00  so  the  start  time  of  the  formation  is  30  minutes.  In 
the  simulation  model  "MOVEMENT",  described  in  Appendix  D,  the 
start  time  is  prompted  and  it  must  be  entered  in  minutes . 

2.  Initial  Position 

The  point  that  the  attacking  aircraft  use  to  cross 
the  Forward  Edge  of  Battle  Area  (FEBA)  is  very  important.  If 
the  enemy  area  is  protected  by  radars,  then  the  initial  point 
must  be  chosen  to  minimize  detection. 

In  this  air  combat  model,  the  FEBA  is  the  line 
defining  the  area  where  attacking  aircraft  are  subject  to 
attrition.  This  line  is  chosen  by  the  model  designer  and  it 
depends  on  the  enemy  interception  capability.  For 
computation  purposes,  this  line  is  represented  by  the  Y  axis 
of  the  Cartesian  coordinate  system.  The  X  axis  is  the 
bottom  boundary  of  the  enemy  area  which  is  to  be  simulated. 

In  Figure  5.2  a  48x30  NM  enemy  area  is  shown.  The  (0,0) 
point  is  the  lower  left  point  of  the  map.  The  point  that  any 
attacking  formation  uses  to  penetrate  the  hostile  area  is 


described  by  its  relation  to  (0,0)  point  of  the  coordinate 
system.  For  example,  in  Figure  5.2,  the  start  point  for  the 
formation  is  (0,5),  which  means,  that  the  formation  crosses 
the  FEQA  5  MM  north  of  the  (0,0)  point.  In  the  same  fashion, 
any  other  point  of  the  formation  is  described  by  relating  to 
(0,0)  point.  For  example,  the  new  point  (39,2)  of  the 
formation  in  Figure  5.2  means  that  the  formation  enters  the 
simulated  battle  again,  after  attacking  the  target,  at  the 
point  39  NM  east  and  2  NM  north  of  the  (0,0)  point.  In  the 
"MOVEMENT"  model  for  each  segment  in  the  flight  path,  the 
initial  point  of  aircraft  penetration  is  prompted  and  it  must 
be  entered  in  NM. 

3.  Heading  of  the  Flight  Path 

The  heading  of  the  flight  path  depends  on  the  terrain 
masking  that  the  pilot  wishes  to  use.  Usually,  when  the 
pilot  wants  to  avoid  radar  detection,  he  tries  to  find 
routes  inside  the  hostile  area  that  gives  him  the  maximum 
protection. 

For  computation  purposes,  the  direction  is  computed 
by  measuring  the  angle  between  Y  axis  and  flight  path.  In 
Figure  5.2  the  angle  between  Y  axis  and  the  first  flight  path 
(OA)  is  60  degrees.  In  the  "MOVEMENT"  model,  for  each 
segment  of  the  flight  path,  the  direction  of  the  a/c  movement 
is  prompted  and  it  must  be  entered  in  degrees. 


4.  Distance 

Like  the  heading,  the  distance  depends  on  the  desired 
terrain  masking.  The  more  terrain  masking,  the  shorter  the 
distance  of  each  flight  path.  In  an  extreme  case,  if  the 
pilot  wants  to  take  full  advantage  of  terrain  masking,  he 
must  follow  the  terrain  form  and  continually  change  course. 
Under  these  conditions,  however,  there  is  a  good  chance  that 
the  pilot  will  lose  orientation.  On  the  other  hand,  the 
pilot  can  neglect  terrain  masking  and  fly  the  entire  mission 
in  one  straight  flight  leg.  However,  his  aircraft  is  much 
more  easily  detected  and  intercepted. 

For  computation  purposes,  the  distance  is  computed  by 
measuring  the  length  of  each  flight  path.  In  Figure  5.2,  the 
distance  of  the  first  flight  leg  is  30  NM.  In  the  "MOVEMENT" 
model,  for  each  segment  in  the  flight  path,  the  distance  is 
prompted  and  it  must  be  entered  in  NM. 

5.  Speed 

The  speed  that  penetrating  formations  use  is  a 
function  of  many  factors.  A  speed  must  be  chosen  to  allow 
for  the  pilot  to  maneuver  and  to  avoid  interceptors.  Fuel 
consumption  and  navigational  errors  must  also  be  considered. 
Additionally,  pilots  use  different  speeds  for  penetrating  the 
enemy  area  and  for  attacking  the  target.  In  Figure  5.2  the 
first  two  flight  legs  are  flown  by  using  speed  equal  to  360 
Knots  and  the  third  one,  after  the  formation  has  attacked  the 
target,  is  flown  by  using  speed  equal  to  420  Knots. 


For  computation  purposes,  the  speeds  are  chosen  by 
the  scenario  designer  for  each  flight  leg.  In  the  "MOVEMENT" 
model  the  speed  for  each  flight  path  segment  is  prompted  and 
it  must  be  entered  in  Knots. 

6.  Altitude 

Like  direction  and  distance,  the  altitude  depends  on 
the  terrain  masking.  A  high  speed,  low  level  flight  will 
minimize  radar  detection  and  interceptions.  On  the  other 
hand,  low  level  flights  use  more  fuel  which  sometimes 
restricts  the  aircraft  range  capability.  High  altitudes 
are  normally  used  during  the  return  route  in  order  to  save 
fuel.  For  computation  purposes,  the  altitudes  are  chosen  by 
the  scenario  designer  for  each  flight  leg.  In  reality,  the 
altitude  is  expected  to  be  the  most  critical  factor,  as  far 
as  the  attrition  rate  of  the  a/c  is  concerned.  In  Figure 
5.2,  the  altitudes  in  each  square  grid  are  the  average 
terrain  altitudes.  Also,  the  three  flight  legs  use  a  height 
of  1000,  500  and  10,000  respectively.  Finally  the  total 
altitude  (average  grid  altitude  +  altitude  of  a/c  above  the 
terrain)  is  shown  inside  the  parenthesis  for  each  flight  path 
and  for  each  square  grid.  This  total  altitude  will  be 
compared  in  the  next  chapter  with  the  smallest  altitude  for 
each  site  (computed  in  Chapter  4)  to  determine  the  capability 
of  the  radar  sites  to  detect  the  attacking  aircraft. 

In  the  "MOVEMENT"  model  the  altitude  is  prompted  and 
it  must  be  entered  in  feet  above  the  terrain. 


The  a/c  movement  model  is  supported  by  the  computer 
program  "MOVEMENT"  described  in  Appendix  0.  It  uses  inputs 
from  other  models  and  its  output  will  be  used  for  "FDETECT" 
and  "PRIORITY"  sub  programs  described  in  later  chapters.  The 
data  necessary  to  run  the  "MOVEMENT"  program  are  shown  in 
Figure  5.2  above.  In  Table  3.  an  output  of  this  program  is 
shown.  The  first  row  is  the  time  step,  the  second  and  third 
rows  are  the  position  of  the  a/c  (X/Y),  the  fourth  row  is  the 
altitude  over  each  position.  The  fifth  row  is  the  direction 
of  the  flight  legs. 

TABLE  3 

OUTPUT  OF  "MOVEMENT"  PROGRAM 


time  steps  : 

L 

1 

3 

A  ... 

N 

1 

x: 

1 

X 

3 

4  .  .  . 

11 

i  ! 

i 

i 

i 

Y: 

6 

8 

.  .  . 

16 

i 

i 

i 

i 

i 

altitudes: 

1000 

1125 

13 -*5 

2  8^5  .  .  . 

10000 

i 

i 

i 

i 

courses: 

60 

60 

135 

135  ... 

300 

i 

■ 

59 
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C.  ADVANTAGES  AND  DISADVANTAGES  OF  THE  A/C  MOVEMENT  MODEL 


! 

|  1 .  Advantages 

a.  The  calculation  of  the  a/c  position  (X,Y)  is 
accurate  (not  an  approximation). 

b.  The  calculation  for  each  target  route  has  to  be 
done  only  once  and  can  be  done  off-line. 

c.  Since  the  movement  model  is  independent  of  the 
air  combat  model,  it  can  be  used  for  other  air  models  and 
simulations. 

2 .  Disadvantages 

a.  The  calculation  of  the  altitude  of  each  aircraft 
over  its  position  is  done  by  summing  the  average  altitude  of 
the  terrain  and  the  altitude  of  the  a/c  above  the  terrain. 
Since,  the  average  altitudes  are  an  approximation,  then  the 
altitude  determination  of  the  of  the  attacking  aircraft  are 
similarly,  an  approximation. 
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TARGET  DETECTION.  INTERCEPTION  AND  PRIORITY  MODELS 
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In  Chapter  4,  the  minimum  altitude  that  a  radar  site  can 


detect  a  target  was  calculated  and  in  Chapter  5  the  flight 


data  (position,  altitude,  course,  etc)  were  generated  for 


each  time  step  of  the  simulation.  In  this  chapter,  the 


analysis  of  the  air  model  will  be  continued  by  combining 


Chapters  4  and  5  and  by  finding  the  points  that  a  particular 


radar  site  can  detect  a  specified  target  given  by  the 


airmodel  simulation  scenario.  Furthermore,  this  Chapter  will 


describe  the  interception  model  and  target  priority  model 


which  will  help  the  radar  site  commander  to  determine  missile 


firing  capability  against  the  target. 


A.  TARGET  DETECTION  MODEL 


The  purpose  of  the  target  detection  model  is  to  determine 


if  there  is  line-of-sight  between  a  particular  radar  site  and 


a  target  and  if  it  is  within  the  radar  site's  maximum  range. 


For  this  purpose,  the  target  detection  model  combines  the 


data  given  by  the  "MINIMUM  ALTITUDE"  model  and  the  data  given 


by  the  "A/C  MOVEMENT"  model.  The  "MINIMUM  ALTITUDE"  model. 


as  described  in  Chapter  4,  gives  the  minimum  altitudes  that 


a  radar  site  can  detect  a  target  over  its  area.  The  "A/C 


MOVEMENT"  model  as  described  in  chapter  5,  gives  the  flight 


data  of  a  target  (a/c  position,  altitude,  course,  speed) 


flying  across  the  enemy  country. 


The  algorithm  the  target  detection  model  uses  the 
following  steps: 

1.  It  takes  the  position  of  the  radar  site  and 
the  position  of  a/c  in  time  step  I. 

2.  It  checks  to  see  if  the  distance  between  the 
radar  site  and  the  target  is  less  than  the 
maximum  range  of  the  radar  site. 

3.  If  no,  set  I  =  1+1  and  go  to  step  1. 

4.  If  yes,  it  compares  the  minimum  altitude  of  the 
radar  site  with  the  minimum  altitude  of  AC  in 
grid  square  for  detection  by  radar  site  1.  If 
the  first  altitude  is  smaller  than  the  second, 
then  a  detection  is  possible.  Set  I  =  1+1  and 
go  to  step  1.  If  the  first  altitude  is  greater 
than  the  second  one,  no  detection  is  possible. 
Set  1=1+1  and  go  to  step  1. 

In  Figure  6.1,  the  minimum  altitudes  of  the  radar  site  1 
(in  parenthesis),  the  route  and  the  altitudes  of  a/c  are 
shown.  By  applying  the  above  algorithm,  a  possible  detection 
occurs  in  grid  squares  (1,7),  (2,7),  (2,6),  (3,6),  (3,5), 

(4.5)  and  (6,4).  On  the  other  hand,  there  is  no  detection  in 
grid  squares  (1,8),  (5,5),  (6,5),  (7,4)  and  (8,4)  because  the 
radar  site  can  not  see  the  target  and  in  grid  squares  (8,5), 

(9.5)  and  (9,6)  because  the  target  is  out  of  the  maximum 
range  of  the  radar  site. 
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Figure  6.1  Radar  Site-Target  Flight  Route  Relation 


In  this  way,  the  detection  state  vector  is  defined  for 
each  time  step  t  and  it  has  the  following  form: 


D(t)  -  dl(t)  ,  d2(t)  ,  d3( t ) 


dn(  t) 


Where  n  *  number  of  targets  in  the  scenario. 
di(t)  *  0,  when  the  site  cannot  detect  target  i. 


di{t;  *  1«  when  the  site  can  detect  target  i, 


The  outcomes  1  or  0  are  only  a  function  of  detection 


geometry  and  depend  only  on  the  maximum  range  and  the 


minimum  altitude  of  the  radar  site  and  the  altitude  of  the 


target. 


All  the  detection  vectors  form  the  detection  matrices. 


Each  radar  site  has  its  own  detection  matrix  for  all  the 


targets.  In  Figure  6.2  such  a  matrix  is  shown: 


TARGETS 


Figure  6.2  A  Detection  State  Matrix  for  One  Site 


For  example,  the  detection  matrix  for  radar  site  1,  in 


Figure  6.1  for  target  1,  for  time  step  1  to  9  is  (in  row  form 


to  save  space): 


^V.V.VV  xVfV  »;*A 


Time  step:  123456789 

Detection  (Y/N):  011100000 

During  the  simulation  at  time  step  k,  if  the  air  model 
needs  to  find  if  the  radar  site  1  can  detect  the  target  1, 
the  only  thing  to  be  done  is  to  check  the  element  (k,l)  in 
the  detection  state  matrix  of  radar  site  one,  to  see  if  the 
value  is  0  or  1. 

The  target  detection  model  is  supported  by  the  computer 
program  "FDETECT"  (File  DETECT)  described  in  Appendix  E.  The 
purpose  of  the  computer  program  is  to  construct  the  detection 
state  matrix  for  each  Radar  site  by  prompting  the  following 
inputs: 

a.  The  minimum  altitudes  of  the  radar  site  given  by  the 
"MINALT"  program,  described  in  Appendix  C. 

b.  The  flight  data  of  the  targets  given  by  the  program 
"MOVEMENT"  described  in  Appendix  D. 

This  computer  program  can  be  run  for  any  radar  site  and 
for  any  number  of  targets. 

B.  TARGET  INTERCEPTION  MODEL 

The  next  step,  after  the  radar  site  has  detected  the 
target,  is  to  start  the  procedure  of  intercepting  it.  In 
real  war,  the  determination  of  whether  the  missile  site  can 
intercept  the  target  or  not  is  given  automatically  by  the 
radar  computer  which  solves  the  geometry  of  the  interception 
considering  all  the  needed  paramaters.  In  simulation,  a 
program  is  required  to  solve  the  geometry  of  the 
interception.  Unfortunately,  different  kinds  of 


radar-missile  sites  use  different  methods  of  interception 
geometry.  In  Figure  6.3,  two  different  methods  are  shown: 
The  target  following  method  and  the  interception  point 
prediction. 


In  first  method,  the  firing  missile  follows  the  target 
until  it  reaches  and  kills  it.  This  method  is  used  by  old 
radar  missile  systems  or  by  missiles  which  use  infrared 
radiation  guidance  system.  The  second  method  is  used  by 


the  most  modern  radar  missile  systems  and  it  is  very 


GO 


efficient  because  it  minimizes  the  interception  time.  In 
this  method  the  missile  is  guided  to  the  point  that  the  radar 
computer  has  calculated,  so  that  the  collision  between  target 
and  missile  will  take  place.  The  method  is  simulated  in  this 
air  model  and  it  is  described  here. 

To  predict  the  point  at  which  a  launched  missile  and  the 
target  will  collide  is  not  difficult.  This  point  depends  on 
the  distance  (D)  between  target  and  radar-missile  site,  on 
the  angle  (W)  between  flight  path  and  radar-missile  site  and 
on  the  speeds  of  the  target  and  the  missile. 

In  Figure  6.4  the  analyzed  geometry  of  the  interception 
point  prediction  method  is  shown: 


TARGET  DETECTION  RADAR  NUSSU.E 

POINT  SITE  POSITION 


Figure  6.4  Interception  Geometry 


The  calculation  of  the  collision  point  is  described  as 
follows: 

Let 

W  *  The  angle  between  flight  path  of  the  target  AC 
and  the  radar  -  missile  site  position  (X,Y). 

D  ■  Distance  between  the  position  of  the  target  at 
the  time  t  and  the  missile  site  position. 

VT  *  The  speed  of  the  target  (KNOTS/h). 

VM  *  The  speed  of  the  missile  (KNOTS/h). 

T  ■  The  time  the  interception  takes. 

The  interception  point  can  be  computed  if  the  distance  AC 
is  found.  To  do  that,  the  law  of  SINES  is  used: 

SIN(W)  SIN(Wl)  SIN(W)  BC 

BC  AC  SIN(Wl)  AC 

BC  »  VM  *  T  and  AC  =  VT  *  T  (2) 

SIN ( W)  VM*T 

From  (1)  and  (2)  - - - -  ==>  SIN(Wl)  = 

SIN(Wl)  VT*T 

VT  VT 

—  *  SIN ( W)  ==>  W1  =  arcSIN  [  —  SIN(W)  ]  (3) 

VM  VM 

The  three  angles  of  one  triangle  sum  to  180  degrees,  i.e. 
W  +  W1  +  W2  =  180  degrees  ==>  W2  *  180  -  W  -  W1  (4) 

and  by  substituting  the  equation  (3)  into  (4)  yields: 


( 


W2  -  180  -  W  -  arcSIN  L  —  SINW  ] 

VM 

By  using  again  the  law  of  SINES,  the  distance  AC  can  be 
found: 

SIN(Wl)  SIN(W2)  AB  *  SIN(Wl)  D  *  SIN(Wl) 


- - - »>  AC  * - - -  ( 

AC  AB  SIN(W2)  SIN ( W2 ) 


Finally  by  comparing  the  equations  (3),  (5)  and  (6)  the 
distance  AC  is  calculated: 


AC 


VM  D  *  SIN ( W) 


VT 

VT  SIN  C  180  -  W  -  arc  (  —  *  SIN(W) )  ] 

VM 


( 


The  time,  that  the  interception  takes,  can  be  calculated 
by  dividing  the  distance  AC  by  the  speed  of  target  VT: 

D  SIN  (W) 

T  =  —  * -  (8 

VT 

VM  SIN  [  180  -  W  -  arc  (  —  *  SIN(W)  ] 

VM 

The  distance  BC,  the  distance  the  missile  travels,  is 
calculated  by  multiplying  the  missile  speed  VM  by  the  time  T 

BC  =  T  *  VM  (9 

In  summary,  if  the  parameters  D,  W,  VT  and  VM  are  known, 
we  can  find  the  distances  that  the  target  and  the  missiles 
travel  until  the  collision  occurs  and  the  time  at  which  the 
interception  takes. 


It  is  obvious  that  an  interception  is  not  always  possible 
for  the  following  reasons: 

When  the  distance  that  the  missile  has  to  travel  is 
greater  than  its  maximum  flight  range. 

When  the  angle  W  is  greater  or  equal  to  90  degrees  and 
the  VM  is  less  or  equal  to  VT.  In  this  case  the  missile  will 
never  catch  the  target. 

The  Target  Interception  Model  is  supported  by  the 
computer  work  space  "INTERCEPT",  described  in  Appendix  E, 
whose  purpose  is  to  calculate,  for  a  given  radar  missile  site 
and  targets,  the  intercept  state  vector.  For  each  time  step 
t,  the  intercept  state  vector  has  the  following  form: 

I(t)  =  (il(t),  i2(t),  i3(t ),....  iN(t))  where: 

N  »  the  number  of  targets  generated  in  the  entire 
simulation. 

I  j  ( t )  **  0,  implies  that  a  missile  fired  at  target  j 
during  this  time  step  will  have  a  predicted 
intercept  point  beyond  the  range  capability  of 
the  missile  system. 

Ij(t)  =»  1,  implies  that  a  missile  fired  at  target  j 
during  the  time  step  will  have  a  predicted 
intercept  point  within  the  capability  of  the 
missile  system  for  one  of  the  reasons  listed 
above. 

All  the  intercept  state  vectors  form  the  intercept  state 
matrix,  whose  form  is  shown  in  Figure  6.5.  The  interception 


state  matrix  and  the  time  that  the  possible  interceptions 
require  will  be  used  later  by  the  air  model  to  determine 
whether  or  not  an  interception  is  successful.  It  will  also 
be  used  to  schedule  the  interception  as  future  events. 


Figure  6.5  An  Interception  State  Matrix  for  One  Site 
C.  TARGET  PRIORITY  MODEL 

Modern  radar  sites  can  track  and  fire  on  more  than  one 
target  and  almost  all  radar  sites  (old  and  modern)  can  detect 
more  than  one  target  simultaneously.  In  such  cases,  each 
radar  site  has  to  know  the  priority  of  target  engagement. 

To  simulate  the  priority  of  the  target,  either  simple  or 
complex  techniques  can  be  used.  In  this  air  model,  the 
following  technique  is  employed.  The  priority  of  each  target 


is  determined  from  its  flight  route  for  each  time  step. 
Analytically,  the  target  priority  model  uses  the  following 
rules: 

Rule  1:  Two  or  more  attacking  targets 

When  two  or  more  targets  are  penetrating  the 
defended  area  at  the  same  time,  the  aircraft  which  heads  more 
directly  towards  the  defenderr's  critical  assets,  has  highest 
priority.  Figure  6.8  shows  this  case.  The  difference 
between  the  critical  heading  of  the  course  and  the  heading 
of  target  1  is  smaller  than  the  difference  between  heading 
critical  course  and  heading  course  of  target  2.  Therefore 
target  1  has  priority  one  and  target  2  has  priority  two. 


Figure  6.6  Two  Attacking  Targets 


Rule  2:  Attacking  and  Returning  Aircraft 


When  the  priority  model  -simulates  two  target8  one  of 
which  ia  attacking  and  the  other  ia  returning  to  its  own 
territory,  it  assigns  a  higher  priority  to  the  first  target 
than  to  the  second.  The  reason  is  that  the  attacking 
targets  are  more  dangerous  than  the  returning  ones.  This 
case  is  shown  in  Figure  6.7.  Target  1  has  higher  priority 
than  target  2. 


Figure  6.7  Attacking  and  Returning  Targets 

Rule  3 i  Two  or  more  Returning  targets 

In  the  case  of  two  or  more  returning  targets,  the 
target  which  heads  more  directly  away  from  the  defender's 
critical  assets  has  the  higher  priority.  The  reason  is  that 
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Pi(t)  ■  1  when  the  target  i  has  priority  1  etc. 

All  the  priority  state  vectors  form  the  priority  state 
matrix  whose  form,  for  5  targets,  is  shown  in  Figure  6.9. 
For  example  in  time  step  1,  in  Figure  6.9,  target  4  has 
priority  one,  target  2  has  priority  two,  target  1  has 
priority  three  and  so  on. 


Figure  6.9  Priority  State  Matrix 


The  target  priority  model  is  supported  by  the  "PRIORITY* 
computer  work  space  described  in  Appendix  E.  The  inputs  of 
this  program  are  the  critical  course  of  the  ground  target 
that  the  attacking  aircraft  are  going  to  attack  and  the 
flight  data  of  each  attacking  formation.  The  output  of  the 


VII.  ADD I T I ON AL  MODELS 


Chapter  seven  describes  the  following  models: 

A.  IDENTIFICATION  FRIEND  OR  FOE  MODEL  (IFF) 

B.  ELECTRONIC  WARFARE  (ECM)  MODEL 

C.  COMMAND  AND  CONTROL  MODEL 

D.  RADAR  PERFORMANCE  MODEL 

E.  PROBABILITY  OF  KILL  MODEL 

A.  IDENTIFICATION  FRIEND  OR  FOE  (IFF)  MODEL 

1 .  Introduction 

Another  problem  that  a  radar  site  commander  has 
when  a  target  is  detected  by  the  radar  site,  is  to 
determine  if  a  target  is  friend  or  enemy.  Since  this 
identification  is  impossible  to  obtain  by  optical  means, 
both  radar  sites  and  aircraft  are  equipped  with  electronic 
devices  which  cooperate  to  give  the  identification  of  the 
aircraft. 

2 .  Description 

The  IFF  process,  during  operations,  is  very 
complicated.  There  are  many  decisions  which  must  be  made,  as 
far  as  the  identification  of  a  target  is  concerned.  For 
example,  what  must  be  decided  if  the  IFF  of  radar  site  is  not 
working?  Will  the  target  be  regarded  as  a  friend  or  a  foe? 

In  this  air  model,  the  IFF  process  is  illustrated  in 


Figure  7.1. 
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Figure  7.1  Logic  Diagram  of  IFF  Model 

The  above  diagram  is  explained  ass 
a«  Zs  IFF  Applied?  Sometimes,  radar  sites  do  not 
apply  IFF  and  they  consider  all  targets  to  be  enemies.  That 


may  happen,  for  example,  during  a  certain  period  of  time, 
when  there  are  no  friendly  aircraft  scheduled  to  fly  within 


range  of  the  particular  radar  site.  In  this  case  radar 
operators  gain  time  during  the  interceptions,  but  they  assume 
that  all  detected  targets  are  enemies  and  they  can  be  killed. 
On  the  other  hand,  if  IFF  is  applied  the  next  question  must 
be  answered. 

b.  Is  IFF  working?  Like  other  electronic  devices, 
the  IFF  equipment  is  not  always  working.  This  may  happen 
because  either  the  IFF  equipment  of  the  aircraft  or  the  IFF 
equipment  of  the  radar  site  is  not  responding  properly. 

If  IFF  is  not  working,  the  user  must  decide  whether  to  treat 
the  target  as  friend  or  as  foe.  If  a  target  is  friend  the 
radar  site  stops  the  interception  process  and  searches  for 
another  target.  If  the  target  is  enemy,  the  interception 
process  continues. 

c.  Is  IFF  correct?  Unfortunately,  IFF  systems  do 
not  always  give  correct  results.  There  is  a  possibility  that 
even  if  the  target  is  an  enemy,  the  IFF  system  will  indicate 
that  the  target  a  is  friend.  On  the  other  hand,  there  is  a 
possibility  that,  even  if  the  target  is  a  friend,  the  IFF 
system  will  indicate  that  the  target  is  an  enemy.  If  the 
first  case  happens,  there  is  no  big  problem,  just  one  enemy 
aircraft  will  pass  without  being  intercepted.  If  the  second 
case  happens,  friends  will  kill  friends  and  a  fratricide  may 
occur. 
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The  air  model  prompts  the  following  inputs  for  IFF 

model : 

1.  If  IFF  is  applied. 

2.  What  is  the  probability  that  IFF  works. 

3.  What  is  the  identification  of  the  target  if  the 
IFF  does  not  work  (friend  of  foe). 

4.  What  is  the  probability  that  the  IFF  gives  a 
correct  answer  (given  that  the  IFF  works). 

B.  ECM  MODEL 

1 .  Introduction 

As  the  technology  of  the  antiaircraft  weapons 
(missiles,  radar  guided  AAA  etc)  has  been  improved,  the 
aircraft  self  protection  means  against  the  antiaircraft 
weapons  have  also  been  improved.  Today,  almost  all  airplanes 
use  electronic  devices  for  antimissile  purposes  such  as 
electronic  countermeasure  means  (ECM),  chaffs,  flares  etc. 
Since  this  airmodel  deals  only  with  radar  sites  which  use 
guided  missiles  to  destroy  the  aircraft,  only  the  electronic 
countermeasure  means  (ECM)  will  be  simulated. 

2 .  Description 

The  electronic  countermeasure  means  are  electronic 
devices  which  first  recognize  that  the  formation  has  been 
detected  and  acquired  by  the  enemy  radar  sites.  Then  they 
try,  by  transmitting  electronic  energy,  to  brake  the 
acquisition  and  if  possible  to  "blind"  the  enemy  radar  for  as 


long  as  the  aircraft  formation  flies  within  its  range. 
Figure  7.2  shows  how  a  radar  site  detects  and  acquires  an 
aircraft  formation. 


Figure  7.3  Aircraft  Formation  Using  ECM 


a.  ECM  Parameters 

The  success  of  the  ECM  employed  by  the  aircraft 
formation  depends  on  many  factors  but  the  most  important  are 
the  power  of  the  ECM  device  (P),  the  transmit  angle  (W)  and 
the  distance  between  aircraft  formation  and  radar  site  (R). 
This  is  shown  in  Figure  7.4. 

(1)  Power.  Like  every  electronic  device,  the 
ability  of  ECM  depends  on  how  much  electronic  energy  is 
transmitted.  The  more  energy  there  is,  the  more  efficient 
the  ECM*  Unfortunately,  large  power  ECM  devices  are  very 
heavy  and  only  special  electronic  aircrafts  can  carry  them. 
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Figure  7.4  Parameters  of  ECM  Device 

Most  fighter  aircrafts  use  small  or  moderate 
power  ECM  and  improve  the  ECM  efficiency  by  narrowing  the 
used  transmitting  angle  w. 

(2)  Angle  of  transmit inq  energy.  As  radar 
sites,  aircraft  which  apply  ECM  uses  the  angle  of 
transmitting  energy,  in  order  to  improve  their  ability  to 
blind  or  break  the  acquisition  of  the  radar  site.  By 
decreasing  the  angle,  less  electronic  energy  is  being  used  to 
achieve  the  same  purpose. 
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Figure  7.5  ECM  Efficiency  -  Distance 

(3)  Distance  between  ECM  and  Radar  Site.  The 
distance  is  very  important.  The  efficiency  of  ECM  is 
inversely  proportional  to  the  distance  square.  How  ECM 
efficiency  changes  by  decreasing  the  distance  between 
aircraft  formation  and  radar  site  is  shown  in  Figure  7.5. 


b.  ECM  Function 


Electronic  countermeasure  means  can  break  the 
acquisition  (lock-on)  or  increase  the  required  time  for 
identification,  track,  fire  and  target  kill  assessment 
events.  To  simulate  the  whole  ECM  process,  two  functions  are 
needed.  One  function  which  increases  the  required  time  for 
the  interception  events  and  another  function  which  changes 
the  probability  of  breaking  the  acquisition  by  the  radar  site 
during  the  interception  process. 

(1)  Interception  Events  Function.  The  function 
which  increases  the  required  time  for  identification,  track, 
fire  and  target  kill  assessment  events  is  assumed  linear  and 
has  the  following  form: 

Time  Events  Matrix  **  (Time  Events  Matrix)  * 

Increase  Factor  ( 1 ) 

where  time  events  matrix  is  a  (N,3)  matrix  shown  in  Figure 
7.6,  and  the  increase  factor  is  a  factor  that  the  events 
matrix  must  be  multiplied  with,  in  order  that  all  time  events 
be  increased.  For  example,  if  the  increase  factor  is  equal 
to  1.5  (this  value  is  used  by  the  air  model)  the  time  events 
matrix  will  be  increased  by  50%. 


aircraft  formation  is  acquired  by  a  radar  site  and  employs 


ECM,  there  is  a  high  probability  that  the  radar  site 


acquisition  can  be  broken.  The  breaking  acquisition  function 


simulates  this  probability  by  using  the  mathematics  model 


shown  in  Figure  7.7. 
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Figure  7.7  Breaking  Acquisition  Model 


This  model  assumes  that  the  probability  of 
breaking  the  radar  acquisition  (Pr)  is  zero  when  the  distance 
(R)  between  aircraft  formation  and  radar  site  is  greater  than 
Rmax,  increases  as  the  distance  R  decreases,  reaches  a 
maximum  value  (Pmax)  when  the  distance  R  is  equal  to  R2, 
remains  constant  for  the  distance  between  R2  and  R1  and  then 
decreases  when  the  distance  R  is  less  than  Rl.  Rl,  R2,  Rmax 
and  Pmax  are  parameters  of  the  breaking  acquisition  function 
specified  by  the  user.  The  equations  which  gives  the 
probability  of  acquisition  breaking  are  as  follows: 

(Pm)x(R) 

Pr  - -  IF  0  <  R  <  Rl 

Rl 
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Pr 


Pm 


IF 


R1  <  R  <  R2 


Pr 


(PKmax) *(R2-R) 
R2-RM 


IF 


R2  <  R  <  RM 


where*  R  is  the  distance  between  aircraft  formation  and  radar 
site. 

Rl,  R2,  Rmax,  Pmax  are  parameters  of  breaking 
aquisition  function  (Figure  7.7). 

Each  aircraft  formation  has  its  own  Rl,  R2,  Rmax 
and  Pmax  and  these  data  are  given  to  the  air  model  in  matrix 
form  as  follows: 


Rl 

R2 

Rm 

Pmax 

a/c 

formation 

1 

5 

30 

80 

0.7 

a/c 

formation 

2 

5 

40 

70 

0.5 

a/c 

formation 

• 

3 

a 

0 

• 

0 

• 

0 

• 

a/c 

formation 

N 

5 

20 

50 

0.5 

The 

first  ■ 

column  is 

the 

Rl  distance 

for  a/c 

formations. 

The 

second 

column  is 

the 

R2  distance 

for  a/c 

formations. 

The 

third 

column  is 

the 

Rmax 

distance  for  a/ 

c  formations 

The 

fourth 

column  is 

the 

Pmax 

for  a 

/c  formations. 

The  fifth  column  is  the  probability  than  ECM  works. 

The  third  row  is  zero  because  the  .third  a/c  formation  does 
not  employ  ECM. 


The  logic  diagram  for  the  ECM  process  is  shown 


in  Figure  7.8. 
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Figure  7.0  ECM  Process  Logic  Diagram 


The  algorithm  of  the  ECM  process  follows  the 


steps: 


Step  1*  Does  the  aircraft  formation  apply  ECM? 

If  NO:  No  changes  to  interception  time  events 
If  YES:  Go  to  step  2. 


Step  2:  Does  ECM  work? 


If  NO:  No  changes  to  interception  time  events. 

If  YES:  Go  to  step  3. 

Step  3:  Is  radar  acquisition  broken? 

If  NO:  Identification,  track,  fire  and  target 
assessment  time  events  will  increase 
according  to  equation  (1). 

If  YES:  Further  intercept  process  is  cancelled 

and  target  goes  again  for  identification. 

C.  COMMAND  AND  CONTROL  MODEL 

1 .  Introduction 

When  two  or  more  radar  sites  cooperate  to  defend  an 
area,  a  centralized  command  and  control  organization  must 
exist  to  control  all  the  radar  sites.  For  example,  if  there 
are  more  than  two  targets  flying  over  this  area,  the 
centralized  organization  will  decide  what  target  radar  site 
one  will  engage,  what  target  radar  site  two  will  engage  and 
so  on.  In  this  way,  more  targets  will  be  destroyed  because 
each  radar  site  will  engage  only  its  own  targets. 

2 .  Description 

Since  a  complex  command  and  control  model  is  very 
difficult  to  simulate,  the  airmodel  uses  the  following  simple 
algorithm:  No  target  will  be  engaged  by  more  than  one  radar 

site.  That  means  that  when  a  particular  target  is  engaged 
by  a  radar  site,  no  other  radar  site  is  allowed  to  engage  it 
until  the  target  is  released  (not  destroyed)  by  the  original 
radar  site.  By  using  this  algorithm  there  is  not  possibility 
for  two  radar  sites  to  destroy  the  same  target.  This  rule 


can  be  revised  and  the  aircraft  formation  can  be  engaged  by 
more  than  one  radar  site  if  the  number  of  aircraft  in  the 
formation  is  greater  than  a  certain  number  N.  This  number 
is  prompted  by  the  air  model  and  it  must  be  entered  by  the 
user.  For  example,  if  N  ■  3  and  the  aircraft  formation  has 
five  aircraft,  then  this  formation  can  be  engaged  by  more 
than  one  radar  site  until  only  3  aircraft  remain  in  this 
formation. 

Figure  7.9  shows  the  logic  diagram  of  command  and 
control  model. 


Figure  7.9  Logic  Diagram  of  Command  and  Control  Model 


The  algorithm  of  the  command  and  control  model 
follows  these  steps: 

Step  1:  Does  command  and  control  model  exist? 

If  NO:  Every  radar  site  can  engage  this 
target. 

If  YES:  Go  to  step  2. 

Step  2:  Is  the  number  of  aircraft  less  than  N  (user 

input)? 

If  NO:  Every  radar  site  can  engage  this 
target. 

If  YES:  Only  one  radar  site  is  allowed  to 
engage  this  target. 

D.  RADAR  PERFORMANCE  MODEL 

1 .  Introduction 

The  performance  of  the  radar  sites  depends  on  many 
factors  such  as  age  of  radar  site,  experience  of  operators, 
weather  conditions,  and  enemy  ECM.  In  addition,  there  is  a 
possibility  that  the  radar  site  will  be  partially  or  totally 
destroyed  during  the  operations. 

2 .  Description 

The  air  model  uses  a  simple  radar  performance  model. 
It  aggregates  all  factors  which  effect  the  performance  of 
radar  site  to  one  factor  called  "performance”.  This  factor 
is  prompted  by  the  airmodel  and  it  takes  one  of  the  two 
values  that  follow: 
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Figure  7.10  Logic  Diagram  of  Radar  Site  s 
Performance  Model 


If  radar  sites  operate  in  full  performance 
(performance  value  is  1)  there  is  no  change  to  airmodel.  If 
radar  sites  work  in  limited  performance  (performance  value  is 
0)  the  identification,  track,  fire  and  target  assessment  time 
are  increased  by  a  factor  input  by  the  user.  The  probability 


that  a  radar  site  acquisition  will  be  broken  due  to  ECM 
can  also  be  increased  by  a  factor  determined  by  the  user, 
the  radar  site  operates  in  limited  performance  situation, 
time  to  identify,  track,  fire  and  assess  a  target  kill 
increases  and  the  chance  that  a  target  acquisition  will  be 
broken  is  also  more  likely. 

E.  PROBABILITY  OF  KILL  MODEL 
1 .  Introduction 

Calculation  of  probability  of  killing  (PK)  an 
aircraft  by  using  a  guided  missile  is  very  important,  but 
unfortunately,  extremely  difficult  to  calculate. 
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Figure  7.11  Probability  of  Kill  Model 


The  problem  is  that  the  probability  of  kill  depends 


on  so  many  factors,  known  and  unknown,  that  it  is  essentially 
impossible  to  create  a  relatively  simple  model  to  estimate  it 
using  all  the  factors. 

2.  Description 

The  airmodel  uses  a  simple,  but  quite  reasonable 
model  to  estimate  the  PK.  It  assumes  that  the  probability 
changes  Linearly  as  the  range  between  the  a/c  and  the  radar 
site  changes.  Figure  7.11  shows  how  the  PK  changes 
relatively  to  R. 

The  PK  is  zero  when  the  distance  (R)  between  aircraft 
formation  and  radar  site  is  greater  than  the  maximum  range 
(Rmax)  of  missile  that  this  particular  radar  site  uses.  The 
propability  of  kill  increases  as  R  decreases  (aircraft 
approaches  the  radar  site)  reaching  the  maximum  value  (PKmax) 
at  R2,  it  stays  constant  for  the  distances  between  R2  and  R1 
and  then  decreases  when  the  distance  R  is  less  than  the 
distance  Rl.  R1 ,  R2,  Rmax  and  PKmax  are  values  prompted  by 
the  airmodel. 

The  equation  used  to  calculate  the  PK,  are  the  same 
as  in  the  ECM  model,  iss 

(PK)  x  (R) 

PK  - -  IF  0  <  R  <R1 

Rl 


PK  -  PKmax  IF 

( PKmax )x(R2-R) 

PK  > -  IF 

R2-Rmax 


Rl  <  R  <  R2 


R2  <  R  <Rmax 
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where:  R  *  the  distance  between  aircraft  target  and  radar 
site. 

Rl«  R2,  Rmax,  Pmax  are  parameters  of  the  probability 
of  kill  function  and  must  be  entered  in  the  air  model  by  the 


each 

radar  site 

in 

matrix 

form  as 

follows: 

R1 

R2 

Rmax 

PKmax 

Radar 

Site  1 

5 

15 

20 

0.50 

Radar 

• 

Site  2 

7 

• 

35 

• 

40 

• 

0.55 

• 

Radar 

Site  N 

13 

50 

55 

0.40 

The  first  column  is  the  R1  distance  for  each  radar  site. 
The  second  column  is  the  R2  distance  for  each  radar  site. 
The  third  column  is  the  Rmax  distance  for  each  radar  site. 
The  fourth  column  is  the  PKmax  for  each  radar  site. 


VIII.  COMPUTER  PROGRAM  "AIRMODEL 


Chapter  VIII  describes  the  "AIRMODEL”  computer  work 
•pace.  It  defines  the  functions,  inputs,  and  outputs  of  the 
program. 

A.  INTRODUCTION 

In  the  previous  chapters,  the  pre-processing  models 
were  described  in  detail.  All  these  models  were 
deterministic  because  the  outputs  are  the  same,  given  the 
same  inputs.  These  models  calculate  all  the  data  that  the 
"AIRMODEL"  work  space  needs  to  run.  The  purpose  of 
"AIRMODEL"  is  to  calculate  the  attrition  rate  of  the 
attacking  aircraft  formation  by  using  the  data  given  by  the 
models  described  in  Chapter  II  through  VII.  This  program 
simulates  the  aircraft  formation  engagement  sequence  by  the 
radar  sites.  An  air  defense  radar  site  follows  a  well 
defined  engagement  sequence  as  it  attempts  to  destroy  enemy 
aircraft.  The  steps  of  the  sequence  ares 

Step  Is  Detection.  The  system/operator  observes  the 

presense  of  an  aircraft  formation  on  the  radar 
screen  in  the  assigned  airspace. 

Step  2s  Identification.  The  operator  identifies  the 
target  as  friend  or  foe. 

Step  3s  Track.  The  site  tracking  radars  obtain  radar 


lock  on  the  target 


Step  4*  Fire.  The  site  launches  a  missile  attempting 
to  destroy  the  target. 

Step  5s  Intercept.  The  missile  arrives  at  a  predicted 

point  in  space.  This  point  is  assumed  to  be  the 


location  of  the  collision  of  the  missile  and  the 
target. 

Step  6:  Aircraft  Kill  Assessment.  The  system/ operator 

decides  if  target  destruction  has  been  achieved. 
These  steps  are  illustrated  in  Figure  8.1. 


Figure  8.1  The  Engagement  Sequence 


If  the  target  leaves  radar  surveillance  (either 
acquisition  or  tracking),  the  sequence  is  terminated 


the  target  returns  to  the  site  detection  envelop*  it  is 
considered  to  be  a  new  acquisition. 


The  time  delays  associated  with  each  step  of  the  sequence 
are  defined  as  follows i 

Td  »  time  between  maximum  possible  detection  and  actual 
detection. 

Tid  ■  time  between  actual  detection  and  identification. 

Tt  ■  time  between  identification  and  track. 

Tf  •  time  between  track  and  fire. 

Ti  ■  time  between  fire  and  interception. 

Ta  ■  time  between  intercept  and  assessment  of  the 
intercept. 

The  flow  chart  of  the  "AIRMODEL"  is  shown  in  Figure  8.2. 
B.  DESCRIPTION 

The  work  space  "AIRMODEL"  consists  of  the  following 
computer  functions! 

1.  MAIN 

2.  SITE 

3.  IDIF 

4.  TRACK 

5.  FIRE 

6.  INTER 

7.  ASSES 

8.  PKILL 

9.  RANDOM 
1 .  Main 

a.  Header 

The  header  of  the  "MAIN"  program  is  ""MAIN"  and 
the  computer  program  is  listed  in  Appendix  F. 
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b.  Purpose 

The  purpose  of  the  main  computer  function  is  to 
prompt  all  the  inputs  from  the  user  and  to  call  the  site 
function  in  order  to  calculate  the  attrition  rate  of  the 
aircraft  formations. 

c.  Input 

The  "MAIN"  computer  program  prompts  the  following 

inputs: 

1.  RLi  The  seed  for  the  random  number  function. 
This  seed  is  needed  to  generate  random  variable  inputs. 

If  the  main  program  runs  the  same  data  twice,  the  same 
sequence  of  the  output  will  be  produced  unless  the  seed  is 
changed . 

2.  Times:  The  number  of  times  that  the  program 

will  run. 

3.  NSITES:  The  number  of  radar  sites  in  the 


particular  scenario. 

4.  DS2:  The  detection  matrix  of  each  radar 

site,  given  by  the  "FDETECT"  computer  work  space.  The  "MAIN" 
program  will  prompt  the  entering  of  the  detection  matrix 
based  on  the  number  of  radar  sites. 

5.  IS2:  The  interception  matrix  of  each  radar 
site,  given  by  the  "INTERCEPT"  computer  work  space.  The 
"MAIN"  program  will  prompt  the  entering  of  the  interception 
matrix  based  on  the  number  of  radar  sites. 

6.  TGT:  The  number  of  aircraft  that  each 


formation  has.  If  the  number  of  aircraft  is  zero  during  the 


simulation,  the  aircraft  formation  is  deleted  from  the 
simulation  process.  The  TGT  input  is  a  vector  whose  number 
of  elements  must  be  equal  to  the  number  of  aircraft 
formations  (each  element  equal  to  the  number  of  a/c). 

7.  CAP 2:  The  number  of  aircraft  that  each 
radar  site  can  engage.  If  a  particular  radar  site  has 
engaged  its  maximum  number  of  aircraft,  it  is  considered  to 
be  satutated  and  is  not  allowed  to  engage  more  targets.  The 
CAP 2  input  is  a  vector  whose  number  of  elements  must  equal 
the  number  of  radar  sites  (each  element  for  each  radar  site). 

3.  MCAP2:  The  number  of  missiles  that  each 
radar  site  possesses.  If  a  particular  site  has  launched  all 
its  missiles,  then  it  is  deleted  from  the  simulation  process. 
The  input  MCAP2  is  a  vector  whose  number  of  elements  must  be 
equal  to  the  number  of  radar  sites  (each  element  for  each 
radar  site). 

9.  PKMs  The  probability  of  kill  matrix,  given 
by  the  "PROBABILITY  OF  KILL"  model. 

10.  IFF:  The  IFF  input  takes  one  of  two  values  as 

follows: 

1  ■  IFF  is  applied  in  the  simulation  process. 

0  =  IFF  is  not  applied. 

If  IFF  is  not  applied,  the  "MAIN”  program  assumes 
that  a  detected  target  is  always  friend. 

11.  IFFWORK:  The  probability  that  the  IFF  works 

12.  IFFFOE:  The  IFFFOE  input  takes  one  of  two 
values  as  follows: 
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1  ■  a  target  is  friend  if  IFF  does  not  work. 

0  ■  a  target  is  foe  if  IFF  does  not  work. 

13.  IFFCOR:  The  probability  that  IFF  is  correct 
given  that  it  works* 

14.  ECMi  The  electronic  countermeasure  matrix, 
given  by  the  NECMM  model. 

15.  CMDs  The  input  CMD  takes  one  of  two  values  as 

follows: 

1  ■  command  and  control  model  is  applied. 

0  *  command  and  control  model  in  not  applied. 

16.  PR:  The  priority  matrix,  given  by  "THE 
PRIORITY"  model. 

17.  AUTO:  The  condition  which  determines  whether 
the  radar  sites  operate  autonomusly,  given  that  command  and 
model  is  applied. 

18.  WTH:  The  input  WTH  takes  one  of  two  values  as 

follows: 

1  *  radar  site  operation  is  normal. 

0  =  radar  site  operation  is  degraded. 

19.  LBR:  The  factor  describing  the  increase  in 
acquisition  time,  given  that  the  operation  of  the  radar  sites 
is  degraded. 

20.  ETSF:  The  factor  describing  the  increase  in 
firing  time  given  that  the  operation  of  radar  sites  is 
degraded. 


21*  TSTEPls  The  (N,4)  event  time  steps  matrix. 

The  columns  of  the  TSTEP1  matrix  are  the  required  time  for 
identification,  track,  fire  and  target  damage  assessment. 

The  number  of  rows  (N)  is  the  number  of  radar  sites, 
d.  Output 

The  "MAIN"  program  gives  the  following  outputs: 

1.  Summary  of  the  User  Scenario. 

2.  Attrition  Matrix  in  details.  The  form  of  the 
attrition  matrix  is  as  follows: 

MISSILE  SITE:  11221221 

TIME  STEP:  10  32  45  102  105  105  14  250 

NUMBER  OF 

TARGETS  KILLED:  13231123 

The  first  row  consists  of  the  number  of  the  radar 
sites  which  have  achieved  a  target  kill.  The  second  row 
documents  the  time  steps.  The  third  row  consists  of  the 
numbers  of  target  kills.  For  example,  the  first  column  means 
that  the  first  radar  site  destroyed  one  aircraft  from 
formation  one  in  time  step  ten. 

3.  Remaining  missiles  for  each  site. 

4.  The  number  of  targets  which  are  destroyed 
but  not  assessed  by  the  radar  site.  In  some  cases,  it  is 
possible  to  intercept  a  target  without  receiving  clear 
indications  of  destruction. 

5.  The  number  of  target  kills  per  aircraft 


formation 


2.  Site 


a.  Header 

The  header  of  “SITE"  function  is  “SS  SITE  I"  and 
the  computer  program  is  listed  in  Appendix  F. 

b.  Purpose 

The  function  "SITE"  combines  the  database  and  the 
other  functions  used  by  the  work  space  "AIRMODEL"  to 
calculate  the  attrition  rate  of  the  attacking  aircraft 
formations  against  ground  targets.  This  function  is  called 
by  the  main  program  once  for  each  radar  site  for  each  time 
step. 

c.  Input 

1.  The  left  argument  (SS)  is  the  number  of  the 

radar  site. 

2.  The  right  argument  (I)  is  the  time  step. 

d.  Output 

The  function  SITE  does  not  provide  any  output, 
instead  it  calls  the  IDIF,  TRACK,  FIRE,  INTER,  ASSES,  RANDOM 
and  PKILL  functions  in  order  to  schedule  future  events  and 
calculate  the  attrition  rate. 

3.  IDIF 

a.  Header 

The  header  of  the  "IDIF"  function  is 
"ID  <-  ISIT  IDIF  R";  the  computer  program  is  listed  in 
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b.  Purpose 

This  function  is  called  by  the  program  "SITE"  to 
perform  the  identification  of  the  targets. 

c.  Input 

1.  The  left  argument  (ISITE)  is  a  vector  whose 
first  element  is  the  number  of  the  radar  site  and  whose 
second  is  the  time  step. 

2.  The  right  argument  (R)  is  the  number  of  the 

target. 

d.  Output 

The  ouput  is  the  future  identification  event. 

4.  Track 

a.  Header 

The  header  of  the  function  "TRACK"  is 
"TT  <-  ISITE  TRACK  T";  the  computer  program  is  listed  in 
Appendix  F. 

b.  Purpose 

The  function  TRACK  is  used  by  the  SITE  program  to 
perform  the  track  target  event. 

c.  Input 

1.  The  left  argument  (ISITE)  is  a  vector,  whose 
first  element  is  the  radar  site  number  and  whose  second  is 
the  time  step. 

2.  The  right  argument  (T)  is  the  number  of  the 

target. 

d .  Output 

The  output  is  the  future  track  event. 


a.  Header 


The  header  of  the  function  "FIRE"  is 
"FR  <-  ISITE  FIRE  F";  the  computer  program  is  listed  in 
Appendix  F. 

b.  Purpose 

The  function  FIRE  is  called  by  the  SITE  program 
to  perform  the  missile  fire  event. 

c.  Input 

1.  The  left  argument  (ISITE)  is  a  vector,  whose 
first  element  is  the  radar  site  number  and  whose  second  is 
the  time  step. 


2.  The  right  argument  (F)  is  the  number  of  the 


target. 


d.  Output 


The  output  is  the  future  missile  Fire  event. 


6.  INTER 


a.  Header 


The  header  of  the  function  "INTER"  is 
"INT  <-  ISITE  INTER  IN";  the  computer  program  is  listed  in 
Appendix  F. 

b.  Purpose 

This  function  is  called  by  the  SITE  program  to 
perform  the  interception  of  the  target. 


ms 


c.  Input 

1.  The  left  argument  (ISXTE)  is  a  vector  whose 
first  argument  is  the  radar  site  number  and  whose  second  is 
the  time  step. 

2.  The  right  argument  (IN)  is  the  number  of  the 

target. 

d .  Output 

1.  The  time  required  for  interception. 

2.  The  distance  that  the  target  flies  until  it 
is  intercepted. 

3.  The  distance  that  the  interceptor  flies  until 
the  interception  is  performed. 

4.  The  future  target  interception  event. 

7.  ASSESS 

a.  Header 

The  header  of  the  function  "ASSES"  is 
"AM  <-  ISITE  ASSWS  AS";  the  computer  program  is  listed  in 
Appendix  P. 

b.  Purpose 

This  function  is  called  by  the  SITE  program  to 
perform  the  target  damage  assessment  event. 

c.  Input 

The  left  argument  (ISITE)  is  a  vector  whose 
first  element  is  the  radar  site  number  and  whose  second  is 
the  time  step. 

2.  The  right  argument  (AS)  is  the  number  of  the 

target. 
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d .  Output 


The  output  is  the  time  required  to  assess  target 
damage.  This  is  used  to-Nschedule  the  assessment  future 
event. 

8.  PKILL 

a.  Header 

The  header  of  the  function  "PKILL"  is 
"PK  <-  PKM  PKILL  R";  the  computer  program  is  listed  in 
Appendix  F. 

b.  Purpose 

This  function  is  called  by  the  SITE  program  to 
calculate  the  probability  of  kill  for  each  radar  site  and  the 
effectiveness  of  ECM. 

c.  Input 

1.  The  first  argument  (PKM)  is  the  PK  of  the  ECM 
matrix  given  by  the  Probability  of  Killing  model  and  ECM 
model  respectively. 

2.  The  second  argument  (R)  is  a  two  element 
vector  whose  first  element  is  the  number  of  the  target  and 
whose  second  is  the  distance  between  the  radar  site  and  the 
target. 

d.  Output 

1.  The  probability  that  a  missile,  launched  by  a 
radar  site,  will  kill  an  aircraft. 

2*  Tha  probability  that  the  acquisition  achieved 
by  a  radar  site,  will  cease  because  of  electronic 
countermeasures . 
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9.  RANDOM 


a.  Header 

The  header  of  the  function  "RANDOM"  is 
"RN  <-  RANDOM"?  the  computer  program  is  listed  in  Appendix  F. 

b.  Purpose 

This  function  is  called  by  the  SITE  program  to 
generate  a  uniform  random  variable  between  zero  and  one. 

c.  Input 
None. 

d.  Output 

The  output  is  a  random  number  used  to  determine 
probabilities  when  a  stochastic  output  is  required. 
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IX.  SIMULATION  OF, A  PARTICULAR  SCENARIO 


A.  INTRODUCTION 

The  technique  of  simulating  aircraft  attacks  against  a 
ground  target  and  their  interception  by  a  missile  site  was 
described  in  Chapter  II  through  VIII.  In  this  chapter,  a 
particular  scenario  is  described  and  simulated  with  various 
combinations  of  inputs  to  calculate  aircraft  attrition. 

It  also  summarizes  all  the  steps  which  must  be  followed  to 
simulate  the  scenario  and  to  calculate  the  aircraft  attrition 
rate.  Furthermore,  the  output  data  given  by  the  airmodel 
using  this  particular  scenario,  will  be  analyzed  and 
subjected  to  a  sensitivity  analysis  of  critical  variables. 

B.  DESCRIPTION  OF  THE  SCENARIO 

The  particular  scenario  used  for  simulation  is  shown  in 
Figure  9.1. 

1 .  Map 

The  map  used  for  the  scenario  design  is  an  Air  Force 
map.  It  is  used  by  pilots  to  plan  aircraft  routes  for 
specific  missions.  The  map  scale  is  1:500,000  and  has 
sufficient  contour  lines  to  read  the  terrain  altitudes  and  to 
plan  the  aircraft  routes  so  as  to  take  advantage  of  terrain 
masking.  A  162  NM  by  120  NM  area  is  chosen  to  represent  the 
enemy  country  which  the  attacking  formations  will  fly  over. 


The  area  is  divided  into  grid  squares.  The  time 


steps  chosen  for  this  simulation  are  10  seconds. 

2.  Aircraft  Formations 

There  are  three  aircraft  formations  simulated  in  this 
particular  scenario. 

a.  Aircraft  Formation  1  (Target  1) 

(1)  Type.  The  aircraft  formation  1  is  an  enemy 
formation  and  it  consists  of  9  bomber  aircraft. 

(2)  Purpose.  The  mission  of  the  formation  is  to 
attack  an  airport  whose  coordinates  are  X  =  90  MM  and  Y  -  9 
NM  inland.  This  target  will  fly  over  the  NIKE  radar  site. 

(3)  Initial  point  and  Start  Time.  Formation  1 
enters  the  enemy  terrain  at  X  *  0  NM,  Y  a  111  NM  point  and  at 
time  zero.  (This  is  the  first  target  entering  the 
simulation) . 

(4)  Flight  Data.  The  flight  routes  of  formation 
1  consist  of  eight  straight  lines  and  the  flight  data  of  each 


flight  route 

are  listed 

below: 

Route 

Course 

(degrees 

Length 
)  (NM) 

Speed 
(knots ) 

Altitude 
( feet ) 

1 

17 

35 

2120 

1000 

2 

120 

44 

420 

1000 

3 

157 

27 

420 

1000 

4 

165 

46 

480 

1000 

5 

334 

38 

480 

2000 

6 

004 

30 

420 

1000 

7 

342 

28 

420 

1000 

8 

270 

69 

420 

10,000 
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(5)  Comments.  The  speed  in  routes  4  and  5  is 
increased  from  420  to  480  knots  because  the  aircraft 
formation  is  in  its  attack  approach  at  this  time.  The 
altitude  in  route  5  is  increased  from  1000  ft  to  2000  ft 
because  the  aircraft  formation  needs  to  be  formed  again  after 
the  airport  attack.  In  route  8  alititude  is  increased 
because  the  aircraft  formation  need  to  save  fuel  in  order  to 
return  to  its  base. 

b.  Aircraft  Formation  2  (Target  2) 

(1)  Type.  The  aircraft  formation  2  is  a  friendly 
formation  and  it  consists  of  four  interceptor  airplanes. 

(2)  Purpose.  The  mission  of  this  formation  is  to 
intercept  enemy  aircraft  and  return  to  its  base  (X  =  163  NM, 

Y  =  2  NM)  for  refueling.  It  penetrates  both  "HAWK"  and 
"NIKE"  radar  sites  and  is  included  here  to  demostrate 

the  accidental  attrition  of  friendly  aircraft  due  to  human 
errors  ( fractricide ) . 

(3)  Initial  Point  and  Start  Time.  Target  2 
enters  the  area  at  X  =  0  NM,  Y  *  81  NM  point  6  minutes  after 
time  zero. 


(4)  Flight  Data.  The  flight  routes  of  Formation 
2  consist  of  three  straight  lines  and  the  flight  data  of  each 
flight  route  are  listed  below: 
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Route 

Course 

(degrees) 

Length 

(NM) 

Speed 
(knots ) 

Altitude 
( feet) 

1 

109 

37 

420 

20.000 

2 

118 

44 

420 

20.000 

3 

131 

64 

420 

10.000 

(5)  Comments .  The  altitude  in  route  3  changes 
from  20,000  ft  to  10,000  ft  because  the  aircraft  formation 
descends  for  landing. 

c.  Aircraft  Formation  3  (target  3) 

(1)  Type .  The  Aircraft  Formation  3  is  an  enemy 
formation  and  it  consists  of  11  bomber  airplanes. 

(2)  Purpose.  The  mission  of  the  formation  is  to 
to  attack  the  airport  whose  coordinates  are  X  »  130  NM,  Y  - 
32  NM  inside  enemy  territory.  This  target's  flight  plan 
takes  it  through  two  overlapping  HAWK  radar  sites.  This  will 
allow  an  investigation  of  the  command  and  control  model. 

(3)  Initial  Point  and  Start  Time.  The  aircraft 
formation  enters  the  enemy  area  at  the  point  whose 
coordinates  are  X  *  0  NM,  Y  *  31  NM  and  at  5  minutes  after 
time  0. 

(4)  Flight  Data.  The  flight  routes  of  aircraft 
Formation  3  consist  of  seven  straight  lines  and  the  flight 
data  of  each  route  are  listed  below: 


Route 

Course 
(degrees ) 

Length 

(NM) 

Speed 
(knots ) 

Altitude 
( feet ) 

1 

088 

40 

420 

1000 

2 

050 

35 

420 

1000 

3 

111 

31 

420 

1000 

4 

082 

28 

480 

1000 

5 

249 

40 

480 

2000 

6 

299 

45 

420 

1000 

7 

280 

49 

420 

1000 

(5)  Comments .  The  speed  in  route  4  and  5  changes 
from  420  to  480  knots  because  the  formation  is  in  its  attack 
phase.  The  altitude  in  route  3  increases  from  1000  ft  to 
2000  ft  because  the  formation  needs  to  regroup  after 
attacking  the  target. 

3 .  Radar  Sites 

In  this  particular  scenario,  there  are  three  HAWK 
radar  sites  and  one  NIKE  site.  The  NIKE  site  covers  almost 
all  the  country;  the  first  two  HAWK  sites  are  overlaping  and 
cover  the  entrance  to  the  defended  area.  The  third  HAWK  site 
protects  a  very  important  industrial  town  inside  the  country. 

The  parameters  of  four  radar  sites  are  shown  on 

Table. 


T 


Radar 

Site 


Position 
(X,  Y) 


Alt. 

( feet) 


Maximum 

Detection 

Range 

(NM) 


1  1 

HAWK1 1 (25,26)  I 
|  | 

1250  1 

j 

25 

1 

I 

20 

1 

1 

5 

1 

1 

1 

HAWK2I (37,58)  1 
j  j 

725  1 

40 

1 

1 

| 

30 

1 

1 

1 

8 

1 

1 

1 

2 

HAWK 3 | (132,62) | 
j 

1550  | 

40 

1 

1 

1 

30 

1 

1 

1 

3 

1 

1 

1 

2 

NIKE  1(85,87)  | 

1  1 

1  1 

2587  | 

1 

! 

25 

1 

1 

1 

1 

87 

1 

1 

1 

1 

10 

1 

1 

1 

1 

1 

Maximum 
Range  of 
Used 
Missile 
(NM) 


Radar 

Missile 

Capacity 


Number 

of 

Engaged 

Targets 


C.  SUMMARY  STEPS  OF  CALCULATING  THE  ATTRITION  RATE  OF 

ATTACKING  AIRCRAFT 

In  this  particular  part  of  chapter  IX,  a  summary  of  the 
steps  which  are  required  in  order  to  calculate  the  attrition 
rate  of  the  attacking  aircraft  formations  are  described.  The 
entire  process  is  divided  into  two  parts:  preprocessing  and 
dynamic  simulation. 

1 .  Preprocessing  Phase 

Step  1:  Initialize  the  scenario.  The  user  must 
specify  a  scenario  similar  to  the  scenarifb  described  above. 

Step  2:  Calculate  the  average  terrain  altitudes  by 
using  the  instructions  of  Chapter  II  and  the  computer  program 
"TERRAIN"  given  in  Appendix  A.  An  output  of  the  Average 
Terrain  Model,  using  the  data  of  the  scenario  described 
above,  is  shown  in  Appendix  G. 

Step  3:  Calculate  the  minimum  detection  altitudes 


for  each  radar  site  by  using  the  instructions  of  Chapter  III 


and  the  computer  program  "MINALT"  given  in  Appendix  B.  An 
output  of  the  Minimum  Altitudes  Model,  using  the  data  of  the 
scenario  described  above,  is  shown  in  Appendix  G. 

Step  4:  Calculate  the  flight  data  of  each  aircraft 
formation  by  using  the  directions  of  Chapter  V  and  the 
computer  program  "CPM"  given  in  Appendix  D.  An  output  of  the 
Flight  Data  Model,  using  the  data  of  the  scenario  described 
above,  is  shown  in  Appendix  G. 

Step  5*  Calculate  the  interception  data  by  using  the 
instructions  of  Chapter  VI  and  the  computer  program  "INTER" 
given  in  Appendix  E.  An  output  of  the  Interception  Data 
Model,  using  the  data  of  the  scenario  described  above,  is 
shown  in  Appendix  G. 

Step  6:  Calculate  the  priority  matrix  by  using  the 
instructions  of  Chapter  VI  and  the  computer  program 
“PRIORITY"  given  in  Appendix  E.  An  output  of  the  Priority 
Matrix  Model,  using  the  data  of  the  scenario  described  above, 
is  shown  in  Appendix  G. 

2 .  Dynamic  Phas .» 

Step  7s  Run  the  "AIRMODEL"  program,  using  the 
instructions  of  chapter  VIII  and  the  computer  program 
“AIRMODEL"  given  in  Appendix  F.  An  output  of  the  "AIRMODEL" 
program,  using  the  data  of  the  scenario  described  above,  is 
shown  in  Appendix  G. 


D.  EXECUTING  THE  AIRMODEL 


1 .  Database 

This  portion  of  Chapter  IX  describes  and  analyzes  a 
complete  simulation  process  following  steps  one  through 
seven.  The  scenario  which  has  been  used  to  generate  this 
output  is  described  in  the  beginning  of  this  chapter.  The 
constant  inputs  of  the  computer  program  "MAIN",  in  the 
"AIRMODEL"  computer  workspace,  are  as  follows: 

a.  seed  *  various  seeds  have  been  used 

b.  times  =  50  times 

c.  NSITES=  4 


d. 

DS2 

a 

listed 

in 

Appendix  G 

e. 

IS2 

a 

listed 

in 

Appendix  G 

f. 

TGT 

= 

9 

4 

11 

g* 

MKAP2 

= 

1 

2 

2 

1 

h. 

MCAP2 

a 

5 

8 

8 

10 

1 

1 

5 

20 

25 

0.5 

1 

1 

1 

5 

25 

30 

0.55 

1 

1 

i. 

PKM 

'= 

1 

1 

1 

l 

5 

25 

30 

0.55 

1 

1 

1 

1 

10 

70 

85 

0.4 

1 

1 

j* 

IFFFOE 

a 

1 

1 

1 

5 

30 

80 

0.7 

0.85 

k. 

ECM 

= 

1 

1 

5 

40 

70 

0.5 

0.8 

l 

1 

10 

50 

80 

0.7 

0.9 

1. 

PR 

= 

listed 

in 

Appendix  G. 

m. 

AUTO 

32 

12 

12  3  2 


are: 


p.  TSTEP1® 


2  2  3  2 


I  1  2  3  2  1 

I  I 

|2  3  3  3  1 

o.  ETSF  =*  50% 

The  inputs  which  will  be  varied 


in  the  airmodel 


IFF 


Probability  that  IFF  works 

Probability  that  IFF  is  correct  given  that  it  works 
Command  and  Control  Model 
Radar  Site  Performance  Model 
Electronic  countermeasure  means  model 
2 .  Result 

a.  Base  Line  Situation 

The  airmodel  has  been  executed  with  the  following 
various  parameters: 

IFF  =  Yes 

Probability  that  IFF  works  =  0.98 
Probability  that  IFF  is  correct,  given  that  it 
works  =  0.96 

Command  and  Control  Model  =*  Yes 
Radar  Operation  =  Normal 
ECM  =  Effective 

The  setting  of  these  parameters  represents  a 
situation  which  is  feasible  during  conflict  operations.  The 
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output  -of  the  airmodel  is  listed  in  Appendix  G.  The  data 
analysis  is  as  follows: 


(1)  Summary  Statistics.  Some  summary  AC 

Attrition  statistics  are  listed  below.  This  data  is  result 

replicting  the  simulation  50  times. 

1 

1  Statistics 

1 

r 

Target  1  I 

1 

Target  2 

Target  3 

I 

1 

1 

1 

!  Number  of 

1 

1 

1  Killed  a/c 

0.64  | 

0.46 

2.2 

1 

1  (average) 

| 

1 

[ 

1 

1 

1  Sample  Variance 

0.72  1 

| 

0  •  25 

1.106 

1 

1 

[ 

1  Minimum  Value 

1 

0  1 

1 

0 

0 

1 

1 

i 

1  Maximum  Value 

1 

3  1 

1 

i 

5 

1 

1 

1 

1  Range 

1 

3  1 

i 

5 

1 

1 

| 

1  Median 

1 

0  1 

1 

0 

2 

1 

i 

i 

1  Low  Quantile 

1 

0  1 

o 

1 

1 

1 

1 

1 

1  Upper  Quantile 

1 

1  1 

1 

i 

3 

1 

1 

.1 

(2)  Distribution  Fitting.  The  most  appropriate 
distribution  for  the  attrition  data  of  the  three  targets  is 
the  binomial  with  the  following  parameters: 


Parameters 


n 

P 

X  square 


Target  1 


9 

0.058 

1.28 


Target  2 


4 

0.115 

3.56 


Target  3 


11 

0.24 

0.87 


2 

Since  the  critical  value  of  X  (alpha  =  0.05) 
to  reject  the  null  hypothesis  (Hos  F(X)  *  Binomial 
Distribution)  is  equal  to  7.81  for  two  degrees  of  freedom,  we 
accept  that  the  distribution  of  the  attrition  rate  data  is 
binomial  with  the  parameters  shown  above. 

(3)  Confidence  Intervals  for  Probability  of  Kill. 
Confidence  intervals  for  the  probability  of  kill  are  given 
below  for  alpha  of  0.05. 


T 

1 

1 

r 

I 

i 

Confidence 

Interval 

Target 

1 

1 

1 

Probability  ( 
Function  1 

1 

Low 

1 

1  Upper 

1 

1 

1 

P(X=0}=0. 56  | 

0.42 

1  0.69 

1 

1 

P { X=1 } =0 . 28  | 

0.15 

I  0.40 

Target 

1 

1 

1 

1 

1 

P{X«2 }-0. 12  1 

1 

0.02 

1 

1  0.21 
i 

1 

1 

1 

1 

P{X«3}=0.04  1 

1 

0.0 

1 

I  0.09 

i 

1 

P(X~0 }=0. 54  | 

0.4 

1  0.67 

Target 

2 

1 

1 

1 

1 

P  [X=l }=0. 46  | 

1 

0.32 

1 

1  0.59 

1 

1 

1 

P{X»0}»0.04  | 

i 

0.0 

1  0.09 

| 

1 

1 

1 

P{X=1 }=0. 24  | 

1 

0.12 

I  0.35 

1 

1 

1 

1 

P{X=2}=0.34  1 

0.208 

1 

I  0.47 

Target 

3 

1 

1 

1 

1 

P{X»3}«0.26  | 

0.13 

i 

1  0.38 

1 

1 

1 

1 

P{X-4}«0.1  | 

i 

0.01 

l 

I  0.18 

1 

1 

1 

1 

1 

P {X=5 ] =0 . 02  | 

1 

0.0 

1 

I  0.05 

1 

b.  Other  Ouputs 

In  this  section,  some  other  outputs  of  "AIRMODEL 
are  given.  These  outputs  are  a  result  of  changing: 

1.  Probability  that  IFF  works  and  probability 
that  IFF  is  correct,  given  that  it  works. 

2*  Command  and  Control  model. 

3.  Radar  site  operation  model. 

4.  Electronic  countermeasure  means  model. 

(1)  Changing  the  IFF  probabilities.  In  this 
case,  the  "AIRMODEL”  computer  program  uses  the  following 
input  variables. 

IFF  *  Yes 

Probability  that  IFF  works  =  0.90 

Probability  that  IFF  is  correct  given  that  it 
works  3,1  0.90 

Command  and  Control. Model  =  Yes 

Radar  Operation  =  Good 

Electronic  countermeasure  means  model  =  Yes 

These  inputs  are  the  same  as  in  the  Base  Line 
Situation,  above,  except  the  probabilities  of  IFF. 

The  means  of  the  output  of  the  target  attrition 


rate  data  are: 


#  of  A/C 


Target  1 


Target  2 


Target  3 


Probabilities 


1  1 

1 

Mean  of  1 
Attrition! 
Rate  1 

0.44 

1 

1 

1 

1.42 

1 

1 

1 

1 

0.64 

1 

1 

1 

1 

0.46 

2.06 


IFFWORK  =0.90 
IFFWORK  =0.90 


2.2 


IFFWORK  =0.98 
IFFWORK  =0.96 


This  table  shows  that  the  number  of  attrited 
aircraft  is  less  in  both  targets  one  and  three  because  more 
enemy  aircrafts  pass  without  interception  (recall  that  the 
IFF  model  assumes  that  if  IFF  does  not  work,  the  detected 
target  is  considered  to  be  friend).  On  the  other  hand,  the 
the  number  of  destroyed  aircraft  of  target  two  is  greater 
decause  the  IFF  has  a  lower  degree  of  reliability. 
Consequently,  there  are.  more  fractricides .  (Recall  that 
target  two  is  a  friendly  aircraft  formation  which  is 
returning  for  a  landing). 

(2)  Changing  Command  and  Control  Input.  In  this 
case,  there  is  no  command  and  control  model  applied.  That 
means  that  the  radar  sites  operate  autonomously.  The  input 
variables  of  "AIRMODEL"  ares 
IFF  =  Yes 
IFFWORK  =  0.98 
IFFCOR  -0.96 

Command  and  Control  Model  =  No 
Radar  Operation  =  Good 
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Electronic  countermeasure  means  *  Yes 


These  inputs  are  the  same  as  in  the  Base  Line 
Situation,  except  that  the  command  and  control  model  is  not 
applied. 

The  table,  below,  shows  that  the  number  of 
destroyed  aircraft  is  greater  when  command  and  control  is  not 
applied.  This  result  is  a  suprise  because  the  opposite  was 
expected. 

The  means  of  the  output  of  target  attrition  are 

given  below: 


I - 

1 

1  #  Of  A/F 

T 

1 

1 

1 

- r 

i 

Target  1  1 

1 

- r 

1 

Target  2  1 

1 

- r 

i 

Target  3  1 

1 

Command  and 
Control  model 

1 

1 

1  Attrition 

1 

1 

1.08  | 

1 

H 

• 

H 

2.36  | 

1 

not  applied 

1  Rate 

1 

1 

1 

J_ 

0.64  | 

1 

1 

0.46  1 

1 

2.2  | 

1 

applied 

However,  the  result  came  about  because  the 
scenario  has  only  three  targets  and  four  radar  sites.  When 
one  target  formation  was  engaged  by  radar  site  one,  radar 
site  two  was  not  allowed  to  engage  this  target  time  (IDLE). 
Therefore,  the  formation  was  intercepted  only  by  radar  site 
one  and  it  suffered  less  attrition.  On  the  other  hand,  when 
command  and  control  is  not  applied,  the  same  target  formation 
is  intercepted  by  more  than  one  radar  site  and  it  suffers 
more  attrition  (both  radar  sites  are  active).  This 
phenomenon  would  not  happen  if  the  scenario  had  more  targets. 


(3)  Changing  the  radar  site  operation  model.  In 
this  case,  the  operation  of  all  the  radar  sites  is  degraded. 
The  variable  inputs  of  "AIRMODEL"  are: 

IFF  *  Yes 

IFFWORK  »  0.98  * 

IFFCOR  -0.96 

Command  and  Control  model  »  Yes 
Radar  Operation  ■  degraded 
Electronic  countermeasure  means  *  Yes 
These  inputs  are  the  same  as  in  the  Base  Line 
Situation,  above,  except  that  all  radar  operations  are 
degraded. 

The  means  of  the  output  of  target  attrition  rate 
data  are  shown  below: 


*  of  A/C 

1 

1 

.1. 

1 

Target  1  1 

1 

1 

Target  2  1 

Target  3  1  Radar  Operation 

1 

Attrition 

1 

1. 

0.60  1 

1 

0.41  | 

1 

1.42  1 

1 

Degraded 

Rate 

1 

1 

0.64  I 

1 

0.46  | 

1 

2.2  1 

1 

Good 

The  table  shows  that  the  number  of  destroyed 
aircraft  for  all  three  targets  is  less,  as  it  is  expected, 
when  the  radar  sites  work  in  a  degraded  mode.  This  attrition 
rate  has  been  reduced  almost  40%.  in  the  case  of  Table  3. 


126 


(4)  Changing  the  electronic  countermeasure  means 


model*  In  this  case,  the  electronic 
countermeasure  means  model  (ECM)  is  not  applied.  That  means 
that  no  target  formation  uses  ECM. 

The  variable  inputs  of  "AIRMODEL"  are: 

IFF  -  Yes 
IFFWORK  *  0.98 
IFFCOR  *  0.96 

Command  and  Control  Model  *  Yes 
Radar  Operation  ■  Good 

Electronic  countermeasure  means  model  =  No 
These  inputs  are  the  same  as  the  Base  Line 
Situation,  above,  except  that  the  electronic  countermeasure 
means  model  is  not  applied. 

The  means  of  the  output  of  target  attrition  rate 
data  is  shown  below: 


1 

1  #  Of  A/C 

1 

i - r 

1  Target  1  1 

1  1 

- p 

Target  2  1 

1 

- r 

Target  3  1 

1 

ECM  Model 

1 

1 

1  Attrition 

1  0.76  1 

1  1 

0.93  1 

1 

4.06  1 

1 

not  applied 

1  Rate 

1 

1 

i  0.64  j 

1  1 

0.46  ! 

1 

2.2  1 

1 

applied 

The  table  shows  that  the  number  of  destroyed 
aircraft  is  greater  when  the  targets  do  not  use  electronic 
countermeasure  means  to  protect  themselves.  The  increase  of 
number  of  the  destroyed  aircraft  when  ECM  is  not  applied,  is 
approximately  100%. 


3.  Sensitivity  Analysis 

By  using  the  scenario  described  above,  sensitivity 
analysis  can  be  performed  by  changing  some  of  the  aircraft 
movement  data.  For  example,  by  changing  the  routes  that  the 
aircraft  formation  use,  one  can  see  how  the  attrition  rate 
will  be  changed.  Other  parameters  that  can  be  changed  for 
sensitivity  analysis  are  the  flight  altitude  and  the  speed  of 
the  aircraft  formation. 

In  this  section,  sensititvity  analysis  is  performed  by 
changing  the  flight  altitude  of  the  aircraft  formations.  A 
change  in  the  flight  altitude  was  performed  to  study  the 
effects  of  terrain  masking  on  aircraft  attrition  rate. 

To  perform  this  sensitivity  analysis,  five  different 
aircraft  flight  levels  were  used  (500,  1000,  2000,  5000  and 
10000  feet).  The  outputs  of  the  "AIRMODEL",  using  these  five 
different  flight  level  data,  are  listed  in  Appendix  G. 

The  number  of  destroyed  aircraft  for  each  different 
flight  level  are  given  below: 


Flight 

Level 

Target  1 

1 - 

1 

1  Target  2 

I 

1 - 

1 

1  Target  3 

1 

1 

1 

1  Total 

1 

500 

0.83 

1  0. 13 

i 

1  0.77 

1 

1  1.73 

1000 

0.64 

1  0.46 

1 

1  2.20 

1  3.30 

1 

2000 

0.82 

1 

1  0.62 

1 

1  2.25 

1 

1 

1  3.69 

5000 

0.94 

1 

1  0.65 

i 

1  2.16 

I  3.75 

1 

10000 

0.96 

1  0.66 

I 

1  2.20 

1 

1 

1  3.82 

i 

I 

1 


The  table  shows  the  increase  in  the  number  of  destroyed 
aircraft  as  the  flight  level  increases. 
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Figure  9.2  above#  shows  the  change  in  attrition  rate  of 
aircraft  as  a  function  of  the  flight  level  for  each  target. 


Figure  9.3  below,  shows  the  change  in  the  total  attrition 
rate  of  the  three  formations  as  a  function  of  the  flight 
level. 

T 


Figure  9.3  Total  Attrition  Rate 

It  is  worth  while  considering  that  the  increase  in 
attrition  rate,  especially  for  target  three,  is  very  large 
when  the  formation  changes  altitude  from  500  feet  to  1000 
feet.  However,  the  attrition  rate  is  essentially  constant 
for  altitudes  above  2000  feet. 


X.  SUMMARY  AND  CONCLUSION 


A.  SUMMARY 

This  thesis  has  presented  a  high  resolution  model  to 
analyze  the  air  attack/defense  combat  situation. 

The  simulation  has  made  efficient  use  of  off-line, 
pre-processing  submodels  and  programs  in  order  to  compute 
the  deterministic  phase  of  the  engagement.  By  pre-flying 
the  aircraft  through  the  entire  scenario,  most  of  the 
detailed  calculations  involving  the  geometry  of  the  air 
engagement  were  performed.  These  calculations  were  used  to 
determine  if  detection,  tracking,  firing,  and  interception 
are  feasible  during  all  segments  of  the  flight.  The 
information  was  summarized  and  stored  in  the  form  of 
matrices  whose  only  elements  are  zero  and  one.  As  a 
result,  the  deterministic  elements  of  the  air  battle  were 
transformed  into  look-up  table  for  each  time  step  of  the 
simulation. 

The  pre-processing  phase  included: 

1.  Generating  the  terrain  models  (Chapters  II  and  III) 

2.  Modeling  the  flight  paths  of  all  aircraft  contained 
in  the  simulation (Chapter  V) 

3.  Completing  detection,  tracking,  firing  and 
interception  calculations  (Chapters  IV  through  VI). 


4.  Modeling  for  the  effects  of  IFF,  degraded  equipment, 
ECM,  target  selection  priorities,  and  command  and  control 
(Chapter  VII). 

The  results  of  the  pre-processing  calculations  were 
input  into  the  dynamic  simulation  phase  of  the  model.  The 
stochastic  portions  of  the  air  attack/defense  encounter 
were  replicated  and  summary  statistics  were  gathered  and 
analyzed.  A  sensitivity  analysis  was  performed  to 
demonstrate  the  relationship  between  aircraft  altitude  and 
attrition  rate.  (Chapter  VIII  and  IX) 

B.  CONCLUSION 

The  technique  of  separating  the  deterministic  and 
stochastic  elements  of  the  air  battle  proved  to  be  very 
efficient.  Using  an  IBM  PC,  it  required  approximately  25 
minutes  to  pre-fly  the  aircraft  and  summarize  the  geometry 
of  the  air  battle  in  matrix  format.  The  stochastic  portion 
of  the  air  battle  required  approximately  8  minutes  of 
executing  time.  Although  these  run  times  are  scenario 
dependent,  the  three-to-one  ratio  is  a  reasonable 
approximation  to  the  apportionment  of  time  required  for 
both  phases  of  the  simulation.  Replicating  the  entire 
simulation  50  times,  without  a  pre-processing  oriented 
model,  would  require  1650  minutes  of  IBM  PC  time.  Using 
the  pre-processing  approach,  the  time  required  for  the  same 
50  replications  is  425  minutes.  This  represents  a  savings 
of  approximately  75%. 


The  submodels  developed  in  this  thesis  can  be  used  in 


several  ways. 

1.  The  attacker  can  use  the  pre-processing  models  to 
design  flight  paths  which  minimize  detection  and  take  more 
advantage  of  terrain  masking. 

2.  The  attacker  can  perform  trade-off  analyses  between 
different  approach  altitudes.  For  a  given  flight  path,  the 
attacker  can  determine  the  highest  altitude  which  will 
offer  deception  and  cover  and  allow  for  efficient  fuel 
consumption. 

3.  The  attacker  can  better  determine  the  missile  sites 
which  offer  the  greatest  threat  of  attrition.  Air  assets 
can  then  be  allocated  to  attack  these  sites  and  degrade 
their  capability. 

4.  The  defender  can  use  the  pre-processing  models  to 
improve  the  positioning  of  air  defense  sites,  given  the 
likely  axis  of  advance  by  attacking  aircraft. 

AIRMODEL  allows  the  analyst  to  investigate  a  wide  range 
of  aircraft  attack/defense  problems.  The  simulation  can  be 
efficiently  replicated  without  re-calculating  the 
deterministic  phases  of  the  engagement.  Consequently,  a 
more  realistic  and  accurate  perspective  of  the  encounter 
can  be  achieved. 

Further  research  and  analyses  can  be  directed  towards 
developing  the  stochasting  sub-models  necessary  to  simulate 
the  engagment  sequence.  This  author  has  developed  a 
general  approach  which  is  efficient  and  realistic. 
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Acquisition,  tracking,  and  firing  time  delays  require  more 
study  and  model  design.  The  user  is  invited  to  employ  this 
modeling  architecture  as  a  specific  framework  for  model 
building.  The  entire  simulation  model  is  modular. 
Therefore,  the  user  can  modify  the  sub  models  to  suit 
specific  design  requirement  and  level  of  fidelity. 


APPENDIX  A 


« 

* 

M 


A.  PURPOSE 

This  program  is  used  to  calculate: 

1.  The  average  altitudes  for  each  grid  square. 

2.  To  construct  a  matrix  named  ATERRAIN  whose  entries  are 
the  average  altitudes  for  grid  squares. 

B.  DESCRIPTION  OP  THE  TERRAIN  SIMULATION  PROGRAM 

The  header  of  the  terrain  program  is  "TERRAIN". 

1 .  Input 

The  "TERRAIN"  program  prompts  one  input: 

A  matrix  A  whose  entries  are  the  altitudes  of  each 
grid  square.  How  these  numbers  are  extracted  is  described  in 
section  one.  Matrix  A  does  not  need  to  be  square.  The 
number  of  rows  is  the  result  of  the  division  of  the  length  of 
X  range  by  the  length  of  the  side  of  the  square  grid.  For 
example  if  the  length  of  the  terrain  is  60  NM  and  the  length 
of  the  side  of  the  grid  square  is  6  NM,  the  number  of  rows  of 
the  matrix  will  be  equal  to  10.  The  number  of  columns  is 
calculated  in  the  same  way. 

2 .  Output 

The  "TERRAIN"  program  calculates  a  matrix  whose 
entries  are  the  average  heights  of  grid  squares  such  as  shown 
previously  in  Table  2. 
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'  •r.  -  . 


/..A  A  • 


avv-j*: 


7  TERRAIN; B; C; I ; A 

m  a 

C2]  A  THIS  FUNCTION  CALCULATE  THE  AVERAGE  OF  THE  HIGHTS  OF  THE  GIVEN  TERAIN 

C3]  A -  INPUT  DESCRIPTION - 

C4]  A  THE  REQUARED  INPUT  IS  A  MATRIX  WHOSE  ENTRIES  ARE  THE  ALTITUDE 
CS3  A  OF  THE  FOUR  CORNERS  OF  EACH  SQUARE  GRID 
ca:  a 

C7]  A - OUTPUT  DESCRIPTION - 

C3I  A  THE  OUTPUT  IS  A  MATRIX  WHOSE  ENTRIES  ARE  THE  AVERAGE  ALTITUDE 
C9]  A  OF  THE  SQUARE  GRID 

CIO!  A - 

Cll]  A - MAIN  PROGRAM - 

C 12]  I-s-0 

C 13]  PSs 'ENTER  THE  MATRIX  OF  THE  SQUARE  GRIND  ALTITUDE' 

C 14]  A*d  4  *< <+/ <pA) )<2> /P3  4  *P4 

CIS]  P3s  THE  INPUT  IS  NOT  A  MATRIX) '  4  +P5 

C  16 ]  A - CREATE  THE  MATRICES - 

C  17]  P4:B*( < (pAC; 1]>-1> , (pAtl; ]> )pO  4  C+<  (  (pAC;  1  ] » —  1 )  ,  (  (pAClj  ]>-l)  >pO 
CIS]  Lis  1*1  +  1  4  *<I> <p <ACs 1])) > /L2 
C  19]  BCI; ]*ACI; ]+AC (1*1)  ;  ] 

COO]  *L1 
C21]  L2s 1*0 

C22I  L3: 1*1+1  4  *<I> (p (BCl! ]) > ) /L3 
C23]  CC; I]+ (BC ; I]*BC ; I+t ]) *4 
C24]  +L3 

C2S1  LS : AHE I GHTS+C 

C26]  THE  MATRIX  ' ' AHEIGHTS  ''' 

C27]  - ■ 

C2S]  AHEIGHTS 
C29] 

7 
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APPENDIX  B 


This  appendix  describes  the  workspace  called  "MINALT" 
(MINimum  ALTitudes).  This  workspace  contains  two  functions 
called  "MALTI"  and  "GRID". 

A.  FUNCTION  "MALTI" 

HEADER  :  The  header  of  the  function  "MALTI"  is  MALTI. 
PURPOSE:  This  function  calculates  the  minimum  altitude 

where  a  given  radar  site  can  "see"  a  target  due 
to  terrain  masking  and  to  limitation  of  the 
maximum  range. 

INPUT:  The  function  "MALTI"  prompts  four  inputs: 

1.  A  three  element  vector  which  contains  the 
position  and  the  altitude  of  the  radar  site. 

(a)  X, Y  cordinates  in  square  grid  position. 
For  example,  if  the  radar  site  is  located 

in  square  grid  (7,3)  then  first  two  elements 
are  7  3. 

(b)  Radar  site  altitude  (in  feet).  Thus  the 
completed  input  is  7  3  1000 

2.  The  maximum  range  of  the  radar  site  in  NM 

3.  A  (N,M)  matrix  containings  the  average 
Altitudes  of  the  grid  squares  calculated 
by  the  program  "TERRAIN"  described  in 
Appendix  A  of  Chapter  2. 


4.  A  scale  factor,  describing  the  size  of  the 


OUTPUT: 


FUNCTION 
HEADER  : 

PURPOSE: 


INPUT: 


OUTPUT: 


grid  squares  desired.  Its  measurement  is  in 
Nautical  Miles  (N.M). 

The  output  is  a  (N  M)  matrix  whose  entries  are 
the  minimum  altitudes  above  which  a  given 
radar  site  can  detect  a  target. 

"GRID* 

The  header  of  the  function  "GRID"  is  SITE  GRID 
CORD. 

This  function  finds  the  grid  squares  which  exist 
between  the  radar  site  and  the  target  grid 
square.  For  example ,  if  the  radar  site  is 
located  in  grid  square  (2,1)  -  see  Figure  3.1  in 
Chapter  3  -  and  the  target  square  is  (4,4),  then 
the  intermediate  grid  squares  given  by  the 
function  "GRID"  are  (2,2),  (3,2)  and  (3,3). 

This  function  supports  the  previous  function 
"MALTI"  but  it  can  be  run  separately,  as  well. 
The  function  "GRID"  prompts  two  inputs: 

1.  A  three  element  vector  which  contains  the 
coordinations  of  the  radar  site  X,  Y  and  its 
altitude.  This  vector  is  the  left  argument 
of  the  function's  header  (SITE). 

2.  The  coordinations  of  the  target  square  grid 
X,Y.  this  input  is  the  right  argument  of  the 
function's  header  (CORD). 

The  output  of  this  program  is  a  (2  N)  matrix 


whose  columns  are  the  coordinates  of  the  square 


grids  which  lay  between  the  radar  site  and  the 
target  square  grid. 


7  SR<-SITE  GRID  CORD;  I  s  A;  XX 1  s  XX2:  Al  5  X 1 5  X2:  XX;  YY;  Y1 ;  Y2;  F;  VEC;  X ;  Y;  R 
m  a 

C2]  *»«—«»»—— —»—«»—  PURPOSE 

C3]  A  THIS  FUNCTION  IS  USED  TO  FIND  THE  SQUARE  GRIDS  WHICH  ARE 
C4]  A  BETWEEN  THE  RADAR  SITE  AND  THE  SQUARE  TARGET  GRID. 

C31  A 

C6]  A  mmmmmmmmmmmmmmmmmmmmmm  INPUT 

C7]  A  t.  THE  CORD I NAT IONS  OF  THE  RADAR  SITE:  X,Y  AND  ALTITUDE  (IN  FEET) 
C3]  A  2.  THE  CORD INAT IONS  OF  THE  SQUARE  TARGET  GRID:  X,Y  . 

C9I  A 

C 10 I  fimmmmmwmmmmmmmmmmmmmmm  OUTPUT  mmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmm* 

Clll  A  A  (2  N)  MATRIX  WHOSE  COLUMNS  ARE  THE  CORD I NATES  OF  THE  SQUARE 
C 101  A  GRIDS  WHICH  LAY  BETWEEN  THE  RADAR  SITE  AND  THE  SQUARE 
C101  A  TARGET  GRID. 

C  14]  A -*=»«-»»—«»*— mmmmm  MAIN  PROGRAM  »*»»■»=»=»«==«=====»»=*»»«==«*=«  + 
CIS]  1-1 

C  16 ]  X 1 -SI TEC  1 1-0. 3  »  XO-CORDC 11-0.3 
C 17]  Y1-SITEC01-0.3  4  Y2-C0RDC21-O.5 

CIS!  A - THE  DISTANCE - - 

C  191  R«-<  (  (X2-Xl)*2)f<  (Y2-Y1 )  *0)  )*0.3 

C201  A - THE  VECTOR - - 

C211  VEC*<  U  (IR)xS)  >*3 
C221  ■*<  CX2-Xl)*0)/L13 

C231  X2-X2+0. 1 

C 24 1  L13:F*C3)o<  (Y2-Y1 ) -XX2-X1  >  ) 

C231  A - X  Y  CALCULATION - * 

C261  X“  (2oF)  mvEC  0  Y*(10F);<VEC 

C271  *(  (  (X2-X1 )  JO)  ,  (  (Y2-YIX0)  )  /L3,L3 

C2S1  XX<-Xi-X  0  YY«rYl-Y  »  -*L6 

C291  L3:  -» <  (Y2-Y1X0)  /L3 

C301  YY*Yl+Y  ♦  XXt-Xi  +  X  0  “L6 

C311  L3 :  YY“Y  1  -Y  0  XX*X1-X 

C321  ■»< ( (X2-Xi)<0> , ( (Y2-YI) >0) ) /L6.L6 

C331  YY+Yi+Y  0  XX«-X1+X 

C34]  L6:A«i-r  (2,  (pVEC)  )pXX,YY 

C3S1  *< (pAClj ]>*0> /RR10 

C361  GR«-  2  1  AC  (SITEC  1 1  ,SITEC2] )  »  -LLL2 
C371  RR10: A1“AC: 11 

C081  A —  REJECT  THE  CASE  WHEN  THE  GRIDCI1  IS  =  GRIDCI  +  1] - ♦ 

C39]  L2:-< (AClj I ] ) - < AC l : 1*11) ) /LI 
C  40 1  -LI  2 

C 4 1 ]  Ll:-< (AC2: 1 1 )  =  ( AC2; I  +  ll) )/Lll 
C421  L12:Al-AltAC:I“l3 
C4.01  LI  1 :  I«-I  +  l 

C 44 ]  -<i  ;<AVEC))/L2  ■>  GR“Q  (  <  ( p A 1 )  -j-2)  ,2)  pAl 

CA51  XXi-GRCl; ] 

C461  XX2H3RC2:  1 

C47]  A -  REJECT  THE  CASE  WHEN  THE  GRID  IS  THE  SAME  WITH  THE  SITE  - * 

C481  -<  (ltXXl)  *  (LSITECU)  >  /LL5 

C49]  -( (ltXX2)A (ISITEC21) > /LL3 

C501  XX1*UXX1  »  XX2-UXX2  »  GR*  (2,  {ft  XX  1 )  )  f>XX  1 ,  XX2 

CS1 1  A - REJECT  IF  THE  CORDINATIONS  ARE  ZERO - * 

C32]  LL3:~< <'itXXl>-0> /LLL3 
C331  +1 ("ltXX2)-0)/LLL3 

C34]  -LLL2 

CS31  LLL3:  XXl-“l  +  XXi  »  XX2-”UXX2 
C361  GR“ <2, ( A  XX  1 ) ) A  <  XX 1 , XX2) 

C37]  LLL2: GR-GR.CORD 


m 


C693  -  9  MALT  I j C0S1 : C0S2; SIN1 :5IN2:L:P:Al;G;GR;DliD2i GGR; GO; GRM; H: I ; J  j  K ; R; SC 

ALEi ALT1 s  ALT2; CORD: D3s  NM; THs TH1 ; TH2; TH3; TH4;  TH3s  SSITE 
Cl]  A 

C23  A  PURPOSE  i««i***ll»*3=ma=czi3BizBZ=z:z3«iat 

C33  A  THIS  FUNCTION  CALCULATE  THE  MINIMUM  ALTITUDE  THAT  A  THE  GIVEN 
C43  A  RADAR  SITE  CAN  SEE  A  TARGET  DUE  TO  MAXIMUM  RANGE  AND  TO  TERRAIN 
C33  A  MASKING 
Col  A 

C73  A  ■=**»»*“»"*«»»*«*  INPUT 

C83  A  I.  A  THREE  ELEMENT  VECTOR  FOR  THE  RADAR  SITE  POSITION 
C93  A  (A) .  X  Y  CORDINATES  IN  SQUARE  GRID  POSITION 

C103  A  (B>.  RADAR  SITE  ALTITUDE  (IN  FEET) 

CUI  A 

C 123  A  2.  THE  MAXIMUM  RANGE  OF  THE  RADAR  SITE  IN  N.M 
C  131  A 

C143  A  3.  A  (N  M)  MATRIX  WITH  THE  AVERAGE  ALTITUDES  (GIVEN  BY  THE  TERRAIN 
C133  A  SIMULATION) 

C  163  A 

C 173  A  3.  THE  SCALE  FACTOR  (THE  LENGTH  OF  THE  SQUARE  GRIDS)  IN  N.M 
C  183  A 

C  193  A mmmmmmmmmmm  OUTPUT  mmmm.mmmmmmmmmmmmmmmmmmmmmmmmmmmmmmm* 

C203  A  THE  OUTPUT  IS  A  (N  M). MATRIX  WITH  THE  MINIMUM  ALTITUDES 
C213  A  ABOVE  OF  WHICH  A  GIVEN  RADAR  SITE  CAN  SEE  A  TARGET 
C223  A 

C233  A mmmmmmmmmmmmmmmmmmm  MAIN  PROGRAM 
C24  3  A 

C233  T2:  ENTER  THE  SITE  DATA  (X  ,Y, ALTITUDE)  ' 

C263  SITE-0  4  SSITE-(r < (2+SITE) -6) > ,*ltSITE  4  SITE  4  SSITE 
C273  •»  (  (pSSITE)  *3)  /T1 

C293  'THE  VECTOR  IS  NOT  CORRECT  !!!  TRY  AGAIN'  4  -T2  4  '  ' 

C293  '  ' 

C303  Til 'ENTER  THE  MAXIMUM  RANGE  OF  THE  RADAR  SITE  (IN  N.M) ' 

C313  RANGE-0 
C323  '  ' 

C333  T4:  ENTER  THE  AHEIGHTS  MATRIX' 

C343  H-0 

C3S3  -(  (pH)>l>  /T3 

C363  THE  MATRIX  DOES  NOT  HAVE  CORRECT  D I MENTIONS  TRY  AGAIN'  4  -T4 
C373  '  ' 

C393  T3i 'ENTER  THE  SCALE  FACTOR' 

C393  SCALE -0 

C403  J-t  4  HI- (pH) pO 

C41 3  PP4: 1-1 

C423  PP3 : CORD- ( I ) , ( J )  4  K-l  a4  'I,J  -  ',*I,J 

C433  GO- ( ( ( I -SCALE) - (SCALE-2) ) , ( ( J:<SCALE) - (SCALE-2) > )  A4  GO-'.JGO 
C443  R-3433  4  NM-6067  4  aSITE- ((< L SSITEC 1 , 23*SCALE> >-( SCALE-2) ), SS ITEC3 3 

C453  L- ( ( (GQC 1 3-5ITEC l 3  >  *2) - < (GOC23-SITEC23 ) *2) >*0.5 
C  463  - (Li RANGE ) /T4 

C47  3  HI C 1 5  J  3-” 1  4  -END 


C453  L- ( ( (GQC 1 3-5ITEC l 3  >  *2)  +  < (GOC23-SITEC23 ) *2) >*0.3 
C  463  -(LiRANGE) /T4 
C47  3  HI C I } J 3-” 1  4  -END 

C  49  3  A - 

C 49 3  A  FIND  THE  GRID  BETWEEN  THE  RADAR  SITE  AND  THE  TARGET  GRID 
C50  3  T4s  GR-SSITE  GRID  CORD 

C313  A - 

CS23  GR-«GR  4  GRM- (p (GRC 1 ; 3 ) ) pO  A4  GR 
C333  PP1 s  GGR-GRC ; K3 

C54  3  G-< ( <  OGR  C : K  3 ) .  -  ( SCALE ) )- (SCALE-2) ) ,HC (GGRC 13) :  (GGRC23) 3  A 4 


*D3 


C 68 3  ■»  (03*0)  /TT2 

C693  Al-'10  4  -»TT  1 

C703  TT2:D3><-<  <  <ALTt*2>*<ALT2*2n-<2*ALTlxALT2xC0Sl)>*0.3  44  '  D3- 

C713  TH3«-ll  (D3*2)*(ALTt*2>  >-<ALT2*2)  >  *  <2*D3>;ALTl  >  44  'TH3»',lTH3 

C723  TH4«-*2oTH3  44  TH4-',*TH4 

C73  3  TH3«-  <oi  >  -  (TH4*TH2)  40  'THS-',*TH5 

C74  3  SIN1-10TH5  4  SIN2*loTH4  44  SIN1.SIN2 

C733  Al«-<  (  <SIN2:<ALTt)+SINn-R)x<NM)  44  AAI-,*A1 

C763  4 - 

C773  TT 1 : GRM C K ] A l  44  At 

C783  K«-K+l  4  ■»  (K>  (pGRC  l ;  3  >  )  /PP2 

C793  ->>PP1 

cao:  - - 

cat:  4  TAKE  THE  LARGEST  ALTITUDE 

C323  PP2:H1CI;  J3--lt  (GRMC4GRM3)  44  ‘FIN  ALT-  ',»H1CI;J3 

C33  3  4  ' - ■ 

C843  END: !♦!♦! 


C3S3 

-*(IS  (pHC ;  13)  )  /PP3 

C  36  3 

1 

C873 

-*<JS  (pHC  1 ;  3)  )  /PP4  4 

HI  C  (SSI TEC  13);  ( SSITEC23  )  3**10 

CSS  3 
C89  3 

THE 

SMALLEST 

ALTITUDE  ARE' 

Hl-fHl 

C90  3 
C91  3 
C«23 

C  93  3 

' NOTE ’ 

C943 
C933 
C963 
C973 
C98  3 

THE  OUTPUT 

IS  STORED 

UNTER  THE  VARIABLE  HI ' 

9 
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APPENDIX  C 


Appendix  C  describes  the  computer  function  named 
“CONTOUR" . 

HEADER  :  The  header  of  the  function  CONTOUR  is  "CONTOUR " 
PURPOSE:  This  function  calculates  if  a  given  radar  site 
can  detect  a  target  flying  at  a  given  altitude 
level  for  all  the  square  grids  which  are 
located  inside  the  radar  site  maximum  distance. 
INPUT  :  The  computer  function  "CONTOUR"  prompts  three 
inputs : 

1.  A  (N  M)  matrix  given  by  the  terrain 
simulation  functions  "TERRAIN"  described  in 
Appendix  A. 

2.  The  altitude  of  the  target  (in  feet). 

3.  A  (N  M)  matrix  given  by  the  minimum  altitude 
simulation  function  "MALTI",  described  in 
Appendix  B. 

OUTPUT  :  The  output  is  a  (N  M)  matrix  the  entries  of 
which  are  -10,  -1,  0  and  1.  The  meaning  of 
these  numbers  are  as  follows: 

-10  »  Position  of  radar  site. 

-1  *  Radar  site  can  not  detect  the  target 

because  it  flies  out  of  the  radar  maximum 
range . 

0  ■  Radar  site  can  not  detect  the  target 


because  there  is  no  line— of-sight# 

1  »  Radar  site  can  detect  the  target. 
The  output  matrix  is  stored  under  the  variable  CON 


9  CONTOUR! ALT 1 s  ALT-! ALT3; ALT4; I ;  J :  H 


CU  * 

C23  H— ..«•*«— »——«»  PERPOSE  . . . . 

cci  *  This  function  calculate  a  table  which  illustrate  if  a 

C43  A  A  PARTICULER  RADAR  SITE  CAN  SEE  THE  TARGET  OR  NOT. 

C33  A 

C63  A  mmmmmmmmmmmmmmmmmmmm  INPUT 

C7 3  A  1.  A  (N  M>  MATRIX  OF  TERRAIN  ALTITUDE (FROM  TERRAIN  SIMULATION) 

C83  A  C.  THE  ALTITUDE  OF  THE  TARGET 

C93  A  3.  A  (N  M)  MATRIX  OF  RADAR  SITE  MINIMUM  ALTITUDE 
C103  A 

CU]  mmmmmmmmmmmmmmmmmmmm  OUTPUT  mmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmm ♦ 

C 121  A  A  (N  M)  MATRIX  WITH  ELEMENT  OF  *10  *1  0  AND  1. 

C  131  A 

C 1 4 3  A  *10  -  RADAR  SITE  POSITION 

C 1 5 3  A  *1  =  RADAR  SITE  CAN  NOT  SEE  THE  TARGET  BECAUSE  OF  THE  MAX.  RANGE 

C  163  A  0  »  RADAR  CAN  NOT  SEE  THE  TARGET  BECAUSE  OF  L INE-QF-S IGHT 

C 173  A  1  -  RADAR  CAN  SEE  THE  TARGET 

C  183  A 

C 193  A————————  MAIN  PROGRAM 

C203  'ENTER  THE  TARGET  ALTITUDE  (IN  FEET) ' 


C213  H«-0  4  H 

C223  A  H«- lOOO 

C233  'ENTER  THE  MATRIX  OF  THE  AVERAGE  TERRAIN  ALTITUDES' 

C243  ALT 1 -Q  A4  ALT 1 
CCS 3  AALTl-MM 

C263  'ENTER  THE  MATRIX  OF  THE  MINIMUM  ALTITUDES  OF  THE  RADAR  SITE' 
C273  ALT2*a  a 4  ALT2 
C283  a  ALT0*H2 

C293  I-J-l  »  CONf ( (pALTlCi 13) , (pALTlCi; 3) )pO 
C303  L3llfl 

C313  L2:ALT3-ALT1CI; J3-H  A4  'ALT1+H  -',*ALT3  4  I.J 

C323  ALT 4-ALT2C I ; J  3  A4  ' ALTSITE- ' ,»ALT4 

C333  -»<  (ALTCCIs  J3)**10)  /L7 

C34  3  CONCI; J3**10  4  -L4 

CCS 3  L7: - < ALT4P " 1 ) /Lt 

C363  CONC  I i J3-* l  4  -L4 

C373  Hi  —  (ALT 4<  ALT3 )  /L6 

C383  C0NCI;J3-0  4  -L4 

C393  L6: CONC I ;  J  3  - 1 

C 40 3  L4: 1-1*1 

C  4 1 3  -(IS (pALTICi 13) ) /L2 

C423  J-J+l 

C  43  3  - ( JS ( p ALT 1 C 1 ; 3 ) ) /L3 

C 44  3  THE  OUTPUT  I Si  ' 

C4S3  -  • 

C  46  3 

C  47  3  CON 

C483  THE  OUTPUT  IS  STORED  UNDER  THE  VARIABLE  CON' 


APPENDIX  D 


Appendix  D  describes  the  computer  function  named 
"MOVEMENT"  and  it  has  two  functions:  CPM  and  ERROR. 

A.  CPM 

HEADER  :  The  header  of  the  function  CPM  is  "CPM". 

PURPOSE:  This  function  calculates  the  time  steps,  the 
position,  the  altitude  over  this  position  and 
the  direction  of  all  flight  legs  of  all  aircraft 
for  each  time  step. 

INPUTS  :  The  computer  function  "MOVEMENT"  prompts  the 
following  inputs: 

1.  Time  Step:  The  time  of  the  time  step  of  model 
(in  seconds)  i.e.  if  the  calculation  has  to 
be  done  every  10  seconds,  the  input  for  time 
step  must  be  10. 

2.  Scale:  The  scale  used  for  the  terrain 
simulation,  i.e.  the  length  of  the  side  of 
the  terrain  square  grid. 

3.  Alt.  Data:  The  average  altitudes  given  by  the 
terrain  model  (in  feet). 

4.  X0 ,  Y0,  T0:  Cordinates  and  time  for  the 
beginning  of  the  flight  path  (in  NM  for 
cordinates  and  minutes  for  time). 

5.  Dir:  The  direction  of  the  flight  path  of  the 


aircraft  (in  degrees) 


6.  Length t  The  length  of  the  flight  leg  of  the 
aircraft  (in  NM) • 

7.  Speed:  The  speed  of  the  aircraft  along  the 
flight  path  (in  NM) . 

8*  Ait:  The  altitude  of  the  aircraft  along  the 
flight  path  (in  feet). 

OUTPUT  :  The  output  is  a  (5,N)  matrix  whose  elements 

are  as  follows: 

1.  First  Row:  Time  steps  of  the  air  combat 
model. 

2.  Second  and  Third  Row:  The  Cartesian 
coordinates  of  the  position  of  all  aircraft 
for  all  time  steps. 

3.  Fourth  Row:  The  altitude 

4.  Fifth  Row:  The  course  of  each  flight  leg. 

5.  N:  The  number  of  the  time  steps. 

The  output  is  stored  under  the  variable  A1 . 

B.  ERROR 

This  function  is  used  in  all  the  models  and  informs  the 
user  about  a  wrong  input. 


7  CPM;DPTI;TSPGDs ALT; DATA; DIR; DISTANCE: I : Cl s QC; SCALE; SPEED; TI ;  X ;  x 1 ; f 1 ;  XO; XYA 
S  YO  s  Y  s  TO  s  T  l ;  V ;  V 1 ;  V2  s Tl ; Z; Z1 s VECTOR; XYT ; FDATA1 
m  a 

CO  3  a - PURPOSE - - 

CC3  A 

C 4 3  A  THE  PURPOSE  OF  THIS  PROGRAM  IS  TO  CALCULATE  THE  TIME  STEPS, 

C33  A  POSITION ( X , Y> ,  THE  ALTITUDE  OF  THE  A/C  FOR  THE  PENETRATION 

t 63  A  ROUTE  ACROSS  THE  HOSTILE  AREA,  THE  COURSE  AND  THE  SPEED  FOR 

C73  A  EACH  FLIGHT  LEG  FOR  ALL  THE  TIME  STEPS. 

C33  A 
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C91  A - INPUT - - 

C10I  A 

C1U  A  1.  TIME  STEP  -  THE  TIME  OF  THE  TIME-STEP  MODEL  YOU  LIKE  (IN  SECOND) 

CIO!  a  I.E  IF  YOU  LIKE  THE  CALCULATION  TO  BE  DONE  EVERY 

C 130  A  10  SECOND  ENTER  FOR  TIME  STEP:  10 

C  143  * 

C 15 3  0.  SCALE  »  THE  SCALE  USED  FOR  THE  TERRAIN  SIMULATION.  I.E  THE 

C 161  A  LENGTH  OF  THE  SIDE  OF  THE  TERRAIN  SQUARE  GRID. 

C 171  A  (IN  N.M) 

C183  A 

C 191  A  3.  ALT. DATA  -  THE  AVERAGE  ALTITUDES  GIVEN  BY  THE  TERRAIN 
COO 3  A  MODEL.  (IN  FEET) 

C01I  A 

C003  A  4.  XO , YO , TO  -  CORDINATES  AND  TIME  FOR  THE  FLIGHT  PATH  BEGGINING. 

C03:  A  (  IN  N.M  FOR  CORDINATES  AND  IN  MINUTES  FOR  TIME) 

C04I  A 

C03I  A  3.  DIR  -  THE  DIRECTION  OF  THE  FLIGHT  PATH  (IN  DEGREE). 

C26I  A 

C07I  A  A.  LENGTH  =  THE  LENGTH  OF  THE  FLIGHT  PATH  (IN  N.N). 

C09I  A 

C29I  A  7.  SPEED  -  THE  SPEED  OF  THE  A/C  ACROSS  THE  FLIGHT  PATH  (IN  N.M) 

C30I  A 

C31I  A  3.  ALT.  =  THE  ALTITUDE  OF  THE  A/C  ACROSS  THE  FLIGHT  PATH 

C30I  A  (IN  FEET) 

C33:  A 

C341  A - OUTPUT - - 

C33:  A 

C36I  A  THE  OUTPUT  IS  A  (4  N)  MATRIX  WHOSE  ROWS  AND  N  ARE  AS  FOLLOWS: 

C37I  A 

C38I  A  1.  THE  FIRST  ROW  IS  THE  TIME  STEPS 

C39I  A  0.  THE  SECOND  AND  THE  THIRD  ROWS  ARE  THE  A/C  POSITION 

C40I  A  IN  X,Y  CORDINATES, 

C 4 1 1  A  3.  THE  FOURTH  ROW  IS  THE  A/C  ALTITUDE  ACROSS  THE  FLIGHT  PATH 

C40I  A  4.  THE  FIFTH  ROW  IS  THE  HEADING  OF  THE  FLIGHT  PATH 

C43I  A  3. THE  SIXTH  ROW  IS  THE  SPEED  OF  A/C  FOR  EACH  FLIGHT  LEG. 

C441  A 

C43I  A - OUTPUT  FOR  CHECK - - 

C46I  A  FLIGHT  DATA  ARE  STORED  UNDER  THE  VARIABLE  FDATA 
C47I  A  AS  FOLLOWS: 

C 491  A  FIRST  ROW  *  DIRECTION  FOR  EACH  FLIGHT  SEGMENT 

C 49 1  a  SECOND  ROW  »  LENGTH  OF  EACH  FLIGHT  SEGMENT 

C30I  A  THIRD  ROW  »  SPEED  OF  A/C  FOR  EACH  FLIGHT  SEGMENT 

C31 I  A  FOURTH  ROW  =■  ALTITUDES  OF  A/C  FOR  EACH  FLIGHT  SEGMENT 

C32I  A 

C33I  A - - - PROGRAM - * 

C34:  A 

C35:  V2-0  4  Al-  6  0  flO  4  FDATAl-  4  l  pO  4  FDAT A<-  4  0  pO 

C3ol  A - ENTER  DATA - - 

C57I  'ENTER  THE  TIME  STEP  (IN  SECOND) ’ 

C33 I  TI-C  4  ' ' 

C5?:  ENTER  THE  AVERAGE  ALTITUDES  FROM  TERRAIN  SIMULATION  ■ TDATA) ' 

C60I  QATA-Q  4 

C61I  ’ENTER  THE  SCALE  OF  THE  TERRAIN  SIMULATION  GRID  •'IN  N  M)  ' 

C62I  SCALE-0  4  '  - 

C63I  LC:  ' - ENTER  FLIGHT  DATA - ■ 

C64] 


CAS I  L7 : ' IF  THE  FLIGHT  PATH  IS  A  NEW  ONE  ENTER  0  ' 

C66I  ’IF  THE  FLIGHT  PATH  CONTINIUS  ENTER  1  ■ 

C67I  Ql-0 

C68I  -( (Ql»0> , (Ql=l) , ( Q 1 >1 ) )/L3,L6,L5 
C69I  LO: 'ENTER  (  XO.YO.TO  )  ' 

C703  XYT-0 

C71 1  -<  (pXYT>-3>  /L12 

C72I  ERROR  4  -LO 

C73I  L12: XO-XYTE 1 1  4  Y0-XYTC2I  4  T0-XYTC3I  4  T1+ <T0x60> +TI 
C7a:  Zl-(6,  (piTl)  )pO  4  Z 1 C  1 ;  3«-  <  V2*  (iTl  )->  4  V2-T1+V2  4  -LIS 

C73:  -L6  4  ’  ' 

C  76 1  L3: ERROR  4  ’  '  4  -L7 

C77I  L6: '  '  4  Zl-  6  0  AO 


C783  L13: 'ENTER  THE  DIRECTION  OF  THE  FLIGHT  PATH  (IN  DEGREE) ' 

C79 I  DIR-0  4  '  '  4  FDATAlCls 3-OIR 
C803  -(DIRS360)  /L20 

can  ERROR  4  -L13 

CS2I  L20s  'ENTER  THE  LENGTH  OF  THE  LEG  OF  THE  FLIGHT  PATH  (IN  N.M) ' 

Ca3I  DISTANCE-0  4  '  '  4  FDATA1  CDs  3-DISTANCE 
C34  3  L21 : 'ENTER  THE  SPEED  OF  A/C  (IN  KNOTS) ' 

C33I  SPEED-0  4  '  '  4  FDATA1C3; 3-SPEED 

C363  'ENTER  THE  FLIGHT  ALTITUDE  OF  A/C  (IN  FEET) ' 

C373  ALT-0  4  '  '  4  FDATA1 C4j 3-ALT  4  FDATA-FDATA.FDATAl 

C883  A - * 

C893  DPTI- (SPEED+3600) xTI  A (DISTANCE  PER  TIME  INCREMENT) 

C90  3  TSPGD-DISTANCE-DPTI  A (TIME  STEPS  PER  GIVEN  DISTANCE) 

C913  VECTOR- U (LTSPGD) )  *OPT I  APER  TIME  INCREMENT) 

C923  V— V2— t. (L ("It VECTOR) )  4  Vl-V  4  V-*ltV 

C933  A - - 

C943  AFIND  X  AND  Y  CORDINATES 

C933  A - - 

C963  X-XO-VECTOR:<(lo(  (  (on-iSO)xDIR)  )  A4  CX 
C973  Y-Y0-VECT0R:<(2o( ( (Ol)-1BO)XDIR> )  A4  TY 

C983  Xl-r (X-SCALE)  A4  Xl-CORD.IN  GRID 
C993  Yl— T (Y— SCALE)  A4  Yl-CORD. IN  GRID 

C  1 00  3  A - - 

C 101 3  A  CREATE  THE  T I ME-POS ITI ON-ALT I TUDE- D I RECT I QN-SF EED  MATRIX 

C  1 02  3  A - * 

C 1033  AA1— <6, (Tl ) )pO  4  AlClj  3-vTl 
C  104 3  1-1  4  XYA— (6,  (/»  X 1  >  >  pO 

C 1053  L1:XYAC; I  3- ( L VI C I  3 ) , (T XC I  3 )  , (T YC I  3 ) , (DATACX 1 C I  3 ; Yl C I  3 3+ALT) , DIR, SPEED 
C 1063  I-I+l 
C  1073  -(IS  (pXl)  )  /LI 

C 1 08 3  Al— A1 , Z1 , XYA  4  XO-'ltX  4  YO-"ltY  4  V2-" It (Al C 1 5 3 > 

C 1 09 3  LlOs  ■ IF  YOU  LIKE  TO  CONTINUE  ENTER  0  * 

Cl 103  OTHERWISE  ENTER  1  * 

C 1 1 1 3  Q2-0 

C 1 123  -( (02-0) , (02-1) , <Q2>1 > ) /L2,L8,L9 

Cl  133  L9: 'WRONG  ENTRY...  TRY  AGAIN'  4  -L10 
C 1 143  L8: 

C1153  A  — (T0*0) /LI 1 
C 1 163  A  Al-  0  1  4A1 
C 1 173  '  '  4  Al  4  '  '  4  '  ' 

C  1 133 
C  1193 
C  1203 


THE  OUTPUT  IS  STORED  UNDER  THE  VARIABLE  Al ' 


APPENDIX  E 


The  purpose  of  this  Appendix  is  to  describe  the  work 
spaces  Of  the  "TARGET  DETECTION",  "TARGET  INTERCEPTION", 
and  "PRIORITY"  models. 


A.  TARGET  DETECTION  MODEL 

The  TARGET  DETECTION  work  space  is  called  “FDETECT"  and 
it  consists  of  two  functions:  DETECT  and  ERROR. 

1 .  Detect 


PURPOSE:  The  purpose  of  the  function  "DETECT"  is  to 
find  out  if  a  particular  site  can  detect 
targets  which  are  flying  over  the  hostile 
terrain  in  specified  altitudes  and  routes. 

INPUT:  This  function  prompts  two  inputs: 

1.  The  minimum  altitudes  of  the  radar  site 
given  by  the  "MINALT"  computer  work 
space. 

2.  The  flight  data  of  the  targets  given  by 
the  "MOVEMENT"  computer  work  space. 

OUTPUT:  The  output  of  this  function  is  a  (M.N+l) 
matrix  where: 

K  •  The  number  of  time  steps 
The  first  column  is  the  time  step 
N  *  The  number  of  target. 


2.  Error 


See  Appendix  D 
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"TARGET  INTERCEPTION"  MODEL 

The  "TARGET  INTERCEPTION"  work  space  is  called  "INTERCPT" 


first,  the  target  2  second  and  so  on 
until  all  targets  are  entered  in  the 
program. 

OUTPUT:  The  output  of  this  program  is  a  (N,  1+(M*3)) 

matrix,  where: 

K  *  The  number  of  time  steps 

The  first  column  is  the  time  steps  from  1  to 

K 

N  ■  The  number  of  targets 
Each  target  has  two  columns: 

First  column  “1,  if  the  interception  is 
possible. 

0,  if  the  interception  is  not  possible. 
Second  column  **  The  number  of  time  steps 
that  a  possible  interception  takes, 
otherwise  zero. 

Third  column  ■  The  distance  between  radar 
site  and  target. 

C.  "TARGET  PRIORITY"  MODEL 

The  TARGET  PRIORITY  model  computer  workspace  is  called 
"PRIORITY"  and  it  consists  of  one  function:  PRIORITY. 

1 .  Priority 

PURPOSE:  The  purpose  of  the  function  "PRIORITY"  is 
to  calculate  the  priority  of  each  target 
when  a  radar  site  detects  more  than  one 
target  simultaneously. 

This  function  prompts  two  inputs: 
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INPUT: 


1.  The  critical  course  to  which  the  target 
flight  courses  will  be  compared. 

2.  The  target  flight  data  given  by  the 
"MOVEMENT"  computer  program. 

OUTPUT:  The  output  of  this  function  is  a  (m,  1+N) 

matrix,  where: 

K  -  The  number  of  time  steps  from  0  to  M 
The  first  column  is  the  time  steps 
N  *  The  number  of  targets. 

The  values  in  the  matrix  represent 
priority  number. 


v  DETECT ; A; A2; X:  SX  ; Y; SY; TGT;  ZERO; M; D: Ki Rs I ; A t  s ALT; KK; HALT; MALT! ; P 
C  1  3  A 

CC3  A - PURPOSE - * 

C33  a 

14  3  A  THIS  PROGRAM  IS  USED  TO  FIND  OUT  IF  A  PARTICULAR  RADAR  SITE 
CS3  *  CAN  'SEE'  TARGETS  WHICH  ARE  FLYING  OVER  THE  HOSTILE  TERRAIN 
C63  A  IN  SPECIFIED  ALTITUDES  AND  ROUTES. 

C73  A 

ca:  a -  inputs - - 

C93  A 

C 103  A  THIS  PROGRAM  PROMPTS  TWO  INPUTS; 

C  1 1  3  A 

C1C3  A  1.  THE  MINIMUM  ALTITUDES  OF  THE  RADAR  SITE  GIVEN  BY  THE  ' MINALT ’ 
C 133  A  PROGRAM. 

C143  A  2.  THE  FLIGHT  DATA  OF  THE  A/C  ( TARGET 1  OR  2  ETC)  GIVEN  BY  THE 
C133  A  PROGRAM  MOVEMENT'. 

C  163  A 

C  173  A - OUTPUTS - * 

C 183  A  THE  OUTPUT  OF  THIS  PROGRAM  IS  A  (M  N>  MATRIX  WHERE: 

C  193  A 

C CO 3  A  M  »  THE  NUMBER  OF  TIME  STEPS 

CC1  3  A  N  *  THE  NUMBER  OF  TARGETS 

CCC3  A 


CC33  A - MAIN  PROGRAM - - 

CC43  A 

CCS  3  D-A  1-1.0 

CC63  ’ENTER  THE  MINIMUM  ALTITUDES  OF  THE  RADAR  SITE  (  MALTISl.C,  ETC  >’ 
CC73  MALT I -0 


CCS 3  L4;  ENTER  THE  TARGET  FLIGHT  DATA  (  TARGET 1 ,2, ETC  )’ 

CC93  TGT-0 

C303  -<  (pTGTC; 1 3) «6> /L3 

CCl 3  ERROR  9  -L4 

C3C3  L3:K>0  0  D-D , (pTGTCC; 3 )  9  SX-f (TGTC2; 3-6)  *  SY-C (TGTC3; 3-6)  9  P-pSX 
C333  L10: K-K-l 

C343  ALT-TGTC4; 3CK3  9  X-SXCK3 

CCS 3  — ( X^O) /L12  A  CHECK  IF  THE  TARGET  HAS  ENTERED  INTO  THE  HOSTILE  AREA 


C363  Al-Al.O  9  *L1Z 

C373  LlC:MALT<>-MALTICX<i-SXCK3;  Y*SYCK33  AO  X  0  Y  9  MALT  9  ALT 

C383  ♦ (MALT^'t  > /L9  A  CHECK  IF  THE  TGT  IS  INSIDE  THE  MAX  RANGE  OF  SITE. 

C393  A1*A1,0  ♦  ^L13 

C 403  L9:-*( ALT? MALT) /Lit  A-  CHECK  IF  THE  ALT  OF  THE  TST  IS  HIGHER  THAN 

C41  3  A - THE  RADAR  SITE  ALTITUDE. 

C 403  A1  *-Al ,  0  9  -»L13 
C A3 3  LL1 : Al-Al , 1 
C443  L13s  •*  (K<P)  /L10 
C453 

C463  'IF  YOU  HAVE  AN  ANOTHER  TARGET  ENTER  O' 

C 473  'IF  YOU  LIKE  TO  EXIT  ENTER  1' 

C483  KK-Q 

C493  *<KK-0,KK»1 ,KK>2> /L4,L7,L6 

C30  3  L6: ERROR  9  *L7 
C313  L7;A2<M  (pD>  ,  (M*lt  (DCtD3  )  )  )  AO 
C323  I<-1 

C333  L8:  A2C  I }  3*  <  (DC  I  3t  A1 )  ,  (M-DC  I  3  >  ftO)  9  A1«-DCI3  +  A1  9  1*1*1 

C343  * ( li (pD) > /LS 

CSS 3  A3M  (  <M,  1)(M«)  ,OA2) 

C363  A3 
CS73  •  '  9  ■  ' 

C383  - - 

[593  '  THE  OUTPUT  IS  STORED  UNDER  THE  VARIABLE  A3' 

C  60  3  ' - - - ' 


W;W2; 

xo 

C13 

H  1 

C23 

a 

C33 

a 

C  4  3 

a 

C33 

a  — 

C63 

a 

C73 

a 

C83 

a 

C93 

a 

C  103 

a 

cm 

a 

C  123 

a 

C  133 

a 

C  143 

a 

C  133 

a 

C  163 

a 

C  173 

a 

C  183 

a 

C  193 

a 

C20  3 

a 

C213 

a 

C223 

a 

C233 

a 

C243 

a 

C233 

a 

C263 

a 

C273 

a 

C233 

a 

7  INTER!  A:  ANS;  D; DF; DF1 ; DIR:  I ;  K;  L:  M;  MM*.  Q;  RMAX  i  RO;  R01 !  SPEED:  T;  TGT;  TI ;  VI ;  V! 
; XR: XS; XT; AC; TEST; AA;OIST 

- PURPOSE - * 


THIS  PROGRAM  IS  USED  TO  CALCULATE  THE  INTERCEPTION  STATE  MATRIX 
FOR  EACH  RADAR  SITE 

- INPUT - - 


A  2.  THE  LOCATION  <X,Y)  OF  THE  MISSILE  RADAR  SITE  (IN  N.M). 


THE  DETECT  FILE  (DETECTS l  OR  2  ETC)  GIVEN  BY  THE  ’ FDETECT ' 
WORKSPACE.  THIS  RADAR  SITE  IS  SUPPOSED  TO  SUPPORT  THE 
INTERCEPTOR  MISSILE  OR  A/C. 


3.  THE  SPEED  OF  THE  INTERCEPTOR  MISSILE  (IN  KNOTS) 


' MOVEMENT '  WORKSPACE . 


NOTE:  THE  ORDER  OF  ENTERING  THE  FLIGHT  DATA  OF  THE  TARGETS 

-  MUST  BE  IN  ASCENT  ORDER.  I.Es  YOU  MUST  ENTER  THE 

TARGET 1  FIRST  THE  TARGET2  SECOND  AND  SO  ON  UNTIL 
ALL  TARGETS  ARE  ENTERED  IN  THE  PROGRAM. 
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C36]  Al+Al.O  4  +L13 

C37]  L10:  MALT+MALTI  CX+SXCK] ;  Y+SYCK]  3  44  X  *  Y  4  MALT  4  ALT 

C38]  -(MALTP'l)  /L9  4  CHECK  IF  THE  TGT  IS  INSIDE  THE  MAX  RANGE  OF  SITE. 

C39]  A1*A1,0  4  -!>L13 

C40]  L9: -MALT)  MALT) /Lll  A-  CHECK  IF  THE  ALT  QF  THE  TGT  IS  HIGHER  THAN 

C41]  4 - THE  RADAR  SITE  ALTITUDE. 

C40]  Al+Al.O  4  -L13 
C43]  LL1:A1-At,l 
C44]  L13:  *  (K<P)  /LIO 
C  43 1 

C46]  ’IF  YOU  HAVE  AN  ANOTHER  TARGET  ENTER  O’ 

C47]  'IF  YOU  LIKE  TO  EXIT  ENTER  f 

C4S]  KK+O 

C49]  *<KK«0,KK«l,KK>0> /L4,L7,L6 

C30]  L6: ERROR  4  -L7 
C31]  L7 : AC* ( (flD) , (M+lt <DCtD] > > ) pO 
C50]  1*1 

C33]  L8;A2CI; ]*( (DCIJtAl) , (M-DCI3) pO)  4  A1*DCI]4A1  4  1*1+1 

C54]  *(I< (pD) > /L9 

C33]  A3*  ( ( (M , 1 )  pt,M>  ,  HAD) 

C361  A3 
C57]  ’  '  4  '  ■ 

CSS]  ’ - ’ 

C391  THE  OUTPUT  IS  STORED  UNDER  THE  VARIABLE  A3’ 

C  60 ]  ’ - - - ' 


7 

W:WO: 

cn 

CO] 

C3] 

C4] 

C3] 

C  6  ] 
C7] 
C8] 
C9] 
CIO] 

C  1 1  ] 
CIO] 

C  13] 

C  14] 
CIS] 

C  16] 

C  17] 
C10] 

C  19] 
COO] 
C01] 
COO] 
C03] 
C04] 
COS] 
C06] 
C07] 
C03] 


INTER: A; ANS: D: DF: DF1 : DIR; I : K: L: M; MM: Q: RMAX ; RO: R01 : SPEED: T: TGT:  TI ; VI ; VO: VR: 
XO: X 1 : XR; XS; XT ; AC: TEST ; AA: DIST 

4  4 - PURPOSE - - 

4 

4  THIS  PROGRAM  IS  USED  TO  CALCULATE  THE  INTERCEPTION  STATE  MATRIX 
4  FOR  EACH  RADAR  SITE  . 

4 - INPUT - - 

4 

4  THIS  PROGRAM  PROMPTS  THE  FOLLOWING  INPUTS: 

4 

4  1.  THE  TIME  OF  THE  TIME  STEP  OF  THE  SIMULATION  (IN  SEC). 


0.  THE  LOCATION  <X,Y)  OF  THE  MISSILE  RADAR  SITE  (IN  N.M) 


4  6. 
4 


THE  DETECT  FILE  ( DETECTS 1  OR  0  ETC)  GIVEN  BY  THE  ' FDETECT ' 
WORKSPACE.  THIS  RADAR  SITE  IS  SUPPOSED  TO  SUPPORT  THE 
INTERCEPTOR  MISSILE  OR  A/C. 


4  4. 
4 


THE  MAXIMUM  RANGE  OF  THE  INTERCEPTOR  MISSILE  (IN  N.M) 
THE  SPEED  OF  THE  INTERCEPTOR  MISSILE  (IN  KNOTS) 


THE  FLIGHT  TARGET  DATA  (TARGET  1  OR  0  ETC)  GIVEN  BY  THE 
’  MOVEMENT ’  WORKSPACE. 


NOTE:  THE  ORDER  OF  ENTERING  THE  FLIGHT  DATA  OF  THE  TARGETS 

-  MUST  BE  IN  ASCENT  ORDER.  I.E:  YOU  MUST  ENTER  THE 

TARGET l  FIRST  THE  TARGETO  SECOND  AND  SO  ON  UNTIL 
ALL  TARGETS  ARE  ENTERED  IN  THE  PROGRAM. 
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C9&3  A -  ANGLE  BETWEEN  TARGET  AND  INTERCEPTOR  -  0  DEGREE 

C97 J  L3:VD-Vt*V2 

C981  A  TIME  FOR  INTERCEPTION  (IN  MIN)  »  '  ,  *  T- ( D+VD>  '*60 
C99I  A  DISTANCE  FOR  TARGET  -  ',lT«(Vl-60> 

C 1003  A  -DISTANCE  FOR  INTERCEPTOR  =  '  ,  *  T  c  (V2-S-60) 

CIO!!  Q-T:<  (V2-J-60) 

C  1021  L l O :  Q  P  ■»  (QiRMAX  )  /LI  1 
C  1033  A— A ,0,0  P  D I ST-D  1ST , 0  P  -LP1 

C1043  H1:DIST-DIST, (L D)  A  CREATES  THE  DISTANCE  VECTOR 
C  1033  AfA.l.r  (  (T*60>-5-TI>  P  -LP1 
C  1 06 3  END: Al-( (pDF) ,2)pA  P  D1ST-  (  (pDF)  ,  l  JpDIST 
C 1073  AA-AA,A1,DIST  P  Al-AA 

C 1003  ’IF  YOU  HAVE  ANOTHER  TARGET  ENTER  O' 

C 1 09 3  'TO  EXIT  ENTER  f 

Cl  10  3  ANS-0 
C  1 1 1  3  -(ANS-O) /L14 

C 1 123  • - 

Cll-~3  the  OUTPUT  IS  STORED  UNDER  THE  VARIABLE  Al ' 

C  1 1 4  3 
V 


V  PRIORITY: TGT 1 j TGT2 ; AN; A ; A 1 ; K; T ; I ; P : p  1  j  D I R 1 : D IRC ; D IR3 ; ALT  1 ;  V; M: TGT ; DEGREE 
C  1  3  A 

C23  A - PURPOSE - - 

C33  A 

C43  A  THIS  PROGRAM  IS  USED  TO  FIND  OUT  THE  PRIORITY  OF  THE  TARGETS 
C33  A  TO  BE  KILLED  BY  THE  ENEMY  RADAR  SITES,  AS  FAR  AS,  HOW  DENGEROUS 
C63  A  THE  TARGET  IS,  IS  CONCERNED. 

C73  A 

C03  A -  INPUTS - * 

C93  A 

C 103  A  THIS  PROGRAM  PROMPTS  THREE  INPUTS: 

C113  A 

C 123  A  1.  THE  CRITICAL  COURSE  THAT  THE  TARGET  FLIGHT  COURSES  WILL  BE 
Cl 03  A  COMPARED. 


C  1 4  3  A 

C 133  A  2.  THE  TARGET  FLIGHT  DATA  (TARGET 1 , 2  ETC)  GIVEN  BY  THE  MOVEMENT' 
C 163  A  PROGRAM. 


C  173  A 

C  193  A - -  OUTPUT  - - 

C  193  A 

C203  A  THE  OUTPUT  IS  A  (M  N+l)  MATRIX  WHERE: 


C213  A 

C223  A  M  «  THE  NUMBER  OF  TIME  STEPS 

C23  3  A  N  -  THE  NUMPER  OF  TARGETS 

C243  A  THE  FIRST  COLUMN  IS  THE  TIME  STEPS 

C233  A 

C 26 3  A - MAIN  PROGRAM - 

C27  3  A 

C293  TGT 1 -TGT2-D IR 1-D IR3-ALT l »V*vO  P  P-0  - 
C 29 3  ENTER  THE  CRITICAL  COURSE' 

CTO 3  DEGREE-0  P  ' ' 

C313  LI: 'ENTER  THE  TARGET  FLIGHT  DATA  (TARGET 1 , 2 , ETC) ' 
C323  TGT -0  P  P-P-l  P  V-V ,pTGTC 1 : 3 
C333  TGT 1 -TGT 1 , TGTC4; 3 
C  34  3  TGT2-TGT2 , TGTC3: 3 
C333  K-p (TGTC  1:3) 

C363  L4:  ' IF  YOU  HAVE  ANOTHER  TARGET  ENTER  O' 

C373  '  OTHERWISE  ENTER  l' 


C383  AN-0 

C393  -< (AN-O) , (AN-1) , (ANiC) >/Ll,L2,L3 

C  40  3  L3: ERROR  P  -L4 

C  4 1 3  L2 : A— A 1  —  ( P ,  ( M- (It (VC  TV 3  )  )  >  >  pO  P  1-1 

C  42  3  L6: ACI:  3— ( ( VC  I  3 ) t  TGT l )  ,  ( M— VC  I  3 ) pO  P  TGT l - ( VC  I  3 ) + TGT 1 
C  43  3  AlCIs 3-( (VCI3)tTGT2>  , (M-VCI3>pO  P  TGT2- ( VC  I  3 ) PTGT2  P  I-I  +  l 
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THE  OUTPUT  IS  STORED  UNDER  THE  VARIABLE  A: 


APPENDIX  F 


This  Appendix  contains  all  the  computer  programs  of  the 
"AIRMODEL*  workspace. 


>7  MA  INS  CMD!  WTH:  STEP;  CAPPC; TK1 : ADSITES; AISITES: NSITES; SDD2: F; MCAPPC:  II;  12;  IC; 
Rl;RC;RC; T; TIMES 
Cl!  A 

CC3  4 - PURPOSE - - 

CC3  ft 

C 4 3  4  THE  PURPOSE  Of  MAIN  PROGRAM  IS  TO  COMBINE  THE  DATA  BASE  WITH 

CS3  ft  the  functions  used  to  calculate  the  attrition  rate  qf  the 

C6  3  ft  AIRCRAFT  FORMATION. 

C73  ft 

C33  4 -  INPUT - - - - 

C93  ft 

C 10  0  ft  THE  MAIN  PROGRAM  PROMPTS  THE  FOLLWING  INPUTS: 

Ill]  ft 

C  1C3  ft  1.  RL:  THE  SEED  FOR  THE  RANDOM  NUMBER  FUNCTION. 

C1C3  4  C.  TIMES:  THE  NUMBER  OF  TIMES  THAT  THE  PROGRAM  WILL  RUN. 

C 143  4  C.  NSITES:  THE  NUMBER  OF  THE  RADAR  SITES. 

C133  ft  4.  DSC:  THE  DETECTION  MATRIX  OF  EACH  RADAR  SITE  (GIVEN  BY 

C 163  4  THE  FDETECT  PROGRAM). 

C 173  4  3.  ISC:  THE  INTERCEPTION  MATRIX  OF  EACH  RADAR  SITE  (GIVEN 

C  183  4  .BY  THE  INTERCPT  PROGRAM). 

C 193  4  6.  TGT :  THE  NUMBER  OF  AIRCRAFT  THAT  EACH  FORMATION  HAS. 

CC03  4  7.  CAPC:  THE  NUMBER  OF  TARGET  THAT  EACH  RADAR  SITE  CAN 

C213  4  ENGAGE. 

CCC3  4  a.  MCAFC:  THE  NUMBER  OF  MISSILES  EACH  RADAR  SITE  HAS. 

CCC3  4  9.  PKM:  THE  PROBABILITY  OF  KILL  MATRIX. 

CC43  4  10.  IFF:  THE  DECISION  IF  IFF  EXISTS. 

CCS 3  4  11.  IFFWORK: THE  PROBABILITY  THAT  IFF  WORKS. 

CC63  4  1C.  IFFFOF:  THE  DECISION  IF  THE  TARGET  IS  FRIEND  OR  FOE,  GIVEN 

CC73  4  THAT  IFF  DOES  NOT  WORK. 

CCB3  4  13.  IFFCOR:  THE  PROBABILITY  THAT  IDENTIFICATION  IS  CORECT, GIVEN 

CC93  4  THAT  IFF  WORKS. 

C303  4  14.  ECM:  THE  ELECTRONIC  COUNTERMEASURE  MEANS  MATRIX  FOR  THE 

C313  4  ECM  FUNCTION  GIVEN  BY  THE  ECM  MODEL. 

C3C3  4  IS.  CMD:  THE  DECISION  FOR  THE  COMMAND  AND  CONTROL  MODEL. 

C333  4  16.  AUTO:  THE  NUMBER  OF  AIRCRAFTS  ABOVE  OF  WHICH  THE  RADAR 

C 34 3  4  SITES  OPERATE  AUTONOMOUSLLY. 

CCS 3  4  17.  PR:  THE  PRIORITY  MATRIX  GIVEN  BY  THE  PRIORITY  MODEL. 

CC63  4  18.  WTH:  THE  DESITION  FOR  GOOD  OR  BAD  RADAR  SITE  OPERATION. 

CC7  3  4  19.  LBR:  THE  FACTOR  THAT  THE  PROBABILITY  OF  AQUASITION 

CCS 3  4  BREAKINK  WILL  INCREASE,  GIVEN  THAT  THE  OPERATION 

CC9 3  4  OF  THE  RADAR  SITES  IS  DISGRACED. 

C 40 3  4  CO.  ETSF:  THE  FACTOR  THAT  THE  EVENT  TIME  STEPS  MATRIX  WILL 

C 4 1 3  4  INCREASE  .GIVEN  THAT  THE  OPERATION  OF  RADAR  SITES 

C4C3  4  .  IS  DISGRADED. 

C 43 3  4  Cl.  TSTEP1:  THE  EVENT  TIME  STEPS  MATRIX. 

C  44  3  4 

C433  4 - OUTPUT - - 

C  46  3  4 

C473  4  THE  MAIN  PROGRAM  GIVES  THE  FOLLOWING  OUTPUTS: 

C483  4 

C 49 3  4  1.  A  SUMMARY  OF  THE  USED  SENARIO. 

CS03  4  C.  IN  DETAILS  ATTRITION  MATRIX. 

CS13  4  C.  REMAIN  MISSILES  PER  RADAR  SITE. 

CSC3  4  4.  THE  NUMBER  OF  TARGETS  WHICH  ARE  KILLED  BUT  NO  ASSESSED. 
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C33J  A  3.  THE  NUMBER  OF  KILLED  TARbfcfS. 

C341  a 

C331  * - COMMENTS - 

C3ol  A 

C371  A  THE  MAIN  PROGRAM  CALLS  THE  ROLLWING  PROGRAMS: 

C3s:  a 

C391  *  1.  SITE 

I  oO  I  A 

'.ail  A - MAIN  PROGRAM - 

Ia21  A 

C631  TGTSTAT-x.0 

C641  'ENTER  THE  NEW  SEED' 

C631  ORL-O 

C661  ENTER  THE  NUMBER  OF  TIMES  THE  PROGRAM  TO  RUN ' 

C671  TIMES-0 
C68I  T-0 

C691  LP3: - (TiTIMES) /TELOS 
C701  T  — T  ♦  1 

C711  ADS  I  TES-A  I S I  TES-s.  0 

C721  fl - RADAR  SITES  DATA - 

C731 

1741  ENTER  THE  NUMBER  OF  THE  RADAR  SITES' 

1731  NSITES-4  9  9  TESTMCAP-s.  NSITES 

C761  P-0 

C  77 1  PRO: P-P+ 1 

C781  - (P, NSITES) /PP1 

179]  ENTER  THE  DETECTION  MATRIX  OF  THE  SITES',  (*P> 

C801  DDS2-Q  9  ADSITES-ADSITES, ( ,DDS2)  9  -PP2 

can  PP1 :  DDS2-  (NSITES ,  (/>DDS2C;  II)  ,  (RDDS2C1;  I)  >pADSITES 

C821  STGT- ( (NSITES) , ( (pODS2Cls Is  3 ) — 1 > >p(v ( (pDDS2Cls Is  3 )  —  1 > ) 

CS31  TEST I NT- ( (NSITES) , ( (pODS2tl: Is  3>-l> )pl 

C341  P-0 

C8S1  PP3: P-P+l 

C361  —  (P'-NSITES)  /PP4 

C871  'ENTER  THE  INTERCEPTION  MATRIX  FOR  THE  SITES', (5P> 

C8S1  IS2-0  9  AISITES-AISITES, (,IS2)  4  -PP3 

C891  PP4 : MIS2— (NS ITES , <ptS2C: 11) , (RIS2C1; I) >pAISITES 

C90I  11-2  9  12-3  o  13-4  0  R1-R2-R3-V0 

C91 1  1 1 1 :  -  <  (pRl ) i  (  (pDDS2Cl;  Is  D-l) )  /1 12 

C92I  Rl-Rl , I t  9  R2-R2, 12  0  R3-R3,I3 

C93I  I1-I1+3  *  I2-I2-3  0  I3-I3+3  9  -til 

C94I  1 12: SEF-MIS2C ; s R1 1  »  STS-MIS2C ; s R2I  0  STD-MIS2C ; ; R31 
C93I  ER9s 'ENTER  THE  VECTOR  NUMBER  OF  A/C  THAT  EACH  TARGET  HAS' 

C96I  TGT-Q 

C97I  -( (pTGT>»( (pDDS2C l; 1; 1 ) — 1 > ) /ER8 
C98I  ERROR  o  -ER9 

C 99 1  ER8: TGT1-TGT  9  ALLTGT-v (pTGT) 

ClOOl  PP6:  ENTER  THE  NUMBER  OF  TGTS  THAT  THE  SITE  CAN  ENGAGE' 

ClOll  CAP2-0  9  CAPP2-CAP2  9  *' 

C  1021  —  (  (pCAF'2 )  "NSITES )  /FP7 

C 1 03 1  ERROR  «  -PP6 

Cl 041  PP7 : 'ENTER  THE  NUMBER  OF  MISSILE  THAT  EACH  SITE  HAS  (VECTOR) ' 
C1051  MCAP2-0  9  MCAPP2-MCAP2  9  ' 

C1061  -(  (pMCAF'2)  -NSITES)  /PP8 

C 1071  ERROR  9  -PP7 

C1081  PP9: 'ENTER  THE  PK  MATRIX  OF  THE  MISSILE  FUNCTION' 

C 1091  PKM-0  9  ' ' 

Cl  101  ER2s  ' IF  YOU  LIKE  IFF  IDIFICATION  ENTER  1' 

Clill  'OTHERWISE  ENTER  O' 

C 1 121  IFF- 1  9  '  ' 

Cl  131  -( (IFF-O) , (IFF-1) ) /ER3,ER10 

Cl  141  ERROR  9  -ER2 

Cl  131  ER1 Os  'ENTER  THE  PROBABILITY  THAT  THE  IFF  WORKS' 

Cl  161  IFFWORK-O  9  ' 

Cl  171  - ( IFFWORKPO) /DD1 

CU81  ERls'IF  YOU  LIKE  ALL  TGTS  TO  BE  FRIENDS  ( IFF  NOT  WORK)  ENTER  1' 
Cl  191  '  IF  YOU  LIKE  ALL  TGTS  TO  BE  FOE  ENTER  O' 

C 1201  IFFFOF-O  9  ' ' 

C 121 1  -< ( IFFFOF-O) , ( IFFFOF-l ) ) /ER3 
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C  1222  ERROR  ■>  -tKl 

Ci::i  DDl :  ENTER  THE  PROB.  THAT  THE  IDIFICATION  IS  CORRECT /IFF  WORKS) 

C124  2  IFFCOR-Q  ■> 

C  1252  ER3: 'ENTER  THE  ECM  MATRIX  OF  ECM  FUNCTION' 

C  1261  ECM+O  »  '  ' 

C 127]  TTls'IF  YOU  LIKE  COMMAND  CONTROL  MODEL  ENTER  1' 

Cl 29]  'OTHETWISE  ENTER  O' 

C  1 29  ]  CMD-0  ■)  '• 

C1302  -(CMD=0)/ER4 

C 131 ]  'ENTER  THE  NUMBER  OF  AIRCRAFTS  ABOVE  OF  WHICH  THE  RADAR' 

C 132]  'SITES  OPERATES  AUTONOMOUSLY' 

Cl 33]  AUTO+O 

C 134]  ER4:  ENTER  THE  PRIORITY  MATRIX' 

C 135]  PR-0 

C 1 36 ]  ER44s ’ IF  YOU  LIKE  GOOD  RADAR  SITE  OPPERATION  ENTER  1' 

C 137]  'IF  YOU  LIKE  DISGRADED  RADAR  SITE  OPPERATION  ENTER  O' 

C 139]  WTH+l  »  '' 

C 139]  -( (WTH-O) , (WTH-l) ) /ER5.ER55 

Cl 40]  ERROR  »  -ER44 

C 1 4 1 ]  ER5:  ENTER  THE  INCEASE  FACTOR  FOR  AQUASITION  BREAK  PROBABILITY' 

C 142]  LBR-M3  4  '  ' 

C 143]  'ENTER  THE  INCREASE  FACTOR  FOR  EVENTS  TIME  STEPS  MATRIX' 

C 144]  ETSF-0  4 

1 1 43 ]  TSTEP1-TSTEP-TSTEP1-M.TSTEP1;<(ETSF>  4  'DISGRADED  OPPERATION 
C  1 46  ]  ECMC  ;  4  2-ECMC  i  42+ECMC  ;  4  2:<L6R 

C  147]  ERS5-.  ENTER  THE  EVENTS  TIME  STEPS  MATRI X ( ID  .TRACK , FIRE ,ASE5S)  ' 

C  1 48 ]  TSTEP-0  »  '•  0  TSTEP 1 -TSTEP- ( (pTGT)  ,  (NSITES)  ,  4)  p  , TSTEP 
C 1 49 ]  DS2-DDS2  4  DS2=THE  DETECTION  MATRIX  OF  SITE  2  +  10  ROWS 
C 150]  ASS-0  »  TFERROR+1 

C151  ]  AS2-( (NSITES, < (pDS2C 1 ; ; 1 ] )  +  10) , 1 )  /»  ( v ( (pDS2Cl; j l]>+10) ) ) , < (NSITES, ( (pDS2C 1 

; ; 13) +10) , ipALLTGT) ) p 0) 

C1523  4 - TIME  STEP  LOOP - - 

C 1531  1+0 

C 154]  4  LP2  x  +  ( I M ) /END 
C 1553  LP2:  +  ( I J  ( <pOS2C 1 ; ; 1 ] ) ) ) /END 
C 156]  I-I-l  P  '  '  »  'I-'  ,*I 

C 157]  4 - RADAR  SITES  LOOP - - 

C 1591  SS-0 

C159]  LP1 :  ■*  (SS> NSITES)  /LP2 
C 1 60 ]  SS-S3+1 

C 16 1 ]  4  - ( MC AP2 CSS] JO) / LLL 1 
C 1621  4  -  ( I  >  TE3TMCAP  CSS]) /LP l 
C 163]  ' . . .SITE  «  ' ,*SS 

C 164]  SS  SITE  I 
C 165]  -LP 1 

Cl 66]  END:  4 - END  OF  THE  SIMULATION - * 

C 167]  SENARIQ' 

C  169]  - ' 

C 169 ]  'NUMBER  OF  TARGET  PER  FORMATION  =  ',*TGTl 
C 170]  'TARGET  CAPACITY  PER  MISSILE  SITE  -  ',(*CAPP2> 

C1712  'MISSILE  CAPACITY  PER  MISSILE  SITE*  ' , (5MCAPP2) 

C 172]  ' IFF  =  ' ,  (5  IFF) 

C 1 73 ]  'PROB  THAT  IFF  WORKS  =  ',(*IFFWORK> 


C  174] 

'PROB  THAT  IDENTIF.  IS  CORRECT/ IFF  WORKS 

■  '  ,  (5 IFFCOR) 

C1752 

'FRIEND/ENEMY  IF  IFF  DOES  NOT  WORKS  *  ', 

( * IFFFOF ) 

C  1762 

'COMMAND  CONTROL  MODEL  =  ',(*CMD> 

C  177] 

'RADAR  OPPERATION (GOOD /BAD)  =  ' , (SWTH) 

C  1732 

'ECM  MATRIX  »  ' 

C  1 79  2 

ECM 

C  1 30  ] 

'PK  MATRIX  =  ' 

C  131  2 

PKM 

C 1322  '  - RESULT - 

C  1332  4 

C 1 94 ]  - ( I FF»o  > / AS 1 

C  1352  TGT  2  HAS  BE  ENOEOED  =  '  .  .  s  ••  TRPRRne-- i  ■.  i 
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C1863  ASlsTKl-S(  (  (,oTK>-3>  ,3>pTK 
C 187]  IN  DETAILS  ATTRITION  MATRIX: 

C1881  'MISSILE  SITE  -  ’  ,  (STK1CI;  I) 

C 1893  TIME  STEP  -  ‘ , C*TK1 CO: 3 > 

C 1 90 3  'KILLED  TGT  -  -,<»TKIC3;3> 

C 191 3  TGTST AT -TGTST AT , ( TGT  1 -TGT) 

C 1933  -LP3 

C  1933  TELOS:TGTSTAT-o<T,  (pTGTl)  >  pTGTSTAT 
C 1943  REMAIN  MISSILES  »-,*MCAP2 

C 1933  NO  ASSESS  TGT  »  ',*ASS 

C 1963  ’KILLED  TGTS  (TGTSTAT) “  '  4  A ,(*< .TGTSTAT) ) 

C 1973  TGTSTAT 

C 1983  A - ERASE  GLOBAL  VARIABLES - 1 - * 

C 1993  QQ-OEX  ‘ STGT  TESTINT  ISC  SEF  STS  STD  TGT  TGTt  ALLTGT  CAP2  MCAP2  IFF  PR  TS 
TEP  TSTEP1  DS2  RD  RN  NT  PK  I  SS  INT  CHECK  ISITE  TK  TK1  ASS  ECM  TFERROR  TK1 
V 


9  SS  SITE  I ; EV; FF; ID; II;Ks  MT  j  NNs  NT ;  Z ; NT 1 ; P; PRIOR; TT ; Z ; EV: TEST ; AM: PK  C 


C  1 3  a 

C23  a - PURPOSE - * 

C33  a 

C43  A  THIS  FUNCTION  IS  USED  TO  CALCULATE  THE  ATTRITION  RATE  OF  THE 
C33  A  ATTACKING  A1CRAFT  FORMATIONS  AGAINST.  GROUND  TARGETS 
C63  A 

C73  A -  INPUT - - 

CS3  A 

C93  A  THIS  FUNCTION  IS  CALLED  BY  THE  MAIN  PROGRAM  AND  USES  THE  FOLLWING 
C 103  A  INPUTS: 

Cl  13  A  1.  THE  LEFT  ARGUMENT  IS  THE  NUMBER  OF  RARAR  SITE 

C 123  A  2.  THE  RIGHT  ARGUMENT  IS  THE  TIME  STEP 

C 1 3 3  A  3.  SOME  OTHER  GLOBAL  VARIABLES 
C  143  A 

C  133  A - OUTPUT - - 


C 163  A 

C 173  A  THIS  FUNCTION  CALLS  THE  IDIF,  TRACK,  FIRE,  ASSES, INTER  .RANDOM 
C 183  A  AND  PK I LL  FUNCTIONS  IN  ORDER  TO  PERFORM  THE  FUTURE  EVENTS  AND 
C 193  A  CALCULATE  THE  ATRITION  RATE, 


C203  A 

C21 3  a - MAIN  PROGRAM 

C223  A 

C233  A -  INITIALIZE  - 


C24  3  -<I>1)/C0NT 

C233  NT1-0  4  TK*\.0  ATIME/KILL  VECTOR 

C263  CONT: ISITE-I ,SS  ACONTAINS  THE  TIME  STEP  *  THE  NUMBER  OF  SITE  CALLED 

C273  -< (+/NN-UAS2CSS; I; 3>«0> /L4  aCHECK  IF  THERE  IS  AN  EVENT 

C233  P-0  A4  'THERE  IS  EVENTS’  4  -NN-',*NN 

C293  LPl:-<P>  </>NN)  )  /L7  A 4  LOOP  FOR  EVENTS 

C303  P-P+l  A 4  -P»-,*P 

C31 3  EV-NNCP3 

C323  - <EV=0) /LPl  a  THE  END  OF  THE  EVENTS  A  EVENT=-,*EV 

C333  EV-EV-4 (*3S) , -00'  A  FORM  THE  EVENTS.  4  'EV-',*EV 

C34  3  -<  <EV< 19)  ,  (EVS29) ,  (EVS39) , (EVS 49) , (EVS59) > /EV1 ,EV2 ,EV3 ,EV4 ,EV5 

C333  EVl sNT-EV-10  4  TT-ISITE  TRACK  NT  4  -LPl  A  CALL  THE  TRACK  FUNCTION 

C363  EV2 : NT-EV-20  4  FF-ISITE  FIRE  NT  4  -LPl  A  CALL  THE  FIRE  FUNCTION 

C373  EV3: NT-EV-30  4  INT-ISITE  INTER  NT  4  -LPl  A  CALL  THE  INTER  FUNCTION 

C383  EV4: NT-EV— 40  4  AM- I SITE  ASSES  NT  A  CALL  THE  ASSES  FUNCTION 

C39  3  A 

C 40 3  A - RESULT  OF  INTERCEPT - - 

C  4 1 3  - ( TGT C NT 3 .0) /TT1 

C 423  A  -IF  ALL  A/C  OF  FORMATION  ARE  KILLED  THE  TGT  GO  OUT  OF  SIMULATION' 
C433  - ( ( (p DSC CSSs ; (NT-1 )  3  )  -  < I-ASES) ) JO) /TT2  A  NO  DETECTION  STOP  PROCCESS 

C  44  3  -LPl 

C433  DSC CSS; ( ( v ( ( A DSC CSSs : (NT- l ) 3 ) -< I-ASES) ) > -I ) ! (NT+1 ) 3-0  4  -LPl 
C463  A  ALL  DETECTION  COLUMN  FOR  THIS  TGT  BECAMES  ZERO 

C473  TTl : — (MCAP2CSS3*0> /LPl  A  NO  REMAIN  MISSILE  SITE  OUT  OF  SIMULATION 
C4Q3  - (TEST I NT CSS; NT 3*0) /LPl 
C 49 3  A  IF  NO  INTERCEPT  NO  KILL  COUNT 
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C303  A - PK  FUNCTION - 

C313  PK-PKM  PKILL  SS.STDCSS: I;NT3  A  CALL  THE  PK  FUNCTION 
CSC!  A  'PK  «  ',*PK  0  A  'R-  ' ,*STDCSSs I:NT3 


C333  RN-RANDOM  ft*  RN  «  ',*RN  ft  CALL  THE  RANDOM  FUNCTION 
C34  3  - (RN’PK ) /RES31 

C353  149  OPOKE  161  0  'KILL  11 '  ONE  TARGET  , (SNT)  0  C6  OPOKE  161 

C563  TGTCNT3-TGTCNT3-1  0  TK -TK , SS , I , NT 
CS7  3  A KILL  SUBSRUCT  ONE  TGT  FROM  FORMATION 
C58  3  MCAP2  C  SS  3  -MCAPC  C  S3 1  - 1 

C593  SKILL  SUBSTUCT  ONE  MISSILE  FROM  MISSILE  STOCK 
C603  TSTEPCNT  j  SS;  3-TSTEP1  CNT;  SSs  3 
1613  — (TGTCNT3SO) /TTC 

C623  -LP1  ft  IF  ALL  A/C  OF  FORMATION  ARE  KILLED  THE  TGT  GO  OUT  OF  SIM. 
C633  TT2:DS2CSS:  (  U  (  (pDSCCSSs ;  (NT+1)  3) -I)  >-I>  s  (NT+1)  3-0  «  -LP1 
C  64  3  RESS 1 :  MCAP2 C  SS 3  -MC APC  C  SS 3 - 1  ft  NO  KILL  SUBSTRUCT  ONE  MISSILE 
C633  TSTEPCNT ;  SS;  3-TSTEP1  CNT s  SS:  3  0  -LP1 
C66  3  ft 

C673  ft - RESULT  OF  ASSESS - 

C6S3  EV5: 

C  69  3  NT-EV-30 

C703  ESI:  A  ’ THE  TGT  ’ ,  <*NT)  ,  '  ENTER  AGAIN  THE  DETECTION  LIST' 

C71 3  A  INITIALIZE  AGAIN  THE  INTERCEPTION  PROCCESS 

C723  STGTCSS: NT  3 -NT  0  ALLTGTCNT3-NT  0  CAPCC SS 3-CAPCCSS 3 + 1 

C733  A 

C 74 3  ft - DETECTION  EVENT - 

C7S3  L4:-L8  AO  'NO  CURRENT  EVENTS' 

C763  L7:  A  END  OF  EVENTS' 

C77  3  LS: Z-l-DSCCSS;  I;  3 

C783  -  (CAP2CSS3  >0)  /L3  ft* — CAP2ACITY  CONTITON 

C79 3  -END  AO  'THE  SI TEC  IS  FULL- ' , * CAP2CSS 3 
C803  L3:MT— /Z 

C813  -( (MT-O) , (MT-l) , (MT>2> ) /L0,L1,L2 

CSC 3  LO: -END  A  NO  TARGET  IN  TARGET  DETECTION  LIST' 

C833  A - CASE  THAT  RADAR  SITE  CAN  SEE  ONLY  ONE  TGT - 1 

C843  L1:NT-ZU  AO  'TGT  TO  BE  ID1«',»NT 
C833  —  ( Z CNT 3P0 >  /D1 

C863  -END  A  0  'THE  TGT- ' ,  (SNT)  ,  '  IS  NOT  IN  DETECTION  LIST' 

C873  D1 :- (STGTCSS; NT3P0) /AL1 

C883  -END  AO  'THE  TGT  ALLREADY  ID1  =',*NT  0  -END 

C893  A  - COMMAND  AND  CONTROL  MODEL - 

C903  ALls-(CMD-O) /L3 
C913  -(ALLTGTCNT3K0) /L3 

C9C3  — (TGT CNT 3  >AUTO) /L5  A  NUMBER  OF  TGT  GREATER  THAN  ONE 

C933  -END  A  THE  TGT '  ,  (SNT) , '  IS  ALLREADY  ENGAGED  BY  AN  OTHER  SITE' 

C943  L5: CAPCCSS3-CAPCCSS3-1  0  STGTCSS:NT3-0  0  ALLTGTCNT3-0 
C953  ID-ISITE  IDIF  NT  0  -END  aCALL  THE  IDIF  FUNCTION 

C 96 3  A - CASE  THAT  RADAR  SITE  CAN  SEE  TWO  OR  MORE  TGT - 

C973  LC:PRI0R-1*PRCI« 3 

C98  3  LPC :  NT—  1 1  PR  I  OR  0  PRIOR-UPRIOR  • 

C99  3  — (NT— 0) /END 

C 1003  - (ZCNT3K0) /DC 

C 1 0 1 3  -LPC  A  'THE  TGT— ' , (SNT) ,  IS  NOT  IN  DETECTION  LIST'  0  -LPC 
C 1023  DC: - (STGTCSS: NT3*0) /ALC 

C 1 03 3  -LPC  A  'TARGET  ALLREADY  ID'  0  'THE  TGT  VECTOR-' , S STGTCSS; NT 3 

C  1 04  3  A - COMMAND  AND  CONTROL  MODEL - 

Cl  OS 3  ALC: - (CMD-O) /L6 
C  106  3  -(ALLTGTCNT3P0) /L6 

C  1073  - < TGT CNT 3 .AUTO) /L6  A  IF  THE  NUMBER  OF  TGT  GREATER  THAN  VAR.  AUTO 

C 1083  A  THEN  THIS  TGT  IS  ENGANGED  BY  MORE  THAN  ONE  RADAR  SITE 

Cl 09 3  -LPC  A  THE  TGT', (SNT),'  IS  ALLREADY  ENGAGED  BY  AN  OTHER  SITE' 

C 1 103  L6: CAPCC33  3-CAPCCSS3-1  0  STGTCSS: NT3-0  0  ALLTGT CNT 3  — 0 
C  1 1 1  3  A 0  'THE  TGT  VECT  = ' ,  <S STGTCSS: NT 3  > 

CUC3  ID-ISITE  IDIF  NT  AO  'THE  IDC  TGT— '  ,  SNT 
Cl  133  -LPC 


7  1D-I5ITE  IDIF  ft 5  CHECK ;  1 1 :  RN s  RD ;  I DD 


cn  « 

cci  ft - purpose - 

coi  ft 

C41  ft  THIS  FUNCTION  IS  USED  BY  THE  SITE  FUNCTION  TO  PERFORM  THE 
C31  ft  ID INTIF ICATIQN  OF  THE  TARGETS. 

Col  ft 

C71  ft -  INPUTS - 

cai  ft 

C91  ft  l.  THE  LEFT  ARGUMENT  IS  A  VECTOR  WHOSE  FIRST  ELEMENT  IS  THE 
CIO!  ft  NUMBER  OF  RADAR  SITE  AND  THE  SECOND  ONE  IS  THE  TIME  STEP. 

C113  ft  C.  THE  RIGHT  ARGUMENT  IS  THE  NUMBER  OF  TARGET. 

CIC1  ft  3.  THIS  FUNCTION  USES  SOME  OTHER  GLOBAL  VARIABLES. 

C101  ft 

C 141  ft - OUTPUT - 

C  131  ft 

C 161  ft  THE  OUTPUT  IS  THE  FUTURE  I DINT IF I CAT ION  EVENT 
C 171  ft 

C  181  ft - MAIN  PROGRAM - ; - 

C 191  ft 


CCOl  I*ISITEC11  »  RDMSITECC1  9  IDD-TSTEPCR:  RD;  1 1 
CC11  ft  'TGT=*',sR  9  'IDD*',*IDD 

CCC1  -  ( MCAPC C RD 1 >0 ) / TT 1  ft  NO  MISSILE  NO  FUTURE  EVENT 
CCOl  -TTC 

CC41  TTt:  ft 'TARGET  TO  BE  TD=  ',*R 

CCS1  -< IDD< i (A (DSCCRD; sR+11) >-I) ) /TTO  ftSIMULATION  OF  ID  STOP  WHEN  THE 
CC61  ft  REMAINING  TIME  STEPS  ARE  LESS  THAN  THE  NUMBER  THAT  THE  ID  PROCESS 
CC71  ft NEEDS 

CCS1  ID-0  o  DS2CRD;  Ult  ( (p  (DS2CRD;  s  <R*1)  3)  >-t>  > +I  s  Rfl  l«-0  9  -END 
C291  TTC: CHECK-*/ < ,DS2CRD:  sR-ll)  CI+\.IDD1  ft  CHECK  IF  THERE  IS  DETECTION 
COOl  ft  'CHECK  OF  ID“ ' ,* CHECK 
COU  -*( CHECKS 0 ) /L6 

C021  ft  NO  DETECTION  ...  NO  IDIFIC. * 

COOl  DSCCRDs  <1-1 1-M.  <IDD-1>  jR-ll-O 

C041  TTC :  STGT C RD ;  R 1  -R  9  ALLTGTCftl-R 

COS1  ft  'TARGET  ' , <*R) , '  AGAIN  FOR  DETECTION' 

C061  ID-0  9  CAPCCRD1-CAPCCRD1-1  9  -*END 

C071  ft - TEST  FOR  IDIDIFICATIQN - * 

COS1  L6i  •*  < IFF*0) /LI  ft  IFF  NO  IFF  IS  APLLIED  ALL  TARGET  ARE  ENEMIES 
C091  RN-RANDOM  A9  RN=',*RN 
C 40 1  -<  ( RN > I FFWGRK )  ,  (RN<  IFFWOR'K)  )  /L0,L4 

C4I3  L0:-<IFFF0F-O> /LI  9  ft  IF  IFFFOF  *  0  ALL  TARGETS  ARE  ENEMIES 
C4C3  ft  IF  IFFFOF  =  1  ALL  TARGETS  ARE  FRIENDS 

C 40 1  10-0  9  STGTCRD-.Rl-R  9  ALLTGTCR1-R  9  CAP2CRD1-CAPCCRD1-1  9  -END 

C 44 1  ft4  'TARGET' ,  (SR)  ,  ENTER  AGAIN'  9  -END 

C431  L4: RN-RANDOM  ft 9  RN»',»RN 

C4&1  - ( <RN> IFFCQR) , <RNS IFFCOR) ) /L7,LS 

C471  ft  IF  IFF  IS  WRONG  THE  FRIEND  TARGET  IS  ENEMY,  AND  THE  ENEMY  TARGET 
C 481  ft  IS  FRIEND, IF  THE  IFF  IS  CORRECT  THE  OPPOSITE. 

C 49 1  L7:-<R«C!/L9  ft  NUMBER  C  TARGET  IS  FRIEND  BY  THE  SENARIO 
C 30 1  ft  'TARGET  ' ,  <sR)  ,  '  ENTER  AGAIN' 

coil  ID-0  9  STGTCRD: Rl-R  9  ALLTGTCR1-R  9  CAPCCRD 1-CAPCCRD 1+ 1  9  -END 
C3C1  L9: TFERRCR-TFERROR-1  9  -LI  ft  NUMBER  A  FRIEND  TARGET  IS  ENGAGED 
C501  L3:-<Rj«C)  /LI  ft  NUMBER  C  TARGET  IS  FRIEND  BY  THE  SENARIO 
CS4 1  ft  'TARGET  '  ,  (*R>  ,  '  ENTER  AGAIN' 

C 3S 1  ID-0  >  STGTCRD; P 1—R  9  ALLTGTCR1-R  9  CAF'CC RD 1 -CAPCC RD1 -1  9  -END 
C06l  LI: II-' .ASCCRD;  < I-IDD) s 1)  1.0  ftFINDS  WHICH  ELEMENT  ASC  MATRIX  IS  ZERO 
C37I  ASC CRD: <1*1 DD) : III —2 <  JPD  > , ' I ' , <*R)  ft  GENERATE  EVENTS < t*IDIFIC) 

C361  ID-ASCCRD: I-IDD;  1 

CS’l  ft  The  EVENTS  ROW  FOP  IDIT  IS  =  ' ,s ID-ASCCRD; I+IDD; 1 
C601  END: 
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7  TT-ISITE  TRACK  T;  CHECK;  IIRDsTRl 


m  a 

C2]  A - PURPOSE - 

CO]  A 

C4]  hi  THIS  FUNCTION  IS  USED  BY  THE  SITE  FUNCTION  TO  PERFORM  THE 
C53  *  TRACK  TARGET  EVENT. 

C6]  * 

C7]  A -  INPUT - 

C9]  A 

C9]  A  1.  THE  LEFT  ARGUMENT  IS  A  VECTOR  WHOSE  FIRST  ELEMENT  IS  THE 
CIO]  A  RADAR  SITE  NUMBER  AND  THE  SECOND  ONE  IS  THE  TIME  STEP. 

CU]  A  2.  THE  RIGHT  ARGUMENT  IS  THE  NUMBER  OF  TARGET. 

C 12]  A  3.  THIS  FUNCTION  USES  SOME  OTHER  GLOBAL  VARIABLES. 

CIO]  A 

C  1 4  ]  A - OUTFUT - 

CIS]  A 

C 16]  A  THE  OUTPUT  IS  THE  FUTURE  TRACK  EVENT. 

C  17]  A 

C  18]  A - MAIN  PROGRAM - 


C  19]  A 

COO]  I-I3ITEC 1 ]  4  RD-ISITEC2]  4  TR1-TSTEPCT; RD;  2]  4  A  ' TR= ’ ,STR1  0  ’TGT=’,5T 
C21 ]  - (MCAP2CRD] >0) /L2 

C221  -L3 

C27]  A - ECM  MODEL - * 

C24 ]  L2:-(ECMCT;2]=0> /Lii  A  ECM  DOES  NOT  EXIST 
C2S]  RN-RANDOM 

C2s:  - (RN, ECMCT; 51 > /L12  A  ECM  DOES  NOT  WORK 

CC7]  PK>ECM  PKILL  T.STDCRD; I;T] 

C28]  RN-RANDOM 

CC9]  -<RN<PK)/L33  ABREAK  LOCK  QN 

C30]  TSTEPCT ; RDs  0  3] -TSTEPCT ; RD;  2  31+T (TSTEPCT; RD;  2  3 ] )  <  < PK ) 

C31  ]  TR1-TSTEPCT ;  RD ;  21 

C321  -LI  A4  'NO  LOCK  BREAK  NEW  TSTEP= '  4  ’PK* ' , *PK  *  TSTEP  4  -LI 
C33]  L12S-L1  4  A  ECM  DOES  NOT  WORK'  4  • ECM= ’ , SECMCT; 5]  4  RN=’,5RN  4  -LI 
C34]  LUs-Lt  A 4  ' ECM  DOES  NOT  EXIST’  4  -LI 
C3S1  L 1  j  —  ( TR  l  ’  (  (DS2CRD;  ;T-11)  )-I>  ) /LS 

C361  TT-0  4  DS2CRD;  ti.lt  (  <p  (DS2CRD-.  s  (T+l)  ])  >  -  I  >  >-I;T*l]-0  4  -END 
C371  L5:  CHECK -'</  (, DSC CRD:  ;T-l])  CI-1.TR1] 

C33]  A  ’TEST  »  ’ , STE3T  4  CHECK  FOR  TRACK* ’ , 5 CHECK 
C39  ]  -< CHECK* 0 > /L4 

C40]  DSC  CRD;  ( 1-1 )  *>.  (TR1+1 )  ;  T-l  ]-0  4  -L3 
C 4 1 3  a’TGT  =  ’  ,  (5T)  ,  ’  AGAIN  FOR  DETECTION’ 

C42]  L33s  ’BREAK  LOCK  ON’  4  ’PK=’,SPK  4  ’RN=’,5RN  4  DS2CRD; I ; (T+l ) ]-0  4  -l3 

C401  L3 : TT-0  4  STGTCRD; T]-T  4  ALLTGTCT]— T  A4  'TARGET  =  ' , (5 TARGET) 

C  44  ]  CAP2CRD 1-CAP2CR0  ]  + 1  4  -END 

C  4S  ]  L4:  1 1- (AS2CRD;  I-TR1 ;  ]  )  i,0  AIT  FINDS  WHICH  ELEMENT  OS  AS2  MATRIX  IS  ZERO 
C46]  AS2CRD; (I-TRl) ; II]-4 (SRD) , '2' , <*T>  aGENERATE  EVENTS  (2= TRACK) 

C  47 ]  TT-AS2CRD;  I-TRl ;  ] 

C 43 ]  a  the  EVENTS  ROW  FOR  TRACK  IS  =  ’, STT-AS2CRD; I-TRl ;  ] 

C 49 ]  END; 
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7  FR-ISITE  FIRE  F; CHECKRD; RD; FT 
Cl]  A 


•  V  ’ 


v.v 

y.y 


CC]  A - PURPOSE - 

CO]  A 

C 4 ]  A  THIS  FUNCTION  IS  USED  BY  THE  SITE  PROGRAM  TO  PERFORM  THE 
CS1  A  FIRE  EVENT. 

C  6  ]  A 

C  7 1  A -  INPUT - 

C33  A 

C 9 ]  A  1.  THE  LEFT  ARGUMENT  IS  A  VECTOR  WHOSE  FIRST  ELEMENT  IS  THE 
CIO]  A  RADAR  SITE  NUMBER  AND  THE  SECQNT  ONE  IS  THE  TIME  STEP. 

CU]  A  2.  THE  RIGHT  ARGUMENT  IS  THE  NUMBER  OF  TARGET. 

CIC]  A  3.  THIS  FUNCTION  USES  SOME  OTHER  GLOBAL  VARIABLES 
C  13]  A 


162 


ifivt  w  .wvr 


"V"  I  If*1"  WK' 


■  ’.t  w  v  *  v  wv 


T«*«  W  ^  I* IWTfl  ||n  w»-  ■  f4  1*  V 


C  143  - - OUTPUT - ; - - 

C153  4 

C 163  4  THE  OUTPUT  IS  THE  FUTURE  FIRE  EVENT. 

C  173  4 

C  183  4 - MAIN  PROGRAM - + 

C  193  4 

COO  3  I  - 1  SITE  C  1  3  0  RD-ISITEC03  0  FT-TSTEPCF; RD; 33 
C013  — (MCAPOCRD 3  >0)  /L3 

C023  -L4 

C033  L3:  —  (FT< (  if (DS2CRD; ;F+13> )  —  I > ) /L5 

C043  FR-0  0  DS2CRD;  <1.  (  <p  (DS2CRD;  ;  (F-l )  3)  >  —  I)  >-*-!;  F-+-1  3—0  0  -END 
COS  3  L3«  CHECK  «•*/.  (  .DS2CRD;  ;F-13>  CI-S.FT3  44  ’TEST  FOR  FIRE-  '  ,  5CHECK 
C063  - C CHECK# 0) /LI 

C073  DS2CRD;  ( I  - 1 )  +’.  (FT+1 )  ; F+ 1 3-0 

COS  3  L4:FR-0  4  STGTCRD:F3-F  4  ALLTGTCF3-F  4  TSTEPCF;  RD;  3-T3TEP1  CF;  RDs  3 
C 09 3  4  'TARGET  ',(|F),'  ENTER  AGAIN  TO  DETECTION  LIST' 

C303  CAPO CRD 3 -CAPO CRD 3+ l  4  -END 

C313  Lis II-(AS2CRD; I+FT; 3) vO  4FINDS  WHICH  ELEMENT  OF  AS2  MATRIX  IS  ZERO 
C323  AS2CRD: (I-FT) ; 113-4 (*RD> , '3' , («F>  4 GENERATE  EVENTS (3-FIRE) 

C333  FR-ASOCRDs I -FT ; 3 

C343  LO:  4  ’THE  EVENTS  ROW  FOR  FIRE  IS  ■  ' , SFR-AS2CRD; I+FT; 3 
C353  END; 


<7  INT-ISITE  INTER  IN;  CHI ;  CH;  NNN;  I  I ;  RD 
C  1  3  4 

C03  4 - PURPOSE - 

C33  4 

C 43  4  THIS  FUNCTION  IS  USED  BY  THE  SITE  PROGRAM  TO  PERFORM  THE 

C53  4  INTERCEPT  OF  THE  TARGET. 

C63  4 

C73  4 -  INPUT - 

CS3  4 

C93  4  1.  THE  LEFT  ARGUMENT  IS  A  VECTOR  WHOSE  FIRST  ELEMENT  IS  THE 

C 103  4  RADAR  SITE  NUMBER  AND  THE  SECOND  ONE  IS  THE  TIME  STEP. 

CU  3  4  0.  THE  RIGHT  ARGUMENT  IS  THE  NUMBER  OF  TARGET. 

C 103  4  3.  THIS  FUNCTION  USES  SOME  OTHER  GLOBAL  VARIABLES. 

C  133  4 

C  1 4  3  4 - OUTPUT - 

C1S3  4 

C 163  4  1.  THE  TIME  THAT  THE  INTERCEPTION  TAKES  TO  BE  DONE. 

C 173  4  0.  THE  DISTANCE  THAT  THE  TARGET  FLY  UNTIL  GETS  INTERCEPTED. 

C 183  4  3.  THE  DISTANCE  THAT  THE  INTERCEPTOR  FLY  UNTIL  THE  INTERCE 


C 193  4  PTIQN  IS  OQNE. 

COO 3  4 

C01 3  4 - MAIN  PROGRAM 


COO 3  I-ISITEC 1 3  0  RD-ISITEC03 
C033  - (MCAP0CRD3 )  0) /L6 

C043  -LI 

COS 3  L6:  4  NN1-IS0CRD:;  3  6  93 

C06 3  NNN-STSCRD; I:  IN3  40  'NUMB.  OF  TIME  STEPS- ',«NNN 
C073  - ( NNN< ( ( ASEFCSS; : 1 3 ) -I ) ) /L4 

COS 3  INT-0  0  DSOCSS; (I+U ( (pSEFCSS; ; 1 3 ) — I ) ) ) ; IN-13-0 
C  09  3  -END 

C 30 3  L4; CH-SEFCSS; I; IN3  40  'CHECK  FOR  INTERCEPTION- ' ,5CH 
C31 3  CH  1  —  •</  (  ,  3EF  CSS;  ;  I N  3  )  C  I  -  vNNN  3  40  CH1-',*CH1 
C303  - < (CH-O) , (CHl-O) ) /LI ,L0 

C333  -L3 

C 34 3  Lis INT-0  4  0  ’NO  TRACK  ...NO  INTERCPTN ’ 

C3S3  STGTCRDs IN3-IN  0  ALLTGTC IN3-IN  0  TSTEPC IN; RD; 3-TSTEP 1 C IN; RD: 3 
C363  40  ' TGT  *  '  ,  (S IN)  ,  '  ENTER  AGAIN  THE  DETECTION' 

C37  3  CAPO C RD  3 -CAPOC RD 3-1  0  -END 

CCS 3  LO: 'LAUNCHING  WAS  SUCCESSFUL  BUT  ...NO  INTERCEPTION  DUE  TO  MASKING' 
C39 3  TE3TINTCSS:  IN3-0  0  MCAPOCRD 3 -MCAPOCRD 3-1  0  -L5 
C  40  3  L3: TEST I NT CSS; IN3-1 
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C413  L3: II+(AS2CRD: I+NNNs 3U0  4 FINDS  WHICH  ELEMENT  OF  AS2  MATRIX  IS  ZERO 
C423  AS2CRD:  (I+NNN)  :  III--*  (5RD)  ,  '4 ' ,  (*  IN) 

C433  INT+AS2CRD: I+NNN;  3 

C 44 1  a  'THE  EVENTS  ROW  FOR  INTERCFT  IS  =  ,iINT 
C433  END: 

9 


9  AM+ISITE  ASSES  AS; CHECK; 1 1 ; RD; ASES 


CU  A 

C23  4 - PURPOSE - 

C33  A 

C 4 3  h  THIS  FUNCTION  IS  CALLED  BY  THE  SITE  PROGRAM  TO  PERFORM 
C53  A  THE  TARGET  DAMAGE  ASSESSMENT  EVENT. 

C63  A 

C73  A -  INPUT - 

CS3  A 

C 93  A  1.  THE  LEFT  ARGUMENT  IS  A  VECTOR  WHOSE  FIRST  ELEMENT  IS  THE 
C103  A  RADAR  SITE  NUMBER  AND  THE  SECOND  ONE  IS  THE  TIME  STEP. 

Cl  13  A  2.  THE  RIGHT  ARGUMENT  IS  THE  NUMBER  OF  TARGET. 

C 123  A  3.  THIS  FUNCTION  USES  SOME  OTHER  GLOBAL  VARIABLES. 

C  133  A 

C 143  A - OUTPUT - 

C  133  A 

C 163  A  THE  OUTPUT  IS  THE  FUTURE  TARGET  DAMAGE  ASSESSMENT. 

C  173  A 

C 183  A - MAIN  PROGRAM - 


C 193  I+ISITEC l 3  4  RD+ISITEC23  0  ASES+TSTEPC AS; RD; 43 

C203  A  IF  THE  RADAR  DOES  NOT  HAVE  TIME  TO  ASSES  THE  TGT  (THE  TGT 

C21 3  A  ENTER  TO  EACH  OWN  TERRAIN) 

C223  +<ASES<  < (p (DS2CRD; ; AS+13) )-I) ) /L3 

C233  AM+O  0  DS2CRD;  U ( ((XDS2CRD; ? <AS+t> 3) >-I) )*I;AS*l3+0 
C243  AS3+ASS+1  0  +END 

C233  a  TEST  IF  THERE  IS  DETECTON  DURING  THE  ASSES  PROCCESS 
C263  L3:CHECK+x/TEST+( ,DS2CRDs  sAS+13) CI+VASES3 
C273  AO  'TEST  -  ',5 TEST  0  ' CHECK  FOR  ASSE5»  CHECK 
COS 3  +< CHECK* 0) /LI 

C.293  AM+O  0  ASS+ASS+ 1  0  STGTCRDs  AS3+AS  0  ALLTGT C  AS 3  +AS 
C303  0  CAP2CR03-CAP2CRD3+1  0  +END 

C313  A  CREATE  ESSES  EVENTS 

C323  Lis II+(AS2CRD; I+ASESs 3U0  FINDS  WHICH  ELEMENT  OF  AS2  MATRIX  IS  ZERO 
C333  AS2CF.D:  (I+A3E3)  s  II3+*  (SRD)  ,  '3 '  ,  (*A3) 

CC43  AGENERATE  EVENTS  (5=TGT  ASSESS) 

C333  AM+AS2CRD: I+ASESs 3 

C36 3  A 'THE  EVENTS  ROW  FOR  ASSESSEMENT  IS  =  ',*AM 
C373  END: 
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7  Pk-PKM  Pk ILL  R;A;BsCsXl:X2;Yls Y2;P 
C  1  3  A 

C23  A  - PURPOSE - 

C33  A 

C 4 3  A  THIS  FUNCTION  IS  USED  TO  CALCULATE: 


C33  A  1 . THE  PROBABILITY  THAT  A  MISSILE  SITE  KILLS  A  TARGET  BY  USING 
C63  A  THE  APPROPRIATE  FUNCTION  (DATA  IN  PKM  MATRIX) 

C73  A  2. THE  EFFECTIVENESS  OF  ECM  BY  USING  THE  APPROPRIATE  FUNCTION 
C83  A  (DATA  IN  ECM  MATRIX) 

C93  A 

C 10  3  A - INPUTS  - 

CU1  4  !•  THE  FIRST  ARGUMENT  IS  THE  PK  MATRIX  OR  THE  ECM  MATRIX 
C123  A  2.  THE  SECOND  ARGUMENT  IS  A  TWO  ELEMENTS  VECTOR  (SITE  AND 
C  133  A  DISTANCE 

C  1 4  3  A 
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CIS]  A - OUTPUT - 

C  1 6  i  A 

C 173  4  1 . THE  PROBABILITY  THAT  A  MISSILE  SITE  KILLS  A  TARGET 

CIS]  A  0. THE  EFFECTIVENESS  OF  ECU 
C19]  hi 

COO  I  4 - MAIN  PROGRAM - 

CC1I  A^PKMCRC l I ; 1 3 
CCC]  B-PKMCRC1];C] 

CC3]  C-PKMCRCl]s3] 

CC4]  P‘PKMCRCl]s4] 

CCS I  ~<RCC]>A>/L1 

CC6]  Xl*Yl«rO  4  XC*A  4  Y2*P  4  -L3 
CC7]  L1:-<RCC]>B> /LC 

CCS]  XI “A  4  Y1*P  4  XC«-B  4  YC*P  4  -L3 

CC9]  LC:  ■*  (RCC1 >C)  /L4 

C30]  Xl-B  4  Yl^P  4  XC-C  4  YC*-0 

C311  A - THE  USED  FUNCTION - 

C3C1  L3«PK* < ( CRC21-X1)x(Y2-Y1> >-s-(X2-X1> >+Yl  4  -H_5 
C331  L4tPK>0 
C34 1  LS: 
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7  RNtRANDOM 


C  I  1  4 

CC1  4 - ' - PURPOSE - 

C31  rt 

C4]  A  THIS  FUNCTION  IS  USED  BY  OTHER  PROGRAMS  TO  GENERATE  A  RUNDOM 
C51  4  NUMBER 

C61  4 

C71  4 -  INPUT - 

C91  4 

C9]  4  NONE 

CIO!  4 

Cll]  4 - OUTPUT - - - 

ClCl  4 

C 131  4  THE  OUTPUT  IS  A  RUNDOM  NUMBER. 

C  1 4  ]  4 

CIS]  4 - MAIN  PROGRAM - 


C  16]  RN«-  (7100)  +100 
7 


This  Appendix  contains  all  data  o-f  the  pre-calculating 


phase  of  the  air  model.  It  also  lists  the  data  of  the 
attrition  rate  of  the  three  aircraft  formations  by  running 
the  dynamic  phase  of  the  air  model. 


ATTRITION  RATE  OF  THE  THREE  AIRCRAFT  FORMATIONS 


1.  NORMAL  CASE 
TARGET 1 

0  1  0  0  0  2  0  2  0  2  0  0  0  0  0  1001030  0  0101  12103  0  0  1  0  0  2  1 
110  10  10  2  10  0 

TARGET2 

100111000001111010001010010100100001010 
l  1  1  0  1  0  1  0  1  1  0 

TARGET' 

1232232134134222233222124100532121  1  1  1  13 
32333412234 

2.  CASE  WITH  DIFFERENT  IFF  PROBABILITIES 
TARGET! 

20  0  1  1010121  122143323001  1000122021  1  1  1  100 
021112110  10 

TARGET2 

1  101  10223232023023241  00031  1  1000002001  10 

4  0  0  1  1  0  1  0  1  0  1 

TARGETS 

2  0  0  0  0  01  1  13  0  2131324441  143121344434  3  24  0  66 

33243123124 


ATTRITION  RATE  OF  THE  THREE  AIRCRAFT  FORMATIONS 


l .  NORMAL  CASE 
TARGET  1 

0100  0  2  0  2  0  2  9  0  0001001  0  3  0  001  0  1  12103  0  0  l  0  0  2  1 
110  10  10  2  10  0 

TARGET2 

1  0  0  1  1  1  0  0  0  0  0  1  1  1  1  0  1  0  0  O  1  0  1  0  0  1  0  1  0  0  l  0  0  0  0  1  0  1  0 
1110  10  10  110 

TARGET3 


1341342222332221241  0  0  532121  1  1  1  1  3 
3  4  1  2  2  3  4 


2.  CASE  WITH  DIFFERENT  IFF  PROBABILITIES 


TARGET 1 

2  0  0  11  0  101211221433230011  0  00122  0  211111  0  0 
0  2  1112  110  10 

TARGET2 

1  101  1022323202302324100  0  31  1  1  0  0  0  0  0  2  0  0  1  10 

40  0  110101  0  1 

TARGETS 

2  0  0  0  0  0  1  1  13  0  2131324441  143121344434324  0  66 

3324512  3  124 

3.  CASE  WHEN  COMMAND  AND  CONTROL  MODEL  IS  NOT  APPLIED 

TARGET 1 

1  0  0  0  1  0  0  0  1  0  0  0  0  1  0  1  1  2  0  1  0  1  1  0  0  0  1  0  0  0  1  0  0  0  0  1  0  1  1 

2  0  1  1  0  1  0  0  0  1  0 

TARGETS 

1  2  0  2  1  2  1  1  0  1  3  0  3  1  4  1  0  1  2  2  1  2  1  2  0  2  1  2  1  1  0  1  3  0  3  1  4  1  0 

1  2  2  2  2  1  2  0  2  1  2 


251212542334  0  1  15123512251212542334  0  1  151 
255  0  4251212 


4.  CASE  WHEN  RADAR  OPERATION  IS  DEGRADED 


TARGET l 

1  1  0  0  0  1  0  1  1  1  0  0  1  2  0  0  0  11021010 


*/!  >/  tv'C 


target; 


0  1  0  2  2  0  0  0  1  0  0  0  1  0  0  0  0  1  0  0  O  1  1  -  l 

target; 

oooi  13231221011  1  3  2  ;  o  02212 

5.  CASE  WHEN  ECM  MODEL  IS  NOT  APPLIED 
TARGET 1 

oi  i  i  o  o  o  l.  i  o  t  2  i  o  o  o  i  i  i  i  i  i  10021  i  i  c 

target; 

1301  i  o  i  ;  i  o  i  i  o  o  ;  i  i  i  i  i  o  3  ;  o  ;  o  i  o  i  o 
target; 

553357155444431613  4'  13336557334 

ATRITION  RATE  OF  THE  THREE  AIRCRAFT  FORMATIONS 


1.  FLIGHT  LEVEL  *  500  FEET 

TARGET 1 

1  0  1  2  0  0  1  1  1  1  1  1  0  3  1  1  2  0  1  1  1  ;  0  2  1  0  0  0  1  0  1  1  1  1  1  0  1  1  O 

;  i  i  o  o  ;  o  o  ;  o  i  1  1  ;  o  ;  1  o  o  o 

TARGET* 

0  0  0  2  0  0  0  1  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  1  0  0  0  0  1  0  0  1  0  0  0  0  0  0  0 

0  0  1  0  0  0  0  0  0  0  0  0  0  0  0  1  0  0  0  0 

target; 

1211111111010101  0  00221200011111010  1  0  0  12 

ooioi  i  o  i  o  o  o  ;  ;  t  ;  o  o  o  t  1 

;.  FLIGHT  level  *  1000  FEET 

TARGET! 

0  1  0  '".i  0  2  0  2  0  2  0  0  0  0  0  1  0  0  1  0  3  0  0  0  1  0  1  1  2  1  0  3  0  0  1  0  0  2  1 

110  10  1  0  2  10  0 

target; 

1  0  0  1  1  1  0  0  0  0  0  1  1  1  1  0  1  0  0  0  1  0  1  0  0  1  0  1  0  0  1  0  0  0  0  1  0  1  0 

1  1  10  10  10  1  1  0 

target; 

12322321341342222332221241  0  0  532121  1  1  1  13 
3233  3  412234 


16  8 


H 


,s 


c 


1 

•to 


3.  FLIGHT  LEVEL  *  2000  FEET 


TARGET 1 

1  1  1  0  1  1  0  0  0  0  0  1  1  0  1  1  1  1  0  0  1  1  0  1  0  1  1  0  1  2  0  4  1  0  0  1  1  1  1 

1  1  0  2  0  0  0  1  0  1  1 

TARGET2 

0  1  1  t  1  1  0  00  0  010  0  101  1  0  0  3  0  1  0  0  0  0  0  0  0  0  0  1  0  1  1  0  1  1 

0  2  1  0  1  0  0  0  0  1  0 


12  12  2  3 


TARGET" 

1  1222212222212214451  1  0  324313  0  3 
54220221  122 


4.  FLIGHT  LEVEL  »  5000  FEET 


TARGET 1 

0  0  0  3  0  1  2  1  0  1  0  0  1  0  1  0  0  1  3  0  1  1  2  1  2  1  0  1  1  1  2  2  1  0  1  0  1  1  0 
1  1  1  1  1  0  0  3  2  1  3 

TARGET2 

00100101000100110111  0  100220100  0  020101  0  2 
0  1  1  0  0031  102 


TARGET3 

1  2  1  0  1  5  2  4  1  1  2  2  2  1 
213414  0  2214 


20  0  210422321  12124  3  1  122 


5.  FLIGHT  LEVEL  -  10000  FEET 

TARGET 1 

1  0  0  3  0  2  0  0  0  0  1  0  1  0  1  0  0  1  1  1  2  1  1  0  0  0  0  1  2  2  1  1  0  1  1  2  2  1  0 

1  1  0  2  1  0  0  0  0  0  0  0  2 

TARGET2 

2  0  1  0  1  1  0  0  1  0  0  0  0  0  1  0  1  0  0  0  1  0  0  0  1  1  1  1  1  0  0  1  2  1  1  0  0  1  0 

0  1  0  0  1  1  1  0  1  0  0  0  1 

TARGET3 

0  42121  132  0  322  0  315221  1  121  1312223131  0  1  12  0 
3131  0  11243311 

ORIGINAL  DATA  OF  THE  SPECIFIC  SENARIO 

GRID  SQUARE  ALTITUDES 
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OUTPUT  OF  TERRAIN  MODEL 
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TERRAIN  GRID  SQUARE  AVERAGE  ALTITUDES 
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83 
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3023  3000  3123  8130  3900  9623  9300 

5  1033  1313  1373  2030  2923  3650  4350 

«>630  10373  9373  7300  7223  7223  6773 

DATA  OF  RADAR  SITE  MINIMUM  ALTITUDE  MODEL 
l.  SITE  1 


*  1  '1  0  1  l  0  - 1  - 1 

'1  *1  *1  -1  -1 

■  1  1  o  0  1  0  1  1 

■1  'l  *1  -1  -1 

'  1  1 937  1  1  1  66  0  229 

'l  '1  *1  ’I  'I 

1161  1841  1587  126  301  479  510  ",9"’ 

'1  *1  '1  'I  ‘I 

1145  1824  1567  1373  *10  814  529  302 

*1  *1  *1  *1  *i 

:493  2890  3283  2623  1626  1001  934  329 

*1  *1  *1  *1  *1 

*1  4704  3977  3488  1977  2083  947  «64 

'1  *1  *1  *1  *1 

*1  "1  4946  3088  2392  2480  2691  - 1 
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"1  *1 

'1 

*1 

"1 

"1 

*1 

*1 

*1 

"1 

”1 

*1 

"1  *1 

"l 

*1 

"1 

*1  "1 

*1  "1 

~1 

*  1 

*1 

*1 

*1 

*1 

*1 

'1 

”  1 

*1 

*1  *1 

"1 

*  t 

*  1 

*  1 

"1  *1 

*1  *1 

'1 

”  1 

*1 

*1 

*1 

*1 

*1 

"  1 

172  • 


iy 

i 

B 

V  “  T7F 

v: 

V'"’V.T 

vrvv-'v 

■V-  '  V- v 

''IV  V',v 

7vr'vTvi  "V1  v  -j-  ^ 

3296 

4047 

4870 

■I 

-r 

“l 

•t 

“1 

*1 

326 

377 

291 

267 

299 

379 

'10 

1210 

1774 

2414 

Vs 

3121 

4100 

3323 

"1 

"1 

'1 

"1 

ft  -t 

*1 

*1 

660 

488 

593 

492 

1257 

883 

1525 

1923 

3026 

2560 

m 

4075 

5875 

7351 

't 

*1 

'1 

'1 

m 

*  1 

*1 

1 028 

2939 

2370 

1875 

1467 

3790 

3390 

4149 

3319 

4594 

y 

5791 

7275 

8773 

*1 

*1 

*1 

'1 

y.  'i 

*1 

"1 

3804 

3184 

2640 

2181 

6539 

5602 

5343 

5850 

6963 

5825 

y 

7031 

7581 

9137 

*1 

*1 

*1 

*1 

V*.  1 

1 

~  1 

"  1 

vSi'vb 

3022 

9450 

8290 

7357 

7364 

7746 

8635 

9911 

8026 

9146 

"1 

■I 

'1 

"1 

'l 

&  *l 

"1 

"1 

'1 

3989 

12498 

1 1198 

10209 

9606 

9450 

9760 

1  049? 

1 1600 

12989 

"1 

"1 

*1 

'1 

'1 

'1 

p  4 

.  SITE 

4 

K%’  20737 

20455 

20246 

20162 

20204 

20372 

20663 

21076 

2 1 609 

22257 

~  1 

'  1 

~  1 

KS; 

"1 

"1 

"1 

*1 

'1 

'1 

■I 

19123 

18779 

18562 

18475 

18518 

18692 

18994 

19423 

19973 

20642 

21425 

22317 

'1 

EX. » 

"1 

"1 

*1 

'1 

“  1 

'1 

'1 

Hll  17727 

17168 

16942 

16851 

16896 

17077 

17393 

17838 

1 8409 

19102 

19910 

20829 

21853 

"1 

*1 

'1 

"1 

"1 

'1 

'1 

HC  15998 

13623 

15386 

15290 

15328 

13523 

15838 

16324 

16919 

17638 

18475 

19423 

8119 

vO 

3684 

"1 

"1 

"i 

"1 

~  4 

.  *1 

s»  A  1  4540 

14144 

13894 

13793 

13844 

1 4045 

14392 

14881 

13504 

16254 

17123 

18103 

7505 

K> 

6594 

7142 

"1 

'1 

'  1 

"1 

'1 

W  13153 

12734 

12467 

12360 

12413 

12627 

12997 

13514 

14169 

14954 

15858 

6452 

5630 

S^4 

6114 

6656 

7257 

'1 

*1 

'1 

'1 

fyC  ua4i 

11392 

11106 

10990 

11048 

11278 

11674 

12224 

12918 

13743 

5527 

4792 

5212 

€ 

5690 

5855 

6428 

7061 

'1 

'l 

“1 

KFT  10607 

10122 

9810 

9683 

9747 

9998 

10427 

1 1019 

11737 

4727 

4081 

4437 

4851 

»  » 

.1  -* 

4990 

5498 

6065 

6693 

'1 

*1 

*1 

A  *458 

8927 

8581 

8440 

8511 

8790 

9262 

9904 

4054 

3496 

3788 

4137 

4250 

/.I 

4692 

3193 

5735 

6378 

7061 

"1 

"1 

VV.  3404 

7814 

7422 

7260 

7342 

7659 

8187 

8893 

3038 

3266 

•->  1  s 

3635 

4011 

£>*• 

4446 

4941 

5498 

6116 

2726 

*1 

*1 

jiy  7462 

6794 

6 --.-•8 

6144 

6242 

6616 

7219 

2706 

2869 

2894 

3144 

3453 

3821 

■ 

4250 

1353 

1199 

1496 

1853 

'1 

"1 

Hr  3650 

5891 

5338 

5092 

3217 

5679 

2501 

2598 

2594 

2776 

3017 

3318 

631 

737 

431 

616 

672 

965 

"1 

'1 

ft.  3442 

3167 

4454 

4105 

4287 

2352 

2359 

2415 

2528 

695 

608 

242 

103 

& 

142 

242 

418 

643 

931 

1282 

'1 

3035 

2814 

3188 

2439 

2370 

1148 

784 

275 

176 

145 

126 

114 

[ft 

146 

241 

399 

620 

903 

1252 

*1 

2673 

2431 

2463 

'10 

1714 

1026 

539 

264 

169 

140 

145 

134 

140 

190 

289 

452 

677 

965 

1317 

'1 

J?7  2478 

2435 

2189 

2101 

846 

478 

239 

150 

131 

149 

154 

170 

ijv,‘ 

280 

393 

370 

807 

1 107 

1650 

'1 

KJv1  2324 

2013 

2069 

1813 

1778 

289 

189 

155 

130 

126 

179 

203 

241 

rW 

400 

699 

905 

1173 

1570 

'1 

*1 

5N  1963 

1969 

1532 

1588 

1569 

1505 

169 

158 

167 

1 80 

239 

273 

403 

779 

1228 

1501 

1820 

2575 

"l 

*1 

yC*  1 989 

1967 

1407 

1426 

1416 

1382 

1339 

236 

267 

275 

305 

426 

703 

SR 

1073 

1650 

2301 

2775 

3250 

*1 

'1 

2059 

2016 

1319 

1327 

1323 

1309 

1293 

1292 

325 

319 

368 

768 

1 1 08 

y 

1291 

2100 

3100 

3340 

3677 

'1 

'  l 

V-  3181 

1294 

1291 

1291 

1291 

1292 

1 300 

1323 

336 

389 

711 

1245 

1559 

1938 

Va.*J 

-■  v  75 

4470 

’1 

'  1 

£  =358 

1340 

1324 

1318 

1321 

1333 

1358 

1402 

1470 

551 

1210 

1 650 

2050 

■# 

2823 

4100 

5-.»25 

5896 

*1 

*i 

'i 

!^«l  3593 

1 445 

1419 

1409 

1414 

1434 

1473 

1533 

1621 

1 740 

1926 

1875 

2373 

1668 

4073 

5876 

7350 

*1 

*1 

”1 

*1 

1610 

1376 

1562 

1369 

1396 

1646 

1720 

1824 

230 1 

2850 

2801 

3400 

\y 

3023 

7273 

'1 

'1 

'1 

'1 

'1 

fV  1904 

1836 

1793 

1779 

1787 

1 820 

1878 

1965 

2082 

2675 

4075 

4950 

5826 

c> 

6130 

6876 

*1 

*1 

'1 

'1 

"1 

C“-“  3200 

2124 

2078 

2059 

2069 

2106 

2363 

2650 

3275 

4450 

too  75 

7250 

8025 

K- 

"1 

"1 

'1 

'1 

"  l 

"  l 

'l 

Vi  3557 

2474 

2423 

2402 

2413 

2925 

3650 

4350 

3525 

6875 

9350 

9975 

'  1 

i 

"i 

*i 

■  i 

'1 

'1 

"1 

*1 

s? 

174 

Cv 

« 

•i 


DATA  OF  MOVEMENT  MODEL 


.WVjJ 

(QO 

vQ 

ijfl 


CH 


mi 


1.  TARGET 1 


TIME  STEPS 


2 

3 

4 

5 

6 

7 

8 

9 

10 

1 1 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

23 

24 

— -1 

26 

27 

28 

29 

30 

31 

32 

33 

34 

.-5 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

5 1 

52 

5*_> 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 

111 

112 

113 

114 

115 

1  16 

117 

118 

119 

120 

121 

122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

isa 

159 

160 

161 

162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 

1 80 

181 

182 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 

20 1 

202 

203 

204 

205 

206 

207 

208 

209 

210 

21 1 

212 

"ti  - 

214 

215 

216 

218 

219 

220 

221 

223 

224 

226 

227 

229 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

266 

255 

256 

257 

258 

259 

260 

261 

262 

263 

264 

265 

X  CORD I NATE 


5 

4 

5 

6 

7 

9 

10 

n 

12 

13’ 

14 

16 

17 

18 

19 

20 

21 

—  •-> 

24 

25 

26 

27 

28 

29 

31 

32 

55 

34 

35  ' 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

5-j 

54 

55 

56 

57 

58 

59 

60 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

73 

74 

74 

74 

75 

75 

76 

76 

77 

77 

78 

78 

79 

79 

79 

80 

30 

81 

81 

82 

82 

83 

83 

83 

84 

84 

84 

85 

85 

86 

86 

86 

87 

87 

87 

88 

88 

SB 

89 

89 

89 

90 

90 

90 

91 

91 

91 

92 

92 

92 

93 

93 

93 

94 

94 

94 

95 

0 

0 

•0 

0 

0 

0 

93 

92 

92 

91 

91 

90 

89 

89 

88 

88 

87 

86 

86 

85 

85 

84 

84 

•83 

82 

82 

81 

81 

80 

79 

79 

78 

78 

77 

77 

77 

77 

77 

78 

78 

78 

78 

78 

78 

78 

78 

78 

78 

78 

78 

79 

79 

79 

79 

79 

79 

79 

79 

79 

79 

78 

78 

78 

77 

77 

77 

76 

76 

76 

75 

75 

74 

74 

74 

3 

73 

73 

72 

72 

/  _ 

71 

71 

71 

69 

68 

67 

66 

65 

64 

62 

ol 

60 

59 

58 

57 

55 

54 

53 

51 

50 

48 

47 

46 

45 

44 

43 

41 

40 

39 

38 

_■  / 

36 

34 

33 

32 

31 

30 

29 

26 

25 

24 

23 

22 

20 

1° 

18 

17 

16 

15 

13 

12 

1  1 

10 

Q 

8 

o 

C 

4 

V  CORDINATE 


110 

109 

109 

109 

109 

109 

108 

108 

108 

107 

107 

107 

103 

102 

101 

101 

1 00 

100 

96 

95 

94 

94 

93 

93 

89 

08 

37 

87 

86 

86 

78 

77 

76 

75 

74 

73 

65 

64 

63 

62 

6J 

60 

108 

108 

1 09 

104 

1  04 

103 

97 

97 

96 

90 

°0 

99 

81 

80 

79 

68 

67 

6e> 

55 

5-7 

SPEED 


4  20 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

4  20 

420 

420 

420 

420 

420 

4  20 

420 

420 

420 

420 

420 

4  20 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

480 

480 

480 

480 

480 

480 

480 

480 

480 

480 

480 

4Q0 

480 

480 

480 

480 

0 

0 

0 

4ao 

480 

480 

480 

480 

480 

480 

480 

480 

480 

480 

480 

480 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

TARGE 

-  2 

TIME 

STEPS 

- 

3 

4 

5 

6 

7 

14 

IS 

16 

17 

18 

19 

26 

27 

28 

29 

30 

31 

38 

29 

40 

41 

42 

43 

SO 

51 

52 

54 

55 

62 

64 

65 

66 

67 

74 

75 

76 

77 

78 

79 

36 

97 

88 

99 

90 

91 

98 

99 

1 00 

101 

102 

1 03 

110 

1  1  l 

112 

1 13 

114 

1 15 

122 

123 

124 

125 

126 

127 

124 

125 

136 

137 

138 

139 

146 

147 

148 

149 

150 

151 

158 

159 

160 

161 

162 

163 

170 

171 

172 

173 

174 

175 

182 

183 

184 

185 

186  ' 

187 

X  CORD  I NATE 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

o 

0 

0 

o 

0 

0 

3 

4 

6 

7 

8 

16 

17 

18 

19 

20 

21 

29 

30 

31 

32 

34 

42 

44 

45 

46 

47 

49 

57 

58 

59 

60 

61 

63 

71 

72 

73 

74 

75 

76 

32 

84 

35 

36 

87 

98 

96 

97 

98 

99 

1 00 

101 

108 

109 

110 

1  1  1 

112  , 

113 

119 

120 

121 

121 

122 

123 

129 

130 

131 

132 

133 

134 

1  4(5 

141 

143 

143 

143 

1  44 

150 

151 

152 

153 

154 

155 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

480 

480 

480 

480 

480 

480 

480 

480 

480 

480 

480 

480 

480 

480 

480 

480 

480 

480 

0 

0 

0 

480 

480 

480 

480 

480 

480 

480 

480 

480 

480 

480 

480 

480 

480 

480 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

420 

8 

9 

10 

1  1 

12 

13 

20 

21 

24 

25 

^*i 

33 

34 

35 

36 

37 

44 

45 

46 

47 

48 

49 

56 

57 

58 

59 

60 

61 

68 

69 

70 

71 

72 

73 

80 

81 

82 

83 

34 

85 

92 

93 

94 

95 

96 

97 

104 

105 

106 

107 

108 

109 

1 16 

117 

118 

119 

120 

121 

128 

129 

1 30 

131 

132 

133 

140 

141 

142 

143 

144 

145 

152 

153 

154 

155 

156 

157 

164 

165 

166 

167 

168 

169 

176 

177 

178 

179 

180 

181 

188 

189 

190 

0 

0 

0 

0) 

0 

0 

0 

0 

o 

0 

o 

0 

0 

0 

o 

o 

o 

2 

9 

10 

12 

14 

15 

23 

24 

25 

26 

27 

28 

36 

37 

38 

39 

40 

42 

50 

51 

c - 

_ •  _ 

54 

56 

64 

65 

66 

67 

68 

70 

77 

78 

79 

80 

81 

32 

89 

90 

91 

92 

93 

<54 

102 

1 03 

1 04 

105 

106 

107 

113 

114 

115 

1 16 

117 

1 18 

124 

125 

126 

127 

128 

128 

135 

135 

136 

137 

138 

139 

145 

156 

146 

157 

147 

157 

148 

149 

150 

177 


u  »  m  m  ij  ii  n  i  >  ij 


Y  CQRDINATE 


0 

0 

0 

0 

0 

0 

0 

0 

o 

o 

o 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

o 

0 

o 

o 

81 

31 

ao 

80 

30 

79 

79 

78 

78 

78 

77 

77 

76 

76 

75 

75 

75 

74 

74 

74 

73 

7^Z 

72 

71 

71 

70 

70 

70 

70 

69 

69 

69 

69 

69 

69 

69 

69 

o9 

63 

68 

68 

68 

68 

68 

68 

oS 

67 

67 

67 

67 

67 

67 

67 

67 

66 

66 

66 

6© 

66 

66 

65 

65 

64 

64 

63 

£>••> 

62 

61 

61 

60 

60 

59 

59 

58 

58 

57 

57 

56 

55 

55 

54 

54 

53 

53 

52 

ala. 

51 

50 

50 

49 

49 

48 

48 

47 

47 

46 

46 

45 

44 

43 

42 

42 

41 

40 

39 

39 

38 

37 

36 

36 

35 

34 

33 

32 

31 

30 

29 

29 

28 

27 

26 

26 

24 

22 

21 

20 

20 

19 

18 

17 

16 

1 6 

15 

14 

13 

13 

1  —■ 

1 1 

10 

10 

9 

9 

7 

7 

6 

5 

4 

ALTITUDES 


0770 

0295 

0175 

0150 

0158 

10258 

10190 

10267 

10230 


21128 
20928 
22275 
20770 
20295 
20175 
20 1 20 
20225 
10258 
10190 
10267 
10458 


O 

0 

0 

20775 

21128 

21125 

22275 

20770 

20168 

20140 

20130 

20235 

10230 

10190 

10260 

10458 


20775 

20390 

21125 

22275 

20770 

20168 

20140 

20130 

10235 

10230 

10190 

10260 

10458 


20975 
20390 
21125 
2*^^ *75 
20163 
20 1 68 
20 1 40 
20130 
10235 
10230 
10228 
10260 
10458 


20975 
20390 
21125 
21825 
20163 
20168 
20140 
20130 
10235 
10230 
10228 
1 0260 
10458 


21688 
20390 
21125 
21825 
20163 
20168 
20 1 40 
20130 
10218 
10230 
1 0230 
10287 
10369 


21638 

20390 

21813 

21825 

20163 

20168 

20140 
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